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Orientation 


ABSTRACT 

Thin  film  bulk  acoustic  wave  resonators  (FBAR)  using  piezoelectric  AIN,  ZnO  and  GaN  thin  films  have  attracted  extensive  research 
activities  in  the  past  years.  Highly  c-axis  oriented  (normal-plane  orientation)  binary  semiconductor  piezoelectric  thin  films  are  particularly 
investigated  for  resonators  operating  at  the  fundamental  thickness  longitudinal  mode.  Depending  on  the  processing  conditions,  tilted 
polarization  (c-axis  off  the  normal  direction  to  the  substrate  surface)  is  often  found  in  the  as-deposited  piezoelectric  thin  films,  which  leads 
to  the  coexistence  of  thickness  longitudinal  mode  and  shear  mode  for  the  thin  film  resonators. 

In  this  project,  we  intend  to  deposit  AIN  films  with  tilted  polarization  angles  on  Si  substrate  using  a  modified  DC  magnetron  sputtering 
method.  So  far,  AIN  thin  films  with  different  tilted  polarization  angles  have  been  achieved.  High  quality  AIN  thin  films  on  other  suitable 
single  crystal  substrate  (such  as  sapphire)  have  also  been  deposited  and  characterized  using  a  composite  resonator  approach.  A  theoretical 
model  has  been  developed  for  composite  thin  film  resonator,  which  can  be  applied  to  extract  the  electromechanical  material  properties  of 
piezoelectric  thin  films.  With  AIN  film  deposited  and  patterned  on  sapphire  substrate,  a  laser  liftoff  approach  has  been  developed  to  transfer 
the  thin  film  resonators  to  the  receptor  substrate  (silicon)  for  ultimate  on-chip  integration. 

We  have  also  conducted  analytical  study  of  thin  film  bulk  acoustic  wave  resonators  (FBARs)  using  ZnO,  AIN  and  GaN  films  with  a  c-axis 
tilt  angle  (off-normal)  from  0  to  180  degree.  The  tilted  c-axis  orientation  induces  normal  plane  and  in-plane  polarizations,  which  leads  to  the 
coexistence  of  the  longitudinal  mode  and  shear  mode  in  the  resonator.  Interestingly,  it  was  found  that  for  ZnO,  AIN  and  GaN  FBARs, 
special  tilted  angles  were  found  for  pseudo-shear  and  pseudo-longitudinal  modes.  This  result  opens  new  ways  for  FBARs  design  and 
fabrication. 
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1.  Foreword 

Thin  film  bulk  acoustic  wave  resonators  (FBAR)  using  piezoelectric  AIN,  ZnO  and  GaN  thin 
films  have  attracted  extensive  research  activities  in  the  past  years.  Highly  c-axis  oriented 
(normal-plane  orientation)  binary  semiconductor  piezoelectric  thin  films  are  particularly 
investigated  for  resonators  operating  at  the  fundamental  thickness  longitudinal  mode.  Depending 
on  the  processing  conditions,  tilted  polarization  (c-axis  off  the  normal  direction  to  the  substrate 
surface)  is  often  found  in  the  as-deposited  piezoelectric  thin  films,  which  leads  to  the  coexistence 
of  thickness  longitudinal  mode  and  shear  mode  for  the  thin  film  resonators. 

In  this  project,  we  intend  to  deposit  AIN  films  with  tilted  polarization  angles  on  Si  substrate 
using  a  modified  DC  magnetron  sputtering  method.  So  far,  AIN  thin  films  with  different  tilted 
polarization  angles  have  been  achieved.  High  quality  AIN  thin  films  on  other  suitable  single  crystal 
substrate  (such  as  sapphire)  have  also  been  deposited  and  characterized  using  a  composite  resonator 
approach.  A  theoretical  model  has  been  developed  for  composite  thin  film  resonator,  which  can  be 
applied  to  extract  the  electromechanical  material  properties  of  piezoelectric  thin  films.  With  AIN 
film  deposited  and  patterned  on  sapphire  substrate,  a  laser  liftoff  approach  has  been  developed  to 
transfer  the  thin  film  resonators  to  the  receptor  substrate  (silicon)  for  ultimate  on-chip  integration. 

We  have  also  conducted  analytical  study  of  thin  film  bulk  acoustic  wave  resonators  (FBARs) 
using  ZnO,  AIN  and  GaN  films  with  a  c-axis  tilt  angle  (off-normal)  from  0°  to  180°.  The  tilted  c- 
axis  orientation  induces  normal  plane  and  in-plane  polarizations,  which  leads  to  the  coexistence 
of  the  longitudinal  mode  and  shear  mode  in  the  resonator.  The  equation  for  predicting  electric 
impedance  of  FBARs  was  derived  from  the  basic  piezoelectric  constitutive  equations.  Material 
properties  including  elastic,  dielectric  and  piezoelectric  coefficients,  bulk  wave  properties 
including  acoustic  velocity  and  electromechanical  coupling  coefficient,  and  impedance  of  FBARs 
were  calculated  and  showed  strong  dependence  on  the  tilt  angle.  Interestingly,  it  was  found  that 
for  ZnO  FBAR,  pure  thickness  longitudinal  modes  occur  at  0°  and  65.4°,  and  pure  thickness 
shear  modes  occur  at  43°  and  90°.  For  AIN  FBAR,  pure  longitudinal  modes  occur  at  0°  and  67.1°, 
and  pure  shear  modes  occur  at  46. 1°  and  90°  for  AIN.  In  other  words,  pure  thickness  longitudinal 
and  shear  modes  exist  in  ZnO  and  AIN  FBAR  at  specific  tilted  polarization  angles.  In  addition, 
two  peaks  of  shear  mode  electromechanical  coefficient  are  found  at  33.3°  and  90°  for  ZnO,  and 


34.5°  and  90°  for  AIN.  Therefore,  ZnO  and  AIN  films  with  specific  tilt  angles  may  provide 
options  in  the  design  and  fabrication  of  FBARs,  considering  their  strong  shear  resonance  with 
high  electromechanical  coefficients.  The  use  of  dual  mode  FBARs  for  mass  sensor  is  also 
analyzed,  the  calculated  large  resonant  frequency  shift  due  to  mass  loading  shows  their  good 
promising  in  sensor  application  with  high  sensitivity.  The  simulation  results  agreed  well  with  the 
reported  experiment  results,  and  can  be  used  for  design  and  application  of  FBARs.  GaN  thin  film 
with  tilted  polarization  has  also  been  considered  for  FBARs  in  this  work.  Similar  results  have 
been  obtained,  indicating  that  GaN  is  also  a  promising  material  for  FBARs  and  sensor 
applications. 
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1)  Statement  of  the  problem  studied 

Development  of  advanced  wireless  communication  technology  has  led  to  an  explosive 
growth  of  emerging  military  and  consumer  applications  of  RF,  microwave,  and  millimeter  wave 
devices,  circuits  and  systems.  Handheld  communication  systems  with  desirable  frequency  control 
perfonnance  are  highly  desirable  for  US  Anny.  Future  personal  and  ground  communication 
systems  as  well  as  communication  satellites  for  US  Army  necessitate  the  use  of  miniaturized  on- 
chip  thin  film  resonators  and  filters,  featuring  small  size,  low  weight,  high  performance  and  low 
cost,  for  precise  frequency  control  applications.  However,  On-chip  fabrication  and  integration  of 
thin  film  RF  microelectromechanical  resonators  has  been  recognized  as  a  technological  challenge 
in  the  current  microelectronics  and  MEMS  development.  The  currently  used  frequency  control 
devices  such  as  surface  acoustic  wave  (SAW)  and  bulk  acoustic  wave  (BAW)  resonators  or 
filters  are  in  general  off-chip  devices  with  considerably  large  sizes.  Thin  film  bulk  acoustic  wave 
resonators  using  piezoelectric  AIN  thin  films  have  attracted  extensive  research  activities  in  the 
past  few  years  but  with  limited  success.  This  is  largely  due  to  the  processing  difficulties  in  the 
direct  deposition  and  integration  of  high  quality  thin  films  on  Si  substrate.  So  far,  thin  film 
resonators  studied  are  usually  composed  of  materials  fabricated  by  methods  that  do  not  result  in 
desirable  crystal  forms  and  orientations.  Consequently  the  materials  do  not  necessarily  have 
properties  that  closely  resemble  bulk  materials  in  magnitude  or  stability.  The  objectives  of  this 
project  include:  1)  fabrication  of  AIN  thin  film  with  tilted  polarization  (polarization  direction  off 
the  normal  of  the  thin  film  thickness  direction)  to  enable  the  realization  of  dual  mode  thin  film 
bulk  acoustic  wave  resonators  (FBARs);  3)  theoretical  and  experimental  studies  of  dual  mode 
piezoelectric  thin  film  resonators;  3)  development  of  an  on-chip  integration  technology  of  that 
can  ultimately  overcome  the  technological  bottleneck  in  the  miniaturization  and  integration  of 
RF  resonators  and  filters;  4)  explore  the  use  of  FBARs  for  sensor  applications. 

2)  Summary  of  the  most  important  results 

The  novel  dual  mode  thin  film  resonators  we  proposed  to  develop  and  studied  in  the  past  few 
years  in  this  project  are  based  on  the  deposition  of  AIN  films  with  tilted  polarizations,  or  various 
different  film  orientations,  which  may  lead  to  a  great  opportunity  for  the  design  of  high 
performance  frequency  control  devices.  Therefore,  this  project,  targeting  at  the  fundamental 
issues  on  the  deposition  of  high  quality  AIN  thin  films  with  different  orientations,  and  the 
fabrication  and  onchip  integration  of  multimode  thin  film  bulk  acoustic  wave  resonators,  intends 
to  investigate  a  novel  approach  to  facilitate  the  development  of  high  performance  onchip 
frequency  control  devices. 

In  the  first  half  phase  of  this  project,  we  intended  to  deposit  AIN  films  with  tilted  polarization 
angles  on  Si  substrate  using  a  modified  DC  magnetron  sputtering  method.  Through  this  ARO 
funding,  AIN  thin  films  with  different  tilted  polarization  angles  have  been  achieved.  High 
quality  AIN  thin  films  on  other  suitable  single  crystal  substrate  (such  as  sapphire)  have  also  been 
deposited  and  characterized  using  a  composite  resonator  approach.  A  theoretical  model  has  been 
developed  for  composite  thin  film  resonator,  which  can  be  applied  to  extract  the 
electromechanical  material  properties  of  piezoelectric  thin  films.  With  AIN  film  deposited  and 
patterned  on  sapphire  substrate,  a  laser  liftoff  approach  has  been  developed  to  transfer  the  thin 


film  resonators  to  the  receptor  substrate  (silicon)  for  ultimate  on-chip  integration.  More 
specifically,  we  have  perfonned  the  following  research  tasks: 

-  Deposition  of  piezoelectric  AIN  films  with  various  orientations  on  silicon  and  sapphire 
substrate  by  DC  magnetron  sputtering  and  MOCVD. 

-  Use  of  SEM  and  XRD  to  characterize  the  structural  characteristics  of  the  as-deposited 
films. 

-  Characterization  of  nanomechanical  properties  of  the  thin  film  materials 

-  Modeling  of  the  electric  impedance  of  three  layer  and  four  layer  thin  film  composite 
resonators 

-  Development  of  application  of  a  spectrum  method  for  the  modeling  and  characterization 
of  composite  thin  film  resonators 

-  Use  of  Laser  interferometer  technique  to  characterize  the  piezoelectric  property  of  the 
AIN  thin  films 

-  Sacrificial  layer  GaN  was  used  between  the  AIN  thin  film  and  the  sapphire  substrate  to 
assist  excimer  laser  liftoff  and  transfer  of  the  thin  film  devices  from  sapphire  substrate  to 
receptor  substrate  (Si);  and  ultimate  onchip  integration  of  thin  film  resonators  has  been 
successfully  demonstrated 

-  Characterization  of  frequency-temperature  stability  of  thin  film  resonators 

In  the  second  phase  of  this  project,  we  have  conducted  analytical  study  of  thin  film  bulk 
acoustic  wave  resonators  (FBARs)  using  ZnO,  AIN  and  GaN  films  with  a  c-axis  tilt  angle  (off- 
normal)  from  0°  to  180°.  The  tilted  c-axis  orientation  induces  normal  plane  and  in-plane 
polarizations,  which  leads  to  the  coexistence  of  the  longitudinal  mode  and  shear  mode  in  the 
resonator.  The  equation  for  predicting  electric  impedance  of  FBARs  was  derived  from  the  basic 
piezoelectric  constitutive  equations.  Material  properties  including  elastic,  dielectric  and 
piezoelectric  coefficients,  bulk  wave  properties  including  acoustic  velocity  and 
electromechanical  coupling  coefficient,  and  impedance  of  FBARs  were  calculated  and  showed 
strong  dependence  on  the  tilt  angle.  Interestingly,  it  was  found  that  for  ZnO  FBAR,  pure 
thickness  longitudinal  modes  occur  at  0°  and  65.4°,  and  pure  thickness  shear  modes  occur  at  43° 
and  90°.  For  AIN  FBAR,  pure  longitudinal  modes  occur  at  0°  and  67.1°,  and  pure  shear  modes 
occur  at  46. 1°  and  90°  for  AIN.  In  other  words,  pure  thickness  longitudinal  and  shear  modes  exist 
in  ZnO  and  AIN  FBAR  at  specific  tilted  polarization  angles.  In  addition,  two  peaks  of  shear 
mode  electromechanical  coefficient  are  found  at  33.3°  and  90°  for  ZnO,  and  34.5°  and  90°  for 
AIN.  Therefore,  ZnO  and  AIN  films  with  specific  tilt  angles  may  provide  options  in  the  design 
and  fabrication  of  FBARs,  considering  their  strong  shear  resonance  with  high  electromechanical 
coefficients.  The  use  of  dual  mode  FBARs  for  mass  sensor  is  also  analyzed;  the  calculated  large 
resonant  frequency  shift  due  to  mass  loading  shows  their  good  promising  in  sensor  application 
with  high  sensitivity.  The  simulation  results  agreed  well  with  the  reported  experiment  results, 
and  can  be  used  for  design  and  application  of  FBARs.  GaN  thin  film  with  tilted  polarization  has 
also  been  considered  for  FBARs  in  this  work.  Similar  results  have  been  obtained,  indicating  that 
GaN  is  also  a  promising  material  for  FBARs  and  sensor  applications. 


Mass  sensitivity  of  thin  film  bulk  acoustic  resonator  sensors  based  on  polar  c-axis  tilted 
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The  mass  sensitivity  of  thin  film  bulk  acoustic  wave  resonators  (FBARs)  using  ZnO  and  AIN  thin 
films  with  tilted  polar  c-axis  has  been  theoretically  investigated.  The  tilted  c-axis  orientation  induces 
normal  plane  and  in-plane  polarizations,  which  leads  to  the  coexistence  of  thickness  longitudinal  mode 
and  thickness  shear  mode  in  the  resonators.  The  equation  for  predicting  electric  impedance  of  FBARs 
with  a  mass  loading  layer  was  derived  from  the  basic  piezoelectric  constitutive  equations.  The  mass 
sensitivity  of  ZnO  and  AIN  dual  mode  resonators  was  found  by  calculating  the  resonant  frequency  shifts 
of  the  thickness  shear  mode  and  thickness  longitudinal  mode  due  to  thin  film  mass  loading.  In  the 
calculation,  ZnO  and  AIN  thin  film  has  a  c-axis  tilt  angle  from  0°-90°,  2pm  thickness  and  300pm  by 
300pm  electrode  area;  four  different  materials  (Al,  Si02,  Au  and  Pt)  were  considered  as  the  thin  film 
mass  loading.  It  was  found  that  both  longitudinal  resonance  frequency  and  shear  resonance  frequency  for 
different  c-axis  angle  have  a  significant  shift  due  to  the  mass  loading;  the  mass  sensitivities,  defined  as 
Af/(fAm),  of  the  longitudinal  and  shear  mode  for  this  four  mass  loading  materials  are  very  close,  and  do 
not  change  much  with  c-axis  tile  angle  with  a  value  rang  around  -900  cm2/g  for  ZnO  FBARs  and  -1550 
cm2/g  for  AIN  FBARs.  The  results  are  applicable  for  design  and  application  of  ZnO  or  AIN  FBARs. 

KEYWORDS:  tilted  c-axis;  thin  film  bulk  acoustic  resonator;  acoustic  wave  sensor;  mass  sensitivity. 
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I.  INTRODUCTION 


Due  to  its  high  sensitivity,  simple  structure,  and  easy  interconnection  with  electronic 
measurement  systems,  bulk  acoustic  wave  (BAW)  resonators  have  been  studied  for  mummeries  sensing 
applications  in  physical  and  chemical  processes.1  Meanwhile,  with  the  improvement  of  device  fabrication 
and  material  growth  techniques,  now  the  resonator  can  be  made  with  very  small  size.  Especially  thin  film 
bulk  acoustic  wave  resonators  (FBARs)  based  on  ZnO  and  AIN  can  be  achieved  with  thickness  from 
several  micrometers  down  to  tenth  of  micrometers.  It  is  well  known  that  the  resonant  frequency  is 
inversely  proportional  to  the  resonator  thickness;  hence,  the  resonant  frequencies  of  this  kind  resonators 
can  be  very  high,  which  usually  has  a  range  from  hundreds  of  MHz  to  tens  GHz.  Commonly,  a  BAW 
resonator  consists  of  one  piezoelectric  layer  with  two  electrodes  deposited  on  both  sides.  For  mass  sensor 
application,  special  material  sensitive  to  target  is  coated  on  the  resonator  surface.  Surface  change  due  to 
physical/chcmical  adsoiption  and  absorption  will  induce  resonant  frequency  shift,  and  the  absolute 
frequency  shift  is  proportion  to  the  square  of  the  operating  frequency  (resonant  frequency).2’ J  In  other 
words,  the  higher  sensor  sensitivity  can  be  accomplished  by  increasing  the  resonator’s  resonant  frequency. 
Therefore,  adopting  FBAR  structures  for  mass  sensor  application,  can  greatly  improve  the  sensor 
sensitivity  considering  their  high  resonant  frequency.  In  fact,  FBARs  have  been  experimentally 
investigated  as  gas  sensors  for  H2,  CO  and  ethanol,4  chemical  sensors  for  metal  ion  Hg+  detection  in 
liquid  environment,5  and  biosensors  for  DNA  and  protein  molecules,6’  7  which  showed  much  higher 
sensitivity  than  quartz  microbalance  (QCM). 

It  should  be  noticed  that  ZnO  and  AIN  thin  film  used  for  FBAR  fabrication  are  usually  with 
normal-plane  c-axis  orientation,  which  makes  FBARs  operated  in  thickness  longitudinal  mode.  For  sensor 
application  in  liquid,  acoustic  wave  energy  of  longitudinal  mode  is  radiated  into  liquid  through 
compressional  motion,  inducing  significant  damping  of  resonators;  while  the  damping  of  shear  mode 
acoustic  wave  is  not  severe.  Hence,  FBARs  operated  in  shear  mode  is  very  promising  for  liquid 
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application.  Fortunately,  it  has  shown  that  thickness  shear  mode  can  be  excited  besides  longitudinal 
mode  when  AIN  and  ZnO  thin  films  has  a  tilted  c-axis  orientation  ,8 13 

However,  the  theoretical  analysis  of  resonator  based  on  c-axis  tilted  ZnO  and  AIN  is  less.  To  our 
knowledge,  there  has  not  been  a  throughout  study  on  mass  sensitivity  of  dual  mode  FBARs  based  on  c- 
axis  tilted  ZnO  and  AIN.  For  this  purpose,  we  derived  the  electric  impedance  expression  of  c-axis  tilted 
ZnO  and  AIN  FBARs  with  a  mass  loading  layer  by  the  basic  piezoelectric  equations,  and  examine  the 
resonant  frequency  shift  of  longitudinal  and  shear  mode  for  different  mass  loading;  further  the  mass 
sensitivity  of  dual  mode  resonators  will  be  calculated  and  discussed. 

II.  Theory 

A.  Dual  mode  ZnO  and  AIN  FBARs  with  a  mass  loading  layer 

Fig.  1(a)  shows  schematic  of  normal  c-axis  ZnO  or  AIN  FBARs  with  a  mass  loading  layer,  and 
(Xj,  x2,  x3)  is  the  original  material  coordinate;  Fig.  1(b)  shows  the  schematic  of  c-axis  tilted  ZnO  and 

AIN  FBARs  with  a  mass  loading  layer.  Compared  with  piezoelectric  layer  and  mass  loading  layer,  the 
acoustic  impedance  of  electrodes  is  very  small,  hence  for  simplifying  the  problem,  the  two  electrodes  are 
ignored.  A  rectangular  Cartesian  coordinate  system  (  xl ,  x2 ,  x3 )  is  chosen  with  top  electrode  on  x3  =  h 

and  bottom  electrode  on  x3  =  0  .  ZnO  or  AIN  thin  film  has  a  tilt  angle  9  with  x3  direction.  As  shown  in 
Fig.  1  (c),  ( Xj ,  x2 ,  x3 )  can  be  treated  as  the  result  of  rotation  of  ( x, ,  x2 ,  x3 )  about  x2  with  an  angle  0 . 
For  the  characterization  of  c-axis  tilted  ZnO  and  AIN  FBARs,  the  material  properties  (dielectric 
permittivity,  piezoelectric  stress  coefficient,  elastic  stiffness)  or  (s  ,  e ,  c  )  in  (x,,  x2,x3),  need  to  be 
known,  and  they  can  be  computed  through  corresponding  value  ( £  ,  e  ,  c  )  in  original  coordinate 
system(  xn  x2,x3)  with  the  aids  of  matrix  algebra14: 
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where  0  is  the  tilt  angle,  a1  is  the  transpose  of  a  ,  and  M'  is  the  transpose  ofM  . 

B.  Impedance  equation  of  oaxis  tilted  ZnO  and  AIN  FBARs  with  a  mass  loading  layer 

Following  a  similar  procedure15,  the  impedance  equation  for  c-axis  tilted  ZnO  and  AIN  FBARs 


with  a  mass  loading  layer  is  derived.  In  the  coordinate  system  ( x, ,  x2 ,  x3 ),  we  have  the  following 


acoustic  field  equations  and  piezoelectric  constitution  equations14: 
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where  Tp,  Sq,  Di  and  Ek  are  the  components  of  stress,  strain,  electric  displacement  and  electric  field 
intensity,  cpq  are  the  elastic  stiffness  constants  under  constant  electric  field  intensity,  ekp  are  the 
piezoelectric  stress  constants,  sfk  are  the  permittivity  constants  under  constant  strain,  and  u;  is  the 

displacement  in  the  direction  of  xt  (i,k  =  1,  2,  3  and  p,q  =  1,  2,  3,  4,  5,  6 ).  The  resonator  vibration  can 
be  treated  as  one  dimension  problem  considering  its  high  ratio  of  lateral  dimensions  to  thickness  of 
resonator.  Hence,  we  can  assume: 

ST^=dTL  =  o 
cbq  dx2 

El  =  E2  =  0,  E3  0  (6) 
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For  a  sinusoidal  excitation,  the  voltage  V  and  current  1  across  FBAR  are  expressed  by 

v  =  \[e& 

I  =  jcoAD,  (7) 


where  A  is  the  area  of  the  electrodes,  CO  is  the  angle  frequency.  For  traction  forces  at  the  boundary 
(x3  =  0,// )  we  have: 
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z’ml  and  Zfnl  are  the  longitudinal  and  shear  acoustic  impedance  of  mass  loading  layer,  which  can  be 
expressed  by 
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where  v'ml  and  vsml  are  longitudinal  and  shear  acoustic  velocity  of  mass  loading  layer;  and  ysml  are 
longitudinal  and  shear  phase  delay  in  mass  loading  layer;  pml  and  lml  are  the  density  and  thickness  of 
mass  loading  layer;  c33  and  c33  are  the  elastic  constants  of  mass  loading  layer.  Here,  we  assume  the  mass 
loading  layer  is  isotropic. 
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From  (1)  to  (10),  the  impedance  of  FBAR  can  be  solved: 
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Here,  f  is  frequency,  v(L)  and  v(S)  are  the  acoustic  velocities  for  longitudinal  mode  and  shear  mode, 
k\  and  k\  are  defined  as  the  electromechanical  coupling  coefficient  of  longitudinal  and  shear  mode, 
respectively. 

III.  RESULTS  AND  DISCUSSION 

A.  Material  properties  of  oaxis  tilted  ZnO  and  AIN 

Based  on  Eq.  (1)  and  material  property  values  of  ZnO16'17  and  AIN18  in  (x,  ,  x2,x3 )  shown  in  Fig 
1(d)  and  (e),  the  material  properties  of  oaxis  tilted  ZnO  and  AIN  i  n  ( x, ,  x2 ,  x3 )  can  been  calculated.  Fig. 

2  and  3  shows  the  nonzero  components  of  elastic  constants  cE  and  cD  ,  where  cE  is  the  elastic  stiffness 
constant  under  constant  electric  field  intensity,  and  c  ’  is  the  elastic  stiffness  constant  under  constant 
electric  displacement.  For  0°<#<180°,  both  cE  and  cD  are  found  to  be  symmetric  as  6x6  matrixes  with 


b  b  b  b  b  b  b  b  r\  j  U  U  u  U  U  L)  U  U  r\ 

C14  —  C2A  ~  C34  —  C16  —  C26  —  C36  —  C56  —  C45  —  ^  aI1<^  C14  —  C24  —  C34  —  C16  —  C26  —  C36  —  C56  —  C45  —  ^  • 

As  shown  in  Fig.  2  and  3,  the  most  nonzero  elastic  constants  are  strongly  dependant  on  the  tilt  angle. 
However  c22  and  c2n  do  not  change  with  the  tilt  angle,  which  is  due  to  the  fact  that  and  x2  are  the 


same  axis. 


Fig.  4  and  5  plot  the  nonzero  components  of  dielectric  constant  under  constant  strain  (es  )  and 
dielectric  constant  under  constant  stress  ( sT  ),  which  both  keep  symmetric  as  3  x  3  matrixes  when  6 
changes  from  0°to  180°.  As  same  as  dielectric  constants,  s22  and  s22  keep  constant  because  x2  does  not 
change  with  0;  the  other  dielectric  constants  show  great  dependency  on  c-axis  tilt  angle.  For  example,  the 
change  of£n  or  c] ,  with  0  has  an  opposite  trend  with  c33  or  s33  ;  at  6=  90  su  or  sn  has  a  maximum 


value  ,and  e33  or  s33  has  a  minimum  value  ;  while,  at  0=0",  s33  or  s33  reaches  maximum,  and e33  or 


0  „S' 


£3l  reaches  minimum,  which  can  be  explained  by  the  fact  that  x,  and  x3  coincide  with  xl  and  x3 


respectively  when  0  is  0°,  and  x,  becomes  x3  and  x3  coincides  with  negative  x,  when  0  is  90°. 


The  strong  oaxis  tile  angle  dependency  of  material  properties  also  can  be  seen  in  piezoelectric 
strain  constant  d  ( 3x6  matrix)  and  piezoelectric  stress  constant  e  (3  x  6  matrix),  and  their  non-zero 
constants  are  shown  in  Fig.  6  and  7. 

B  Simulation  of  impedance  spectra  of  oaxis  tilted  ZnO  and  AIN  FBARs  with  a  mass  loading  layer 

Using  the  impedance  expression  of  FBAR  and  the  calculated  material  properties  including  elastic, 
dielectric  and  piezoelectric  coefficients,  the  resonator  spectrum  with  mass  loading  can  be  simulated.  The 
parameters  used  in  the  simulation  are  listed  in  Table  1.  The  calculated  impedance  spectra  of  FBARs  with 
different  mass  loading  are  shown  in  Fig.  8  and  9.  For  most  angles,  it  can  be  seen  that  there  are  two 
resonance  peaks.  The  first  peak  is  shear  resonance  with  resonant  frequency  around  0.7  GFIz  for  ZnO 
FBAR  and  1.5  GFIz  for  AIN  FBAR,  and  the  second  one  is  longitudinal  resonance  with  resonant  frequency 
around  1.5GFIz  for  ZnO  FBAR  and  2.7  GFIz  for  AIN  FBAR.  While  For  some  specific  angles,  there  is 
only  one  peak  in  the  spectrum.  This  is  because  the  resonator  could  be  operated  in  a  single  mode  when  a 

tilt  angle  makes  k]  =  0  (pure  shear  mode)  or  k2  =  0  (pure  longitudinal  mode)  according  to  the 
impedance  expression  of  FBAR.  Through  solving  k\=  0  and  k2s  =  0 ,  it  is  found  that  pure  longitudinal 

mode  occurs  at  0°  and  65.4°,  and  pure  shear  mode  occurs  at  43°  and  90°  for  ZnO  FBAR;  pure  longitudinal 
mode  occurs  at  0°  and  67.1°,  and  pure  shear  mode  occurs  at  46.1°  and  90°  for  AIN  FBAR.  As  mentioned 
before,  the  appearance  of  strong  shear  resonance  for  some  certain  tilt  angles  makes  ZnO  or  AIN  FBARs 
desirable  for  sensor  application  in  liquid. 

C.  Mass  sensitivity  of  oaxis  tilted  ZnO  and  AIN  FBARs 

From  Fig.  8  and  9,  we  can  see  that  both  longitudinal  peak  and  shear  peak  shift  from  right  to  left 
due  to  the  mass  loading.  To  have  an  overall  understanding  of  dual  mode  FBAR  as  mass  sensor,  we 
gradually  increase  the  mass  loading  by  increase  the  thickness  of  mass  loading  layer,  and  examine  the 
change  of  the  resonant  frequency  fs .  fs  is  the  frequency  corresponding  to  the  maximum  conductance  of 
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resonator,  which  is  preferred  as  sensor  output.  For  a  FBAR  based  on  AIN  or  ZnO  with  oaxis  tilt  angle  9 , 
we  adopt  fsL  (9)  and  fsS  (9)  to  express  the  longitudinal  and  shear  resonant  frequency,  respectively. 


Fig.  10  and  1 1  shows  the  change  of 


'  and  ' 


with  the  thickness  of  mass  loading  layer 


A/J9)  AfsS{9) 

for  some  specific  c-axis  tilt  angles,  where  the  materials  for  mass  loading  layer  are  Al,  Si02,  Au  and  Pt 
with  thickness  from  0  to  500nm,  and  the  material  parameters  for  simulation  are  shown  in  Table  1.  It  can 

f  ($)  f  ( 9 ) 

be  seen  from  Fig.  10  and  11  that  both  — — — —and -AU - are  linearly  proportional  to  the  thickness  of 

A 4 4(60 

mass  loading  layer  when  the  thickness  is  small.  However,  with  further  increased  thickness,  the  decrease 


of 


fAO)  fAO) 


-and- 
A fAO)  A fsS{9) 


become  fast  for  low  acoustic  impedance  (lower  than  piezoelectric  layer  ZnO  or 


AIN)  mass  loading  material  Al  and  Si02,  and  slowly  for  high  acoustic  impedance  (higher  than 
piezoelectric  layer  ZnO  or  AIN)  mass  loading  material  Au  and  Pt,  which  is  believed  to  be  due  to  the 
mismatch  of  acoustic  impedance  between  mass  loading  layer  and  piezoelectric  layer20  In  addition  we  can 
easily  find  that  the  resonant  frequency  shift  for  Au  and  Pt  is  higher  than  Al  and  Si02,  and  this  can  be  well 
explained  by  that  Au  and  Pt  have  higher  acoustic  loading  because  of  higher  densities. 

Following  the  mass  sensitivity  definition  of  sing  mode  FBARs19,  the  mass  sensitivity  of  dual 
mode  resonator  is  defined  as  below: 


Sl(9) 


A f AO)  If  AO) 

Am 


,  Ss(9 ) 


A fsS(9VLs(9) 

Am 


(12) 


where  SL  (9)  is  the  longitudinal  mode  mass  sensitivity  of  ZnO  or  AIN  FBAR  with  c-axis  tilt  angle  9 ,  and 
Ss{9)  is  the  shear  mode  mass  sensitivity;  AfsL(9)  and  AfsS(9)  are  resonant  frequency  shift  for  mass 
loading  Am  =  pmtm ;  pm  and  tm  are  the  density  and  thickness  of  mass  loading  layer. 

Using  (12),  the  mass  sensitivity  of  unload  FBARs  is  calculated  through  measuring  resonant 
frequency  shift  for  lnm  mass  loading.  Table  11  and  111  list  the  mass  sensitivity  of  FBAR  based  on  2pm 
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ZnO  or  AIN  FBARs  with  c-axis  tilt  angle  in  the  range  of  0°-90°.  Form  Table  11  and  111,  we  can  see  that 
both  SL(0)  and  Ss(8)  are  almost  same  for  different  mass  loading  material,  and  their  minor  difference 
can  be  understood  since  the  piezoelectric  layer  and  mass  loading  layer  is  not  completely  non  lossy  in  our 
simulation.;  we  also  can  see  that  SL  (0)  has  similar  value  with  Ss(0) ,  and  they  do  not  change  much  with 

tilt  angle:  SL(ff)  and  Ss(0)  of  ZnO  FBAR  are  around  -900(cm2/g),  and  SL(0)  and  Ss{6)  of  AIN 
FBAR  are  around  -1550(cm2/g).  The  reason  is  shown  below,  according  to  (1 1),  the  longitudinal  and  shear 
mode  vibration  are  found  to  interact  each  other  through  parameter  a  ;  for  any  9 ,  cf5  is  very  small 


compared  with  cf3  and  c$5  (shown  in  Fig.  2  and  3),  which  induces  small  a  less  than  5°.  Hence  the 
coupling  between  longitudinal  and  shear  mode  is  thought  to  be  weak,  and  the  resonator  can  be 
approximate  to  a  simple  combination  of  two  single  mode.  The  mass  sensitivity  ( Sm  )  of  single  mode 
resonator  is  found  only  dependant  on  the  density  ( p  )  and  thickness  ( d  )  of  the  piezoelectric  layer19. 


S„  - 


p  d 

r  p  p 


(37) 


Hence,  for  ZnO  or  AIN  dual  mode  FBARs,  SL(0)  and  Ss(9)  are  close,  and  do  not  change  much  with  tilt 

angle  although  the  material  properties  are  greatly  affected  by  the  c-axis  tilt  angle. 

It  should  be  pointed  out  that  the  mass  sensitivity  (shown  in  Table  11  and  111)  with  -876  (cm2/g)  to  - 
925  (cm2/g)  for  2pm  ZnO  FBAR  and  -1534  (cm2/g)  to  -1571  (cm2/g)  for  2pm  AIN  FBAR,  are  pretty  high 
compared  with  the  mass  sensitivity  of  a  6M  QCM,  -14(cm2/g)  19.  This  high  mass  sensitivity  with  the 
existence  of  pure  shear  mode  makes  ZnO  and  AIN  FBARs  promising  for  sensor  application  in  liquid. 


IV.  CONCLUSIONS 

In  summary,  mass  sensitivity  of  FBARs  based  on  c-axis  tilted  ZnO  and  AIN  thin  films  have  been 
theoretically  studied.  It  is  found  that  the  material  properties  and  impedance  spectra  of  FBARs  are  greatly 
dependant  on  the  c-axis  orientation.  Pure  shear  resonant  peak  is  found  at  43°  and  90°  for  ZnO,  and  46.1° 
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and  90°  for  AIN.  The  mass  sensitivity  of  longitudinal  and  shear  mode  for  2|im  ZnO  or  AIN  FBAR  is  very 
close,  and  do  not  change  much  with  c-axis  tile  angle  with  a  value  rang  around  -900  cm2/g  for  ZnO  and  - 
1550  cm2/g  for  AIN.  This  high  mass  sensitivity  and  the  existence  of  pure  shear  mode  make  ZnO  and  AIN 
FBARs  desirable  for  sensor  application  in  liquid. 
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Caption  of  Tables  and  Figures: 

TABLE  I.  Materials  properties  and  data  used  in  theoretical  calculation 
TABLE  II.  Mass  sensitivity  of  FBAR  based  on  oaxis  tilted  ZnO 
TABLE  II.  Mass  sensitivity  of  FBAR  based  on  c-axis  tilted  AIN 

FIG.  1.  The  schematic  of  FBARs  based  on  c-axis  tilted  ZnO  or  AIN  thin  film  and  coordinate  systems:  (a) 
FBAR  based  on  ZnO  or  AIN  thin  film  with  normal  c-axis  and  coordinate  system  (Xi,  X2,  X3),  (b)  FBAR 
based  on  ZnO  or  AIN  thin  film  with  tilted  c-axis  and  coordinate  system  (X!  ,  X2 ,  X3  ),  (c)  the  relation  of 
(X1;  X2,  X3)  and  (Xi  ,  X2 ,  X3  ),  (d)  ZnO  material  properties  in  (Xl5  X2,  X3),  and  (e)  AIN  material 
properties  in  (X3,  X2,  X3). 

FIG.  2.  The  elastic  constants  of  ZnO  thin  film  as  a  function  of  tilt  angle:  (a)  cE  ,  (b)  cD 

FIG.  3.  The  elastic  constants  of  AIN  thin  film  as  a  function  of  tilt  angle:  (a )cE  ,  (b)  c" 

FIG.  4.  The  dielectric  constants  of  ZnO  thin  film  as  a  function  of  tilt  angle:  £33  ,  £33 

FIG.  5.  The  dielectric  constants  of  AIN  thin  film  as  a  function  of  tilt  angle:  £33  ,  s33 

FIG.  6.  The  piezoelectric  constants  of  ZnO  thin  film  as  a  function  of  tilt  angle:  (a)  e  ,  (b)  d  . 

FIG.  7.  The  piezoelectric  constants  of  AIN  thin  film  as  a  function  of  tilt  angle:  (a)  e  ,  (b)  d  . 

FIG.  8.  Simulation  of  impedance  spectrum  of  FBAR  based  on  c-axis  tilted  ZnO  thin  film(S:  Shear  mode, 
L:  Longitudinal  Mode). 

FIG.  9.  Simulation  of  impedance  spectrum  of  FBAR  based  on  c-axis  tilted  AIN  thin  film(S:  Shear  mode, 
L:  Longitudinal  Mode). 

FIG.  10.  Resonant  frequency  shift  of  FBAR  based  on  c-axis  tilted  ZnO  film  with  different  mass  loading. 
FIG.  11.  Resonant  frequency  shift  of  FBAR  based  on  c-axis  tilted  AIN  film  with  different  mass  loading. 
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TABLE  I.  Materials  properties  and  data  used  in  theoretical  calculation 


Piezoelectric 

layer 

Density 

(kg/m3) 

Thickness 

(pm) 

Electrode  area 
(pmxpm) 

Mechanical 
quality  factor  of 
longitudinal  mode 

Mechanical 
quality  factor  of 
Shear  mode 

ZnO 

5700 

2 

300x300 

350 

400 

AIN 

3260 

350 

400 

Mass  loading 
layer 

Density 

(kg/m3) 

Thickness 

H 

Young's 

modulus 

(GPa) 

Poisson's 

ratio 

Mechanical 
quality  factor  of 
longitudinal  mode 

Mechanical 
quality  factor  of 
Shear  mode 

Si02 

2200 

0-500 

70 

0.17 

1000 

1000 

A1 

2700 

70 

0.35 

Au 

19300 

79 

0.44 

Pt 

21450 

168 

0.38 
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TABLE  II.  Mass  sensitivity  of  FBAR  based  on  c-axis  tilted  ZnO 


c-axis 

0  (°) 

S,  (d)  (cm2/ g) 

ss(&)  (cm2/ g) 

Si02 

A1 

Au 

Pt 

Si02 

A1 

Au 

Pt 

0 

-908 

-908 

-907 

-904 

15 

-899 

-899 

-896 

-895 

-900 

-900 

-896 

-900 

30 

-883 

-883 

-881 

-881 

-925 

-925 

-922 

-923 

43 

-915 

-915 

-912 

-914 

60 

-881 

-881 

-881 

-881 

-880 

-880 

-881 

-876 

65.4 

-881 

-881 

-881 

-881 

75 

-880 

-880 

-876 

-876 

-885 

-885 

-886 

-886 

90 

-901 

-901 

-902 

-900 
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TABLE  III.  Mass  sensitivity  of  FBAR  based  on  c-axis  tilted  AIN 


c-axis 

oC) 

Sl{6)  (cm2/ g) 

Ss(0)  (cm2/ g) 

Si02 

A1 

Au 

Pt 

Si02 

A1 

Au 

Pt 

0 

-1571 

-1571 

-1565 

-1566 

15 

-1561 

-1561 

-1554 

-1557 

-1551 

-1550 

-1544 

-1543 

30 

-1542 

-1542 

-1539 

-1539 

-1573 

-1573 

-1570 

-1566 

46.1 

-1562 

-1562 

-1560 

-1557 

60 

-1537 

-1537 

-1534 

-1534 

-1548 

-1548 

-1544 

-1543 

67.1 

-1536 

-1536 

-1534 

-1534 

75 

-1537 

-1537 

-1534 

-1534 

-1536 

-1536 

-1534 

-1529 

90 

-1537 

-1537 

-1534 

-1534 
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(d) ZnO  material  properties 
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(e)  AIN  material  properties 
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FIG.  1.  The  schematic  of  FBARs  based  on  c-axis  tilted  ZnO  or  AIN  thin  film  and  coordinate 
systems:  (a)  FBAR  based  on  ZnO  or  AIN  thin  film  with  normal  c-axis  and  coordinate  system  (X1? 
X2,  X3),  (b)  FBAR  based  on  ZnO  or  AIN  thin  film  with  tilted  c-axis  and  coordinate  system  (X,  , 
X2 ,  X3 ),  (c)  the  relation  of  (Xi,  X2,  X3)  and  (Xi ,  X2 ,  X3 ),  (d)  ZnO  material  properties  in  (X3,  X2, 
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FIG.  2.  The  elastic  constants  of  ZnO  thin  film  as  a  function  of  tilt  angle:  (a)  cE  ,  (b)  cD 
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FIG.  4.  The  dielectric  constants  of  ZnO  thin  film  as  a  function  of  tilt  angle: 
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FIG.  5.  The  dielectric  constants  of  AIN  thin  film  as  a  function  of  tilt  angle:  £,s,  , 
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FIG.  10.  Resonant  frequency  shift  of  FBAR  based  on  c-axis  tilted  ZnO  film  with  different  mass 
loading. 
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FIG.  11.  Resonant  frequency  shift  of  FBAR  based  on  c-axis  tilted  AIN  film  with  different  mass 
loading. 
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Analysis  of  dual-mode  thin  film  bulk  acoustic  resonators  based  on  polar 
c-axis  tilted  wurtzite  gallium  nitride 
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We  present  the  analytical  study  of  dual-mode  thin  film  bulk  acoustic  resonators  (FBARs)  based  on 
wurtzite  gallium  nitride  (GaN)  thin  films  with  tilted  c-axis  orientation.  The  normal-plane  and 
in-plane  electric  polarization  vector  components  caused  by  the  titled  c-axis  orientation  induce  the 
coexistence  of  the  thickness  longitudinal  and  thickness  shear  modes  in  the  thin  film  resonator.  Due 
to  the  c-axis  orientation  dependency  of  material  properties  including  elastic  stiffness,  dielectric 
constants,  and  piezoelectric  coefficients,  bulk  wave  properties  of  FBARs  including  acoustic  velocity 
and  electromechanical  coupling  coefficient,  and  electrical  impedance  spectra  are  greatly  affected  by 
the  tilt  angle.  It  is  found  that  the  pure  longitudinal  mode  occurs  at  0=0°  and  66.4°,  and  the  pure 
shear  mode  occurs  at  0=42.8°  and  90°.  In  addition,  the  electromechanical  coupling  coefficient  for 
thickness  shear  mode  has  a  maximum  value  3.18%  at  0=34.2°,  which  is  larger  than  the  maximum 
value  for  pure  thickness  longitudinal  mode  2.44%  at  0=0°.  The  calculation  results  show  that  FBARs 
based  on  tilted  c-axis  wurtzite  GaN  thin  films  are  promising  for  frequency  control  device  design  and 
applications.  ©  2010  American  Institute  of  Physics.  [doi:10. 1063/1. 3427434] 


I.  INTRODUCTION 

In  the  past  years,  thin  film  bulk  acoustic  wave  resonators 
(FBARs)  have  attracted  much  attention  for  application  in  rf 
and  microwave  frequency  control  and  signal  processing 
fields  such  as  cellular  phones,  navigation,  satellite  communi¬ 
cation,  and  other  data  communication  systems  like  wireless 
local  area  networks.1  ~5  Commonly,  FBAR  consists  of  one 
thin  piezoelectric  layer  with  two  metal  electrodes  deposited 
on  both  sides.  Aluminum  nitride  (AIN)  and  zinc  oxide  (ZnO) 
are  the  mainstream  piezoelectric  materials  for  FBARs  fabri¬ 
cation,  and  their  thickness  usually  ranges  from  several  mi¬ 
crometers  down  to  tenth  of  micrometers,  which  results  in 
high  resonant  frequency  from  hundreds  of  megahertz  (MHz) 
to  over  tens  of  gigahertz  (GHz).  In  addition  to  applications  in 
electronic  circuits  as  frequency  control  devices,  FBARs  can 
also  be  used  for  sensor  application.  In  the  case  of  mass  sen¬ 
sor  application  as  the  quartz  microbalance,  it  is  well  known 
that  the  resonant  frequency  will  change  due  to  surface  mass 
change,  and  the  absolute  frequency  is  proportional  to  the 
square  of  the  operating  frequency.1'  Therefore,  using  FBAR 
structures,  with  operation  frequency  in  GHz  range,  can 
greatly  improve  the  sensitivity  of  the  device.  In  fact,  chemi¬ 
cal  sensors  based  on  FBARs  have  been  experimentally  in¬ 
vestigated  for  low  concentration  detection,  which  showed 
much  higher  sensitivity  and  resolution  comparable  to  quartz 
crystal  microbalances. 8 

FBARs  are  often  operated  in  thickness  longitudinal 
mode  since  the  ZnO  and  AIN  thin  film  materials  used  for 
device  fabrication  are  usually  with  normal-plane  c-axis  ori¬ 
entation.  For  sensor  applications  in  liquid  medium,  thickness 
shear  mode  bulk  acoustic  wave  resonators  are  needed.  This  is 
because  liquid  cannot  support  shear  wave,  the  damping 
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caused  by  the  liquid  medium  is  not  severe  for  thickness  shear 
mode  bulk  acoustic  wave  resonator;  while  for  thickness  lon¬ 
gitudinal  mode  resonator,  the  acoustic  wave  energy  is  radi¬ 
ated  into  liquid  through  compressional  motion,  leading  to 
significant  damping  in  the  liquid  medium.  Fortunately,  pre¬ 
vious  research  works  have  shown  that  thickness  shear  mode 
can  be  excited  through  the  use  of  AIN  and  ZnO  thin  films 
with  tilted  c-axis  orientation. 9-14 

Until  now,  GaN  thin  films  have  received  a  great  attention 
for  potential  high  frequency,  high  temperature,  and  high 
power  applications  due  to  the  unique  properties  of  large  for¬ 
bidden  gap,  considerable  hardness,  and  high  thermal 
conductivity.1 1-17  The  piezoelectric  properties  of  wurtzite 
GaN  also  make  it  an  exciting  material  for  manufacture  of 
FBARs.  For  example,  FBAR  with  resonant  frequency  up  to 
6.3  GHz  and  quality  factor  higher  than  1100  has  been 
achieved  using  GaN  thin  film.  The  polar  c-axis  orientation 
of  GaN  thin  film  is  usually  normal  to  substrate  surface, 
which  produces  only  thickness  longitudinal  mode  and  cannot 
be  applied  for  sensor  application  in  liquid  as  mentioned  be¬ 
fore.  Films  with  in-plane  polarization  have  also  been  fabri¬ 
cated  on  certain  substrate  crystal  such  as  tetragonal  ( y) 
LiA102  (100)  for  the  fabrication  of  structures  free  of  electro¬ 
static  fields  for  efficient  white  light-emitting  diodes.1  In  this 
paper,  we  consider  GaN  film  with  tilted  polar  c-axis  from 
normal-plane  to  in-plane  change,  and  examine  the  orienta¬ 
tion  dependence  of  elastic,  dielectric,  and  piezoelectric  prop¬ 
erty  of  wurtzite  GaN  thin  film;  further  the  characteristics  of 
dual-mode  thin  film  resonators  will  be  analyzed. 

II.  THEORETICAL  ANALYSIS 

A.  Material  properties  of  c-axis  tilted  wurtzite  GaN 

The  material  properties  of  wurtzite  GaN  thin  films  with 
tilted  c-axis  can  be  calculated  from  properties  of  wurtzite 
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FIG.  1.  The  schematic  of  FBARs  based  on  wurtzite  GaN  thin  Him  and 
coordinate  systems:  (a)  FBAR  based  on  wurtzite  GaN  thin  film  with  normal 
c-axis  and  coordinate  system  (Xj  ,X2,X3),  (b)  FBAR  based  on  wurtzite  GaN 
with  tilted  c-axis  and  coordinate  system  (XJ.XJ.XJ),  and  (c)  the  relation  of 
(Xj.X2.X3)  and  (XJ.X2.X3). 


GaN  thin  films  with  normal  c-axis  through  coordinate  system 
rotation.  Figure  1(a)  shows  FBAR  based  on  wurtzite  GaN 
thin  film  with  normal  c-axis,  and  ( X3.X3.X3 )  is  the  original 
material  coordinate;  Fig.  1(b)  shows  the  schematic  of  FBAR 
based  on  c-axis  tilted  wurtzite  GaN,  where  the  top  and  bot¬ 
tom  electrode  are  ignored.  A  rectangular  Cartesian  coordi¬ 
nate  system  (xJ.X2.X3)  is  chosen  with  top  electrode  on  x3 


=  h  and  bottom  electrode  on  x3  =  0.  The  c-axis  of  GaN  film  is 
tilted  at  an  angle  9  to  the  x3  direction.  Coordinate  system 
(xJ.X2.X3)  can  be  treated  as  the  result  of  rotation  of 
(x3.X2.X3)  about  x2  with  an  angle  9,  as  shown  in  Fig.  1(c). 
For  the  characterization  of  c-axis  tilted  wurtzite  GaN 
FBARs,  the  material  properties  in  (xj  ,x2  ,x3)  need  to  be 
known.  Assuming  the  dielectric  permittivity,  piezoelectric 
stress  coefficient,  and  elastic  stiffness  are  (e',e',c')  for  co¬ 
ordinate  system  (xJ.X2.X3),  and  (e.e.c)  for  coordinate  sys¬ 
tem  (x3.x2.X3),  the  material  properties  in  new  coordinate 
system  can  be  computed  through  properties  in  original  coor¬ 
dinate  system  with  the  aids  of  matrix  algebra 
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where  9  is  the  tilt  angle,  a'  is  the  transpose  of  a,  and  M'  is  the 
transpose  of  M. 


B.  Impedance  of  FBAR  based  on  c-axis  tilted  wurtzite 
GaN 

91 

Following  a  similar  procedure,”  the  equation  for  c-axis 
tilted  FBAR  impedance  is  derived.  According  to  the  acoustic 
field  equations  and  piezoelectric  constitution  equations,  in 
the  coordinate  system  (xj,x2,xj)  we  have'" 
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Di  =  e\qSq  +  sfkEk,  (9) 

where  Tp,  Sq,  Dh  and  Ek  are  the  components  of  stress,  strain, 

r.  / 

electric  displacement,  and  electric  field  intensity,  cpq  are  the 
elastic  stiffness  constants  under  constant  electric  field  inten- 
sity,  ekp  are  the  piezoelectric  stress  constants,  ejk  are  the 
permittivity  constants  under  constant  strain,  and  u,  is  the  dis¬ 
placement  in  the  direction  of  x\  (i,k=  1,2,3  and  p,q 
=  1,2, 3, 4, 5, 6).  The  thickness  of  GaN  is  much  smaller  than 
the  lateral  dimensions,  which  means  that  the  resonator  vibra- 
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tion  can  be  treated  as  one  dimension  problem;  hence,  we  can 
assume 


TABLE  I.  Material  properties  of  wurtzite  GaN  (Refs.  22-24). 


Property 
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For  a  sinusoidal  excitation,  the  voltage  V  and  current  7  across 
FBAR  are  expressed  by 
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From  Eqs.  (6)— ( 16),  the  impedance  of  FBAR  can  be  solved 
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Here,  /  is  frequency,  v{L)  and  u>Sl  are  the  acoustic  velocities 
for  longitudinal  mode  and  shear  mode,  and  k2L  and  k2s  are, 
respectively,  defined  as  the  electromechanical  coupling  coef¬ 
ficient  for  the  longitudinal  and  shear  mode. 
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(17g)  III.  RESULTS  AND  DISCUSSION 

A.  Dependence  of  wurtzite  GaN  properties  on  the  tilt 
angle 

(17h)  Using  Eqs.  (l)-(5)  and  material  properties  in  Table  I,  the 

material  properties  in  (x[,x2,x3)  can  be  calculated.  We  cal- 
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FIG.  2.  (Color  online)  The  elastic  constants  of  wurtzite  GaN  thin  film  as  a 
function  of  tilt  angle:  (a)  cE  and  (b)  cD  . 

culate  two  type  elastic  constants  cE  and  cD  shown  in  Fig.  2. 
cE  is  the  elastic  stiffness  constant  under  constant  electric 
field  intensity,  and  cD  is  the  elastic  stiffness  constant  under 
constant  electric  displacement.  When  0  changes  from  0°  to 
180°,  6  by  6  matrix  cE  and  cD  keeps  symmetric  with  cE4 


FIG.  3.  (Color  online)  The  dielectric  constants  of  wurtzite  GaN  thin  film  as 
a  function  of  tilt  angle:  (a) £33  and  (b)  £33. 
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L36  “c56  ‘ 


L45  ■ 


14  “c24  _c34 

0.  From  Fig.  2,  it  can  be  seen  that 


cf2  and  C22  do  not  change  with  the  tilt  angle,  which  can  be 
explained  by  the  fact  x2  and  x'2  are  the  same  axis.  However, 
the  other  nonzero  elastic  constants  are  strongly  dependant  on 


f  d'  F*  nr 

the  tilt  angle.  For  example,  cu,  cn,  c33,  and  c33  are  sym¬ 
metric  about  0=90°,  decrease  at  first,  and  increase  when  0  is 
over  45°. 

The  two  calculated  dielectric  constants  es  and  eT  ,  the 
dielectric  constant  under  constant  strain  and  under  constant 
stress,  keep  symmetric  as  3  X  3  matrixes  for  0°  <  180°, 

and  the  nonzero  items  are  shown  in  Fig.  3.  Similar  with 


dielectric  constant,  e22  and  e,2  keep  constant  because  x'2  does 
not  change  with  0.  On  the  contrary,  other  dielectric  constants 


cr 

show  great  dependency  on  tilt  angle.  The  variation  in  e , ,  or 

rjit  i  pi  rjit  pi  T>1 

en  with  6  has  an  opposite  trend  with  e33  or  e33.  or  en 
has  a  maximum  value  at  0=90°  and  a  minimum  value  at  0 
=  0°,  while  s33  or  e33  has  a  same  maximum  value  at  0=0° 
and  a  same  minimum  value  at  $=90°.  This  is  attributed  to 
the  fact  that  axis  x[  becomes  x3  and  x3  coincides  with  nega¬ 
tive  X|  when  0  is  equal  to  90°. 

The  piezoelectric  strain  constant  d'  (3X6  matrix)  and 
piezoelectric  stress  constant  e'  (3X6  matrix)  also  behavior  a 
strong  dependency  on  tilt  angle,  and  the  nonzero  constants 
are  shown  in  Figs.  4  and  5. 
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FIG.  4.  (Color  online)  The  piezoelectric  strain  constants  ( d ')  of  wurtzite 
GaN  thin  film  as  a  function  of  tilt  angle. 

B.  Bulk  acoustic  wave  properties  o  wurtzite  GaN 

Acoustic  velocity  and  electromechanical  coupling  coef¬ 
ficient  are  two  important  parameters  for  FBAR  design.  Fig¬ 
ures  6  and  7  show  the  calculated  acoustic  velocity  and  elec¬ 
tromechanical  coupling  coefficient  of  GaN  thin  films  for 
longitudinal  and  shear  mode.  The  longitudinal  velocity  (v(Lj) 
is  found  in  the  range  of  7598-8145  m/s,  which  is  larger  than 
shear  velocity  (i/5*)  range  4132-4899  m/s.  On  the  contrary, 
longitudinal  electromechanical  coupling  coefficient  ( kL 2) 
varies  from  0%  to  2.44%,  which  is  smaller  than  shear  elec¬ 
tromechanical  coupling  coefficient  ( ks 2)  range  0%  to  3.18%. 
Based  on  the  impedance  expression  of  FBAR,  the  resonator 
could  be  operated  in  a  single  mode  when  a  tilt  angle  makes 
kL=0  (pure  shear  mode)  or  ks= 0  (pure  longitudinal  mode). 
Table  II  summarizes  the  results  for  some  special  angles:  the 
pure  longitudinal  mode  occurred  at  0=0°  and  66.4°,  and  the 
pure  shear  mode  occurs  at  0=42.8°  and  90°.  The  electrome¬ 
chanical  coupling  coefficient  of  shear  mode  has  a  maximum 
value  3.18%  at  0=34.2°,  which  is  larger  than  the  maximum 
value  of  longitudinal  mode  2.44%  at  0=0°.  With  the  devel¬ 
opment  of  modern  mobile  communication  systems,  filters 
with  wider  bandwidth  are  greatly  needed,  which  requires 
higher  electromechanical  coupling  coefficient  of  piezoelec¬ 
tric  material  for  FBAR  based  filters.  Hence,  compared  with 


(b)  e  (“) 


FIG.  5.  (Color  online)  The  piezoelectric  stress  constants  ( e ')  of  wurtzite 
GaN  thin  film  as  a  function  of  tilt  angle. 

(kL2)  up  to  2.44%  (0=0°),  the  higher  value  of  ks2  up  to 
3.18%  (0=34.2°)  can  provide  a  better  choice  for  the  design 
of  higher  bandwidth  filter. 


C.  Simulation  of  electric  impedance  spectrum  of 
FBAR  based  on  wurtzite  GaN 

Based  on  the  one  dimension  impedance  equation,  the 
electric  impedance  spectrum  of  FBAR  can  be  simulated. 
Table  III  lists  the  parameters  used  in  the  simulation.  Figure  8 
shows  the  calculated  impedance  spectra  of  GaN  FBARs  with 
different  tilt  angle.  It  can  be  seen  that  there  are  two  reso¬ 
nance  peaks  for  most  angles.  The  first  peak  is  shear  reso¬ 
nance,  and  the  second  one  is  longitudinal  resonance.  The 
resonance  frequencies  of  longitudinal  and  shear  modes  are 
about  in  the  range  of  1.9  GHz  to  2.0  GHz,  and  1  GHz  to  1.2 
GHz,  respectively.  As  mentioned  before,  the  appearance  of 
strong  shear  resonance  for  some  certain  tilt  angles  makes 
GaN  FBARs  possible  for  sensor  application  in  liquid.  Con¬ 
sidering  that  low  electromechanical  coupling  coefficient  of 
shear  mode  may  lead  to  deterioration  of  resonator  operation 
in  liquid,11  pure  shear  mode  produced  at  42.8°  and  90°  with 
good  electromechanical  coupling  coefficient  is  the  best 
choice  for  FBAR  operation  in  liquid  application. 
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FIG.  6.  Acoustic  velocities  of  wurtzite  GaN  thin  film  as  function  of  tilt 
angle. 


FIG.  7.  (Color  online)  Electromechanical  coupling  coefficients  of  wurtzite 
GaN  thin  film  as  function  of  tilt  angle. 


TABLE  II.  Bulk  acoustic  wave  properties  of  wurtzite  GaN. 


Velocity  (m/s) 

Electromechanical 
coupling 
coefficient  (%) 

Tilt  angle  (°) 

Mode 

vs 

fc)2 

fe)2 

0.0 

Pure 

longitudinal 

8145 

4132 

2.44 

0 

34.2 

Quasi  shear  and 
longitudinal 

7691 

4809 

0.27 

3. 18  (max) 

42.8 

Pure  shear 

7605 

4899 

0 

2.65 

66.4 

Pure 

longitudinal 

7789 

4503 

0.41 

0 

90.0 

Pure  shear 

7963 

4161 

0 

1.40 

0.8  1.0  1.2  1.4  1.6  1.8  2.0  12.  2.4 


(a)  f  (GHz) 


(b)  f  (GHz) 


FIG.  8.  (Color  online)  Simulation  of  impedance  spectrum  of  FBAR  based 
on  c-axis  tilted  wurtzite  GaN  thin  film  (S:  shear  mode,  L:  longitudinal 
mode). 


TABLE  III.  Parameters  used  in  impedance  simulation. 


Piezoelectric  layer 

Density 

Thickness 

Electrode  area 

Mechanical  quality 
factor  of 

longitudinal  mode 

Mechanical  quality 
factor  of 

shear  mode 

Wurtzite  GaN 

6150  kg/m3 

2  /Jim 

300  X  300  /xm2 

1000 

1000 
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IV.  CONCLUSIONS 

In  summary,  FBARs  based  on  wurtzite  GaN  thin  films 
have  been  theoretically  studied.  It  is  found  that  the  acoustic 
velocity,  electromechanical  coupling  coefficient,  and  imped¬ 
ance  spectra  of  FBARs  are  greatly  dependant  on  the  c-axis 
orientation.  Strong  shear  resonant  peak  is  found  at  special  tilt 
angle.  At  0=34.2°,  the  electromechanical  coefficient  of  shear 
mode  has  a  maximum  value  3.18%;  at  0=42.8°  and  90°, 
resonator  can  be  operated  in  pure  shear  mode.  These  special 
tilt  angles  with  high  shear  electromechanical  coupling  coef¬ 
ficient  provide  more  opportunities  for  FBAR  filter  design 
and  sensor  application  in  liquid. 
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Abstract — In  this  paper,  we  present  the  analytical  study  of 
thin  film  bulk  acoustic  wave  resonators  (FBARs)  using  ZnO 
and  AIN  films  with  a  c-axis  tilt  angle  (off-normal)  from  0°  to 
180°.  The  tilted  c-axis  orientation  induces  normal  plane  and  in¬ 
plane  polarizations,  which  leads  to  the  coexistence  of  the  lon¬ 
gitudinal  mode  and  shear  mode  in  the  resonator.  The  equation 
for  predicting  electric  impedance  of  FBARs  was  derived  from 
the  basic  piezoelectric  constitutive  equations.  Material  proper¬ 
ties  including  elastic,  dielectric,  and  piezoelectric  coefficients, 
bulk  wave  properties  including  acoustic  velocity  and  electro¬ 
mechanical  coupling  coefficient,  and  impedance  of  FBARs  were 
calculated  and  showed  strong  dependence  on  the  tilt  angle. 
Interestingly,  it  was  found  that  for  ZnO  FBARs,  pure  thickness 
longitudinal  modes  occur  at  0°  and  65.4°,  and  pure  thickness 
shear  modes  occur  at  43°  and  90°.  For  AIN  FBARs,  pure  lon¬ 
gitudinal  modes  occur  at  0°  and  67.1°,  and  pure  shear  modes 
occur  at  46.1°  and  90°  for  AIN.  In  other  words,  pure  thickness 
longitudinal  and  shear  modes  exist  in  ZnO  and  AIN  FBARs 
at  specific  tilted  polarization  angles.  In  addition,  two  peaks 
of  shear  mode  electromechanical  coefficient  are  found  at  33.3° 
and  90°  for  ZnO,  and  34.5°  and  90°  for  AIN.  Therefore,  ZnO 
and  AIN  films  with  specific  tilt  angles  may  provide  options  in 
the  design  and  fabrication  of  FBARs,  considering  their  strong 
shear  resonance  with  high  electromechanical  coefficients.  The 
use  of  dual-mode  FBARs  for  mass  sensors  is  also  analyzed;  the 
calculated  large  resonant  frequency  shift  caused  by  mass  load¬ 
ing  shows  that  they  have  good  prospects  for  use  in  sensor  ap¬ 
plications  with  high  sensitivity.  The  simulation  results  agreed 
well  with  the  reported  experiment  results,  and  can  be  used  for 
design  and  application  of  FBARs. 


I.  Introduction 

RECENTLY,  in  the  RF  and  microwave  frequency  control 
and  signal  processing  fields,  thin  film  bulk  acoustic 
wave  resonators  (FBARs)  have  received  great  attention 
because  of  their  small  size,  low  insertion  loss,  and  low 
power  consumption  [1]— [5] .  For  example,  filters  based  on 
AIN  FBARs  which  can  be  operated  in  the  low  and  me¬ 
dium  gigahertz  range  have  been  fabricated  for  signal  pro¬ 
cessing  and  communication  devices  [6]— [8] .  Also,  FBARs 
show  good  prospects  for  sensor  applications.  Gabl  et  al. 
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presented  a  biosensor  system  based  on  2-GHz  longitudi¬ 
nal-mode  ZnO  FBARs  for  DNA  and  protein  molecule  de¬ 
tection;  these  sensors  showed  much  higher  sensitivity  than 
quartz  crystal  microbalances  (QCMs)  because  of  substan¬ 
tially  higher  resonant  frequencies  of  FBARs  [9] .  Zhang  et 
al.  investigated  the  effects  of  liquid  nature  and  conductiv¬ 
ity  on  ZnO  FBARs,  and  Ti02-coated  ZnO  FBARs  with 
high  mass  sensitivity  were  fabricated  for  metal  ion  K+  de¬ 
tection  in  a  liquid  environment  [10].  Benetti  et  al.  reported 
a  chemical  sensor  based  on  AIN  FBARs  with  resonance 
frequency  around  1.6  GHz,  which  was  used  to  detect  low 
concentrations  of  H2,  CO,  and  ethanol  with  a  fast  and 
repeatable  response  [11]. 

However,  for  sensor  applications  in  liquid,  the  perfor¬ 
mance  of  thickness  longitudinal-mode  FBARs  will  be  ad¬ 
versely  affected  because  the  acoustic  wave  energy  is  radi¬ 
ated  into  liquid  through  compressional  motion,  resulting 
in  excessive  damping  of  the  resonator  [12].  On  the  other 
hand,  shear  acoustic  waves  do  not  produce  compressional 
motion  and  thus  are  relatively  undamped  by  liquid.  There¬ 
fore,  shear  mode  is  a  preferred  choice  for  liquid  operation 
of  FBAR,  like  the  operation  of  a  QCM  [13].  Now  it  must 
be  determined  what  kind  of  AIN  or  ZnO  thin  films  can 
excite  the  shear  mode. 

In  recent  studies,  it  has  been  shown  that  thickness  shear 
mode  resonance  can  be  excited  in  piezoelectric  thin  films 
with  tilted  c-axis  (polarization  off-normal)  orientations. 
Foster  et  al.  theoretically  and  experimentally  analyzed  the 
excitation  of  longitudinal  and  shear  wave  in  a  ZnO  trans¬ 
ducer,  which  showed  the  shear-wave  excitation  greatly  ex¬ 
ceeded  the  longitudinal  wave  excitation  when  the  c-axis 
was  inclined  at  angles  near  40°  [14].  Wittstruck  et  al.  re¬ 
ported  the  fabrication  of  the  Mg^Zn^^O  through  metal 
organic  chemical  vapor  deposition  (MOCVD);  the  c-axis 
of  the  Mgj.Zn^.j.O  film  was  lying  in  the  plane  surface,  pro¬ 
viding  a  possibility  for  thickness  shear  wave  sensing  in  a 
liquid  environment  [15].  Yanagitani  et  al.  reported  pure 
thickness  shear-mode  FBARs  made  of  (1010)  and  (1120) 
textured  ZnO  films  [16],  [17].  Link  et  al.  fabricated  solidly 
mounted  FBARs  using  18°  c-axis  inclined  ZnO  thin  films; 
mechanical  quality  factors  of  192  were  determined  in  wa¬ 
ter,  indicating  these  devices  are  attractive  for  sensing  ap¬ 
plications  in  liquids  [18]  .Wu  et  al.  theoretically  analyzed 
bulk  acoustic  wave  properties  of  all  crystalline  plane  ori¬ 
ented  sputtered  and  epitaxial  ZnO  films  using  the  Christ- 
offel  equation;  it  showed  that  ZnO  films  with  different 
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crystalline  plane  orientation  formed  different  acoustic 
wave  modes  and  properties  [19]. 

In  the  case  of  AIN,  Bjurstrom  and  Wingqvist  et  al. 
deposited  c-axis  tilted  AIN  films  by  reactive  magnetron 
sputter  deposition  system.  The  variation  of  the  electrome¬ 
chanical  coupling  and  the  quality  factor  of  the  resonators 
with  tilt  angle  were  examined  and  a  dual-mode  FBAR 
with  a  micro-fluidic  transport  system  was  constructed 
to  investigate  the  device  performance  in  liquid  [20]- [23]. 
Martin  et  al.  presented  solidly  mounted  FBARs  based  on 
c-axis  tilted  AIN  films;  the  resonant  frequencies  of  resona¬ 
tors  could  be  high,  up  to  around  8  GHz  for  longitudinal 
mode  and  4  GHz  for  shear  mode  [24],  [25].  Chuang  et  al. 
fabricated  solidly  mounted  dual-mode  FBARs  using  0°, 
15°,  and  30°  c-axis  tilted  AIN  films,  and  experimentally 
investigated  the  effect  of  mass  loading  on  the  resonators 
[26].  Chen  et  al.  analyzed  the  input  impedance  of  a  dual¬ 
mode  AIN  FBAR  by  ignoring  the  coupling  of  longitudinal 
and  shear  mode,  which  showed  the  coexistence  of  dual 
mode  was  dependent  on  the  tilt  angle  [27]. 

In  this  paper,  we  present  the  theoretical  analysis  of 
thickness  longitudinal  and  thickness  shear  dual-mode 
FBARs  based  on  c-axis  tilted  ZnO  and  AIN  films.  Follow¬ 
ing  a  similar  procedure  [14],  the  equation  for  predicting 
electric  impedance  of  FBARS  was  derived  by  the  basic 
piezoelectric  equations.  The  material  properties  including 
elastic,  dielectric,  and  piezoelectric  constants  were  calcu¬ 
lated  with  a  tilt  angle  range  from  0°  to  180°,  which  were 
used  to  determine  the  bulk  wave  properties  of  films  in¬ 
cluding  acoustic  velocity  and  electromechanical  coupling 
coefficient.  Interestingly,  it  is  found  that  pure  thickness 
longitudinal  and  pure  thickness  shear  modes  occur  in  the 
c-axis  tilted  AIN  and  ZnO  thin  film  resonators  at  spe¬ 
cific  tilt  angles.  The  maximum  thickness  shear  mode  elec¬ 
tromechanical  coefficients  are  found  at  33.3°  and  90°  for 
ZnO,  and  34.5°  and  90°  for  AIN  FBARs.  The  response  of 
dual-mode  resonators  were  analyzed  and  compared  with 
reported  experiment  results.  In  addition,  the  sensitivity  of 
FBAR  mass  sensors  are  analyzed  by  considering  the  effect 
of  a  thin  layer  of  mass  loading  on  the  frequency  shift  of  the 
resonators.  The  analytical  and  simulation  results  indicate 
that  high  sensitivity  can  be  obtained  when  the  dual-mode 
FBARs  are  used  for  sensor  applications. 

II.  Theory 

A.  Dual-Mode  FBARs 

Fig.  1(a)  shows  the  schematic  of  tilted  c-axis  ZnO  or  AIN 
FBARs,  where  the  top  and  bottom  electrode  are  ignored.  A 
rectangular  Cartesian  coordinate  system  (aq,  x&x's)  was  cho¬ 
sen  with  top  electrode  on  13  =  h  and  bottom  electrode  on 
x.>,  =  0.  The  c-axis  of  ZnO  or  AIN  films  is  tilted  at  an  angle 
9  to  the  £3  direction.  Fig.  1(b)  shows  the  c-axis  normal  co¬ 
ordinate  (aq,^,^).  Coordinate  system  (aq,  aq,  £3)  can  be 
treated  as  the  result  of  rotation  of  (aq,oq,ag)  about  aq  with 
an  angle  9 ,  as  shown  in  Fig.  1(c). 


X, 


Xj 


Fig.  1.  The  schematic  of  FBAR  and  coordinate  systems,  (a)  Schematic 
of  FBAR  and  (:ci,  £3)  coordinate  system;  (b)  (x^  £2,23)  coordinate 

system;  (c)  the  relation  between  the  (rci,  £2,  £3)  and  (£1,22,23)  coordinate 
systems. 

B.  Material  Properties  of  ZnO  and  AIN  Films 

For  the  characterization  of  tilted  c-axis  ZnO  or  AIN 
FBARs,  the  material  properties  in  (x[,X2,Xs)  need  to  be 
known.  Assuming  the  dielectric  permittivity,  piezoelectric 
stress  coefficient,  and  elastic  stiffness  are  (s',e',c')  for  coor¬ 
dinate  system  (aq,  aq,  aq) ,  and  (e,  e,  c)  for  coordinate  system 
(aq,oq,  aq),  the  material  properties  in  new  coordinate  sys¬ 
tem  can  be  computed  through  properties  in  original  coor¬ 
dinate  system  with  the  aid  of  matrix  algebra  [28] : 

e  =  aeax  (1) 

e  =  aeM *  (2) 

d  =  McMt  (3) 


an 

Ol2 

On 

cos(0) 

0 

sin(0) 

a  = 

0-21 

022 

023 

= 

0 

1 

0 

0-31 

a  32 

a  33 

—  sin(0) 

0 

cos  (9) 

and  (5),  see  next  page,  where  9  is  the  tilt  angle,  at  is  the 
transpose  of  a,  and  M  is  the  transpose  of  M. 
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ail 

a  12 

a  13 

2ai2ai3 

2anai3 

2anai2 

all 

a  22 

a  23 

2a  22a  23 

2a2ia23 

2a2l022 

all 

a  32 

a  33 

2a 32a 33 

2a3ia33 

2a3ia32 

a2i«3i 

a22«32 

a  23a  33 

«22a33  +  a23«32 

a2i«33  +  o,23a3i 

aoia32  +  a22a3i 

ana3i 

aioa32 

a  13a  33 

ai2a33  +  ai3a32 

ana33  +  ai3a3i 

ana32  +  012031 

ana2i 

ai2tt22 

ai3a23 

ai2«23  +  ai3022 

ana23  +  ai3a2i 

ana22  +  ai2tt2i 

C.  Electric  Impedance  of  FBARs 

According  to  the  acoustic  field  equations  and  piezoelec¬ 
tric  constitution  equations  [28],  [29],  we  have: 


\Ti 

d 

d 

d 

d\ 

dx[ 

0 

0 

0 

dx3 

dx2 

t2 

dt2 

d 

d 

d 

t3 

d2uo 

0 

dx'2 

0 

dx3 

0 

dx  1 

Ta 

=  P 

dt 2 

0 

0 

d 

d 

d 

0 

t5 

d2u3 

dx3 

dx  2 

dx  1 

T6 

dt 2 

d 

dx'i 

0 

0 

0 


0 

d 

dx2 

0 

d 

dx3 


0 


d 
dx3 
d  d 
dx2  dxi 


0 

0 

d 

dx3 

d 

dx'2 

d 

dx[ 

0 


Tp  —  Cpq  Sq  ekpEk 

Pi  C-iqS  q  ~\~  £  ik  Efa 


(6) 


(7) 


(8) 

(9) 


dT\ 

dx[ 


dl\ 

dx'2 


dP3 

dx3 


=  0. 


=  0 

(10a) 

E3  0 

(10b) 

). 

(11) 

From  (6)-(10),  we  can  obtain: 

E ,  du2  E  ,du3  E ,  dUl  , 

Ta  =  c  44  ——r  +  C  34  — T  +  C  45  — - —  —  C3aE3  (12) 
OX 3  OX 3  OX 3 


8Ta  d2u2 
dx'3  ^  dt 2 


: '  du3 


dT5  d2Ui 
dx'3  P  dt 2 

:>dui 


.  i  du2 


8T3  =  d~u3 
dx3  ~  P  dt 2 

,  dui  ,  du2  ,  du3 
P3  =  e33  ——f  +  e34  ——f  +  e33  ——p 

OX  3  OX  3  OX  3 


For  any  tilt  angle  9 ,  it  is  found  that 

C34  =  C45  =  e34  =  0. 


(13) 


p  /  du\  p  t  dux  jp  r  duo  /  ,  , 

^5  —  c55  T  7  +  c35  “  7  +  c45  T  T  —  e35-^3  (14) 

<72)3  0X3  OX  3 


(15) 


^3  —  c33  ~T  7  +  c35  ~  7  +  CU  ~  T  ~  e33^3  (16) 

0X3  OX  3  <7X3 


(17) 


+  £33  E's-  (18) 


(19) 


Using  (12),  (13),  and  (19),  we  get  u2  =  0.  Substituting 
(14)  into  (15)  and  substituting  (16)  into  (17),  we  have: 


where  Tp,  Sq,  Pp  and  Ek  are  the  components  of  stress, 
strain,  electric  displacement,  and  electric  field  intensity, 
Cpq'  are  the  elastic  stiffness  constants  under  constant  elec¬ 
tric  field  intensity,  e'kp  are  the  piezoelectric  stress  constants, 
£ik  are  the  permittivity  constants  under  constant  strain, 
and  Ui  is  the  displacement  in  the  direction  of  Xi  ( i ,  k  = 
1,2,3;  p,  q  =  1,2, 3,4, 5, 6).  Because  the  thickness  of  the 
piezoelectric  layer  is  much  smaller  than  the  lateral  dimen¬ 
sions,  the  thickness  vibration  of  the  resonator  can  be 
treated  as  a  one-dimensional  problem;  hence,  we  may 
have: 


E'  d  «i 
C55  _  /  2 
dx32 


E  ’  d  u  3 

C33 


F,d2u3  d~u  1 

C35  -—n;  =  P- 


dx32 


dx32 


e<  d  U\ 

C35  5TTT  =  P 


dx3 2 


dt2 

d2u3 
dt2  ’ 


(20) 

(21) 


where 


-Jf'_cE'Ae33)2 

C33  —  C33  H - AT 
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cl' 

E'  .  e33e35 

=  C35  +  pS  ' 

£33 

(22b) 

cl' 

_  E>  ,  (e35)2 
—  C55  +  g  ,  ■ 

(22c) 

For  a  sinusoidal  excitation,  the  solution  to  (20)  and  (21) 
is  given  by 


u  l 


Ul 

sin  (a)  cos  (a) 

«L 

u3 

cos  (a)  —  sin  (a) 

us 

(23a) 


(  j4l  sin  (wo;  3 /U(L))  +  Bl  cos {uix'z/v^  )ejwt  (23b) 


es  =  e35  cos(a)  -  e33  sin  (a).  (26d) 

The  voltage  across  the  piezoelectric  layer  is  found  from 
(11),  (18),  (26c),  and  (26d)  to  be 

V  =  fhE3dx3' 

Jo 

=  -yrkW  -  «l(0)]  -  -yr[us(ft)  -  «s(0)]  + 

£33  £33  £33 

(27) 

where  A  is  the  area  of  the  electrodes.  The  current  /  can 
be  expressed  by 


u  s  =  (Hg  sin  ( 070:3/ +  Bscos(iox'3/v{S)))ejw\ 

(23c) 


d Q 
d  t 


wAD3. 


(28) 


where 


a 


w  =  27 r/ 


1 

-  arctan 
2 


2c  35 


cl' 


cl' 


(24a) 

(24b) 


,(L) 


c33 


cl' 


2  P 


v®  = 


(24c) 


1/2 


(24d) 


/  is  the  frequency,  tal)  and  u(s)  are  the  acoustic  velocities 
for  longitudinal  mode  and  shear  mode,  respectively,  and 
A^,  5l,  Hg,  Bg  are  the  unknown  amplitudes  determined  by 
the  boundary  conditions.  The  stress  can  be  solved  through 
substituting  (23)  and  (24)  into  (14)  and  (16),  given  by 


The  traction  force  at  the  boundary  ( x3  =  0 ,h)  is  zero, 
thus  we  have 


r5(0)  =  T5(h)  =  T3(  0)  =  T3(h )  =  0.  (29) 


Using  (23),  (25),  and  (29),  we  find  that 


uL(h)  -  uL( 0)  = 


us(h)  -  us{ 0)  = 


^LAstan)^) 

wZle33 

2esD3ta,n(^E) 

wZssf3' 


(30a) 

(30b) 


Hence,  substituting  (30)  into  (27)  and  using  (28),  the  im¬ 
pedance  of  FBAR  is  found  to  be 


Z  = 


V 
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JwU0 
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k? 


tan(7L/2) 

7l/2 


kl 


tan(7S/2) 

7s/2 
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where 


ZL  =  pv{L) 

Zs  =  pv^ 

eL  =  el  sin  (a)  +  e33cos(a) 


(25c) 

(26a) 

(26b) 
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where 
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7L  =  1(U 

(32a) 
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(32b) 

7S  =  ^S) 

cneo 

U) 

II 
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O 

(32c) 

(32d) 

£33  p(u(L)) 

V-  (es)2 

(32e) 

4 7«3)> 

Here,  k l  and  kg  are  respectively  defined  as  the  electrome¬ 
chanical  coupling  coefficient  for  the  longitudinal  and  shear 
modes. 
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TABLE  I.  Material  Properties  of  ZnO  and  AlN. 


Property 

Density  (kg/cm3) 
Elastic  stiffness  (Gpa) 


C 11 

C12 

C13 

0 

0 

0 

C12 

C11 

Cl3 

0 

0 

0 

C13 

Cl3 

C33 

0 

0 

0 

0 

0 

0 

C44 

0 

0 

0 

0 

0 

0 

c  44 

0 

0 

0 

0 

0 

0 

^66 

Piezoelectric  stress  constant  (C/m2) 
0  0  0  0  eis  0 

6  =  0  0  0  eis  0  0 

e3i  e3i  e33  0  0  0 

Dielectric  permittivity  (10-11  F/m) 

eft  0  0 

0  eft  0 
0  0  ef3 


AIN  ZnO 


p 

3260 

5700 

rE 

cll 

345 

210 

cE 

c12 

125 

121 

rE 

c13 

120 

105 

cE 

c33 

395 

211 

rE 

c44 

118 

43.0 

rE 

c66 

110 

44.5 

e3i 

-0.58 

-0.57 

e33 

1.55 

1.32 

eis 

-0.48 

-0.48 

£11 

8.0 

7.61 

.S' 

£33 

9.5 

8.85 

III.  Results  and  Discussion 

A.  Dependence  of  ZnO  and  AIN  Properties  on  the  Tilt 
Angle  6 

The  material  properties  of  ZnO  [30],  [31]  and  AIN  [32] 
in  the  {x\,X\,x{)  coordinate  system  are  listed  in  Table  I. 
Using  (l)-(5),  (22),  (26c),  and  (26d),  the  material  proper¬ 
ties  and  effective  parameters  in  (x\,  x%  xf)  coordinate  sys¬ 
tem  are  calculated  and  shown  in  Figs.  2  and  3,  which  can 
clearly  show  that  elastic  stiffness,  dielectric  permittivity, 
and  piezoelectric  coefficient  are  tilt-angle  dependent.  Ac¬ 
cording  to  (22),  eg’  and  eg'  can  be  regarded  as  piezoelec- 
trically  stiffened  eg  and  eg  ,  respectively;  and  have  larger 
values,  as  shown  in  Figs.  2(a)  and  3(a).  The  value  of  eg  , 
as  shown  in  Figs.  2(b)  and  3(b),  is  much  smaller  compared 

with  eg'  and  egg,  which  makes  the  calculated  angle  a 
through  (24b)  very  small,  in  a  range  of  0  to  5°  shown  in 
Figs.  2(d)  and  3(d).  a  indicates  the  strength  of  the  cou¬ 
pling  between  the  longitudinal  and  shear  modes.  The 
smaller  the  value  of  a  is,  the  weaker  is  the  coupling.  The 
small  value  of  a  makes  the  calculated  e l  and  eg  through 
(26)  very  close  to  e 33  and  e 55,  respectively,  which  can  be 
clearly  seen  from  Figs.  2(e)  and  2(f)  and  Figs.  3(e)  and 
3(f).  The  calculated  eg  ,  as  shown  in  Figs.  2(c)  and  3(c), 
varies  from  7.61  x  10-11  to  8.85  x  10-11  F/m  for  ZnO  and 
8.00  x  10-11  to  9.50  x  10“n  F/m  for  AIN. 

B.  Bulk  Acoustic  Wave  Properties  of  ZnO  and  AIN 

It  is  necessary  to  know  the  acoustic  velocity  and  elec¬ 
tromechanical  coupling  coefficient  of  ZnO  and  AIN  films 
with  different  tilt  angle,  because  they  are  two  important 
bulk  acoustic  wave  properties  for  FBAR  design.  Figs.  4(a) 
and  5(a)  show  the  acoustic  velocity  of  ZnO  and  AIN  films 
for  longitudinal  and  shear  mode  based  on  (24c)  and  (24d); 
Figs.  4(b)  and  5(b)  show  electromechanical  coupling  coef¬ 


ficient  of  ZnO  and  AIN  films  for  longitudinal  and  shear 
mode  based  on  (32d)  and  (32e).  The  longitudinal  and 
shear  velocity  of  AIN  are  found  in  the  range  of  10287  to 
11354  m/s  and  6016  to  6368  m/s,  respectively;  for  ZnO, 
the  longitudinal  and  shear  velocity  are  in  the  range  of 
5935  to  6362  m/s  and  2747  to  3233  m/s,  respectively.  For 
most  tilt  angles,  the  longitudinal  electromechanical  coef¬ 
ficient  (fcL2)of  ZnO,  varying  from  0%  to  8.53%,  is  higher 
than  that  of  AIN,  varying  from  0  to  6.0%;  and  the  shear 
electromechanical  coefficient  (fcg2)of  ZnO,  varying  from  0 
to  13.1%,  is  higher  than  that  of  AIN,  varying  from  0  to 
6.5%. 

Based  on  the  electric  impedance  expression  (31),  it  can 
be  found  that  the  conditions  for  a  single  mode  is  k^  =  0 
(for  pure  shear  mode)  or  kg  =  0  (for  pure  longitudinal 
mode),  which  corresponds  to  to  e^  =  0  or  eg  =  0  accord¬ 
ing  to  (32d)  and  (32e).  The  tilt  angle  for  pure  modes  are 
numerically  solved  and  summarized  in  Table  II  and  III, 
which  also  include  the  special  tilt  angle  where  the  maxi¬ 
mum  shear  electromechanical  coupling  coefficient  occurs. 
For  ZnO,  the  pure  longitudinal  mode  is  excited  at  9  =  0° 
with  kg2  =  8.5%  and  6  =  65.4°  with  kgf  =  1.2%;  the  pure 
shear  mode  is  excited  at  6  =  43°  with  fcg2  =  10.4%  and  9 
=  90°  with  kg2  =  6.6%;  and  the  maximum  of  electrome¬ 
chanical  coupling  coefficient  for  shear  mode  is  13.1%  at  9 
=  33.3°.  For  AIN,  the  pure  longitudinal  mode  is  excited 
at  9  =  0°  with  kg2  =  6%  and  9  =  67.1°  with  kgf  =  0.6%; 
the  pure  shear  mode  is  excited  at  9  =  46.1°  with  kg2  = 
4.6%  and  9  =  90°  with  kg2  =  2.4%;  and  the  maximum  of 
electromechanical  coupling  coefficient  for  shear  mode  is 
6.5%  at  9  =  34.5°. 

C.  Simulation  of  Electric  Impedance  Spectra  of  FBARs 
Based  on  ZnO  and  AIN  Films 

It  is  well  known  that  the  mechanical  quality  factor  (Q) 
of  thin  films  varies  with  the  thin  film  deposition  process¬ 
ing  conditions.  The  reported  Q  value  of  ZnO  and  AIN 
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Fig.  2.  ZnO  material  properties  and  effective  parameters  in  the  (x[,  x'%  X3)  coordinate  system,  (a)  C33,  C33',  c|',  and  (b)  c.33  and  033  ',  (c)  £33';  (d) 
or,  (e)  633  and  q,;  (f)  635  and  es. 


TABLE  II.  Bulk  Acoustic  Wave  Properties  of  ZnO. 


Tilt 

angle  (°) 

Mode 

Velocity  (m/s) 

Electromechanical  coupling 
coefficient  (%) 

Vl 

Vs 

( *L )2 

(h)2 

0.0 

Pure  longitudinal 

6362 

2747 

8.5 

0 

33.3 

Quasi  shear  and 

5992 

3216 

1.1 

13.1  (max) 

longitudinal 

43.0 

Pure  shear 

5936 

3213 

0 

10.4 

65.4 

Pure  longitudinal 

6028 

2916 

1.2 

0 

90.0 

Pure  shear 

6070 

2842 

0 

6.6 

TABLE  III.  Bulk  Acoustic  Wave  Properties  of  AlN. 


Tilt 

angle  (°) 

Mode 

Velocity  (m/s) 

Electromechanical  coupling 
coefficient  (%) 

vl 

vs 

{ky  (**)2 

0 

Pure  longitudinal 

11354 

6016 

6.0 

0 

34.5 

Quasi  shear  and 

10749 

6363 

0.8 

6.5  (max) 

longitudinal 

46.1 

Pure  shear 

10548 

6325 

0 

4.6 

67.1 

Pure  longitudinal 

10380 

6104 

0.6 

0 

90 

Pure  shear 

10287 

6089 

0 

2.4 
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Fig.  3.  AIN  material  properties  and  effective  parameters  in  the  (x(,  A,  x'3)  coordinate  system,  (a)  C33,  CsA  c|',  and  cff/;  (b)  c33  and  c^';  (c)  £33  ;  (d) 


a;  (e)  e'iz  and  e^;  (f)  e)5  and  es. 


is  around  several  hundred  [10],  [18],  [20],  [22].  Here,  for 
simplification,  QL(mechanical  quality  factor  of  longitudi¬ 
nal  mode)  is  assumed  to  be  350  for  ZnO  and  400  for  AIN; 
Q§  (mechanical  quality  factor  of  shear  mode)  assumed  to 
be  350  for  ZnO  and  400  for  AIN.  Thus,  the  new  acoustic 
velocity  up'  (longitudinal  mode)  and  vg  (shear  mode)  are 
adopted  for  impedance  calculation  when  Qg  and  QL  are 
taken  into  consideration  [33]: 


Vl  =  *>L  +  j 


.  VL 

J2Ql 

(33a) 

■  Vs 

J2Qs' 

(33b) 

In  addition,  the  thickness  of  ZnO  or  AIN  films  is  assumed 
to  be  2  pm  and  the  electrode  area  is  assumed  to  9  x 
10-8  m2.  Because  of  symmetry,  tilt  angle  from  0  to  90° 
is  adopted  in  the  simulation.  Figs.  6  and  7,  respectively, 
show  the  calculated  impedance  spectra  of  ZnO  and  AIN 
FBARs  with  different  tilt  angles,  where  the  resonance  fre¬ 
quencies  of  shear  and  longitudinal  mode  are  fairly  separat¬ 
ed  because  of  the  difference  of  the  acoustic  velocities.  For 
ZnO,  the  resonance  frequencies  of  thickness  longitudinal 


and  shear  modes  are  found  in  the  range  of  1.4  to  1.6  GHz 
and  0.68  to  0.81  GHz,  respectively;  the  ratio  of  longitu¬ 
dinal  and  shear  resonance  frequency  varies  from  1.84  to 
2.32.  Link  et  al.  reported  that  the  ratio  of  longitudinal  and 
shear  resonance  frequency  was  around  2.0  for  18°  c-axis 
inclined  ZnO  thin  film  [18],  which  is  consistent  with  our 
theoretically  calculated  value  2.09  ( 6  =  18°).  For  AIN,  the 
resonance  frequencies  of  thickness  longitudinal  and  shear 
modes  are  found  in  the  range  of  2.57  to  2.84  GHz  and  1.50 
to  1.60  GHz,  respectively,  and  the  ratio  of  longitudinal 
and  shear  resonance  frequency  varies  from  1.66  to  1.89. 
For  a  specific  tilt  angle  ( 6  =  28°)  of  AIN,  the  theoretical 
value  found  was  1.726,  which  agrees  well  with  1.72  calcu¬ 
lated  through  the  experiment  result  of  Wingqvist  et  al. 
[23].  For  another  tilt  angle  ( 6  =  30°)  of  AIN,  the  theoreti¬ 
cal  value  was  found  to  be  1.713,  which  was  also  close  to 
1.76  calculated  by  the  experimental  results  of  Chuang  et 
al  [26], 

It  also  can  be  seen  from  Figs.  6  and  7  that  resonance 
peaks  of  the  spectra  are  different  with  the  change  of  tilt 
angle.  Interestingly,  besides  tilt  angle  0°  (i.e. ,  normal  po¬ 
larization),  a  pure  thickness  longitudinal  mode  occurs  at 
65.4°  for  ZnO  and  67.1°  for  AIN,  with  smaller  electrome¬ 
chanical  coefficient.  In  addition,  besides  tilt  angle  90°  (i.e., 
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on 


(a) 


on 


(b) 

Fig.  4.  Bulk  acoustic  wave  properties  of  ZnO.  (a)  Acoustic  velocity  and 
(b)  electromechanical  coupling  coefficient. 

in-plane  polarization)  a  pure  thickness  shear  mode  occurs 
at  43°  for  ZnO  and  46.1°  for  AIN,  with  higher  electrome¬ 
chanical  coefficient.  The  strongest  resonance  peak  of  thick¬ 
ness  shear  mode  for  ZnO  occurs  at  33.3°,  where  the  maxi¬ 
mum  shear  mode  electromechanical  coupling  coefficient  is 
found  to  be  13.1%.  For  AIN,  the  strongest  resonance  peak 
of  thickness  shear  mode  occurs  at  34.5°,  where  the  maxi¬ 
mum  shear  mode  electromechanical  coupling  coefficient  is 
found  to  be  6.5%.  From  Figs.  4(b)  and  5(b),  it  can  also 
be  found  that  two  peaks  for  fcg2  exist  at  tilt  angle  6  = 
33.3°,  90°  for  ZnO,  and  9  =  34.5°,  90°  for  AIN.  It  is  noticed 
that  a  low  electromechanical  coupling  coefficient  of  shear 
mode  will  lead  to  deterioration  of  resonator  operation  in 
liquid  [23].  Hence,  tilt  angle  close  to  33.3°,  90°  for  ZnO 
films,  and  34.5°,  90°  for  AIN  films  are  good  options  for  sen¬ 
sor  applications  of  FBARs  in  liquid  environment.  In  fact, 


eC) 


(a) 


0  (°) 

(b) 

Fig.  5.  Bulk  acoustic  wave  properties  of  AIN.  (a)  Acoustic  velocity  and 
(b)  electromechanical  coupling  coefficient. 


Wingqvist  et  al.  [23]  presented  AIN  FBARs  with  a  tilt 
angle  around  30°;  they  showed  a  strong  shear  resonance 
around  1.2  GHz  with  high  Qg  of  around  150  in  the  water, 
and  the  viscous  load  experiments  indicated  the  potential 
of  FBARs  for  biosensor  bioanalytical  tools  as  well  as  for 
liquid  sensing  in  general. 

D.  Sensitivity  Analysis  of  ZnO  and  AIN  FBAR  Mass 
Sensors 

When  bulk  acoustic  wave  resonators  are  used  for  mass 
sensors  such  as  in  the  case  of  the  commonly  used  QCM, 
it  is  well  known  from  the  Sauerbrey  equation  that  the 
resonant  frequency  will  change  due  to  surface  mass  load¬ 
ing  [34],  [35],  A /=  -(2/o2Aro)/(4(/iqpq)1/2),  where  /0  is 
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PL  =  1 


eL  sin2 1 

(U 

(Sfc°s(7S)  +  jfSin(7S)) 

1  +  jegsin2! 

(U 

lfsin(7L)] 

eL  x  1 

( cos(7L)  c°s  (7  s)  -  sin  (7  L)  sin  (7  s)  +  j  1 

f  zkL 

\  ZL 

cos  (7  L)  sin  (7  s)  +  ^7  sin  (7  L)  cos  (7  s))  j 

(35g) 


pb  =  1  + 


essin2(^)(^gcos(7L)  +  +  jeLsin2(^)§j-sin(7S)] 

es  X  (§^cos(7L)cos(7S)  -  sin (7  L)  sin (7  s)  +  j(^-cos(7L)sin(7S)  +  sin  (7  L)  cos  (7  s) 


(35h) 


the  resonant  frequency  (in  hertz),  A/ is  frequency  change 
(in  hertz),  Am  is  the  mass  change  (in  grams),  A  is  the 
piezoelectrically  active  crystal  area  (the  area  between  elec¬ 
trodes,  in  square  centimeters),  pq  is  the  density  of  quartz 
(pq  =  2.648  g/cm3),  and  pq  is  the  shear  modulus  of  quartz 
for  AT-cut  crystal  (pq  =  2.947  x  1011  g/cm-s2).  The  abso¬ 
lute  frequency  change  is  proportional  to  the  square  of  the 
operating  frequency.  Because  the  operational  frequency  of 
FBARs  can  be  up  to  tens  of  gigahertz,  it  is  believed  that 
sensor  sensitivity  will  be  greatly  improved.  For  analyzing 
sensor  performance  of  the  FBARs,  the  electric  impedance 
expression  of  an  FBAR  with  a  mass  loading  layer  needs  to 
be  derived.  Acoustic  wave  resonator  sensors  usually  con¬ 
sist  of  four  layers:  one  piezoelectric  layer,  two  electrode 
layers,  and  a  surface  loading  (mass)  layer,  as  shown  in  Fig. 
8.  Usually  the  electrode  layer  is  quite  thin  and  its  acoustic 
impedance  is  very  small  compared  with  the  piezoelectric 
and  mass-loading  layers.  Thus  for  simplicity,  the  effect 
of  electrodes  is  neglected  here.  In  other  words,  we  only 
consider  the  effect  of  the  mass  layer  on  the  resonance  fre¬ 
quency  shift  of  the  dual-mode  FBARs.  For  derivation  of 
the  impedance  expression,  (l)-(28)  are  adopted  here,  but 
the  traction  forces  at  the  boundary  ( x 3'  =  0  ,h)  change  to 

U  0)  =  r3(0)  =  0  (34a) 


T5(h) 

-juZl i 

0 

ui(h) 

Hh) 

0 

-juzh 

u3(h) 

Z\  and  Zg j  are  the  longitudinal  and  shear  acoustic  imped¬ 
ance  of  the  surface  mass-loading  layer,  which  can  be  ex¬ 
pressed  by 

Zsi  =  jPsivli  tan  ( 7 si )  (34c) 

ZsSi  =  JpsiWsSitan(7ss1)  (34d) 

^sSi 
/si  ’ 

(34e) 

where  v^\  and  uSi  are  longitudinal  and  shear  acoustic  veloc¬ 
ity  of  mass-loading  layer;  7^  and  7s  are  longitudinal  and 
shear  phase  delay  in  the  layer;  psl  and  4i  are  the  density 
and  thickness  of  the  mass  layer;  c33  and  c|3  are  the  elastic 
constants  of  the  mass  layer.  According  to  the  boundary 


condition  (34)  and  (1)— (28) ,  the  electric  impedance  of  the 
mass-loaded  FBAR  can  be  solved: 


Z  = 


1 


juC0 


l-(kLf  pL  W7l/2)  _(feS)2  g  tan(7s/2) 


where 


7l/2 


C  n  = 


U-L)  = 


£33^4 

h 

(eL)2 


£33  'p(u(L0 


(fcs)2=  (M 


£3S3>(S))2 


7l  = 


7s  = 


uih 

uih 

IS)’ 


7s/2 


(35a) 

(35b) 

(35c) 

(35d) 

(35e) 

(35f) 


with  pL  and  ps  as  defined  in  (35g)  and  (35h),  see  above, 
and 

ZgiL  =  Z] 1  cos2(a)  +  Zss  sin2(a)  (35i) 

^SiS  =  Zscos2(a)  +  ZsIj,sin2(a)  (35j) 

Zg] iL  =  (Zss  —  ^d)sin(a)cos(a).  (35k) 


Through  analyzing  (35),  the  condition  for  pure  longitudi¬ 
nal  mode  is  a  =  0,  es  =  0,  which  is  equivalent  to  =  0, 
e35  =  0;  and  the  condition  for  pure  shear  mode  is  a  =  0, 
eL  =  0,  which  is  equivalent  to  C35  =  0,  e'33  =  0.  Hence, 
from  Figs.  2  and  3,  we  can  conclude  that  the  pure  longitu¬ 
dinal  mode  for  ZnO  and  AIN  is  excited  at  9  =  0°,  and  for 
pure  shear  mode  it  is  6  =  90°.  In  addition,  a  is  small  be¬ 
cause  of  the  small  value  of  C35',  then  longitudinal  mode  is 
predominant  when  es  =  0  and  shear  mode  is  predominant 
when  eL  =  0;  therefore,  at  6  =  65.4°  for  ZnO  and  9  —  67.1° 
for  AIN,  it  can  be  approximated  by  pure  longitudinal 
mode;  at  9  =  43°  for  ZnO  and  9  =  46.1°  for  AIN,  it  can  be 
approximated  by  pure  shear  mode. 

For  a  better  understanding  of  the  dual-mode  FBAR  for 
mass  sensors,  we  gradually  increase  the  mass  loading  and 


abs(Z)  (£2)  abs(Z)  (£2) 
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f  (GHz) 


f  (GHz) 

(b) 

Fig.  6.  Impedance  spectrum  simulation  of  FBAR  based  on  c-axis  tilted 
ZnO  film,  (a)  0  =  0,  15,  33.3,  and  43°;  (b)  0  =  60,  65.4,  75,  and  90°. 


f  (GHz) 


1.0  1.5  2.0  2.5  3.0  3.5  4.0 

f  (GHz) 

(b) 

Fig.  7.  Impedance  spectrum  simulation  of  FBAR.  based  on  c-axis  tilted 
AIN  film,  (a)  0  =  0,  15,  34.5,  and  46.1°;  (b)  0  =  60,  67.1,  75,  and  90°. 
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Fig.  9.  Resonant  frequency  shift  of  FBAR  based  on  c-axis  tilted  ZnO 
film  vs.  thickness  of  mass  loading,  (a)  0  =  0,  43,  65.4,  and  90°;  (b)  8  = 
15  and  30°. 


investigate  the  change  of  the  resonant  frequency  /s,  which 
corresponds  to  maximum  conductance  and  is  usually  ad¬ 
opted  for  sensor  output.  Here,  /s l(9)  and  /sg(0)  are  used 
to  express  the  longitudinal  and  shear  resonant  frequency 
of  FBAR  based  on  AIN  or  ZnO  with  c-axis  tilt  angle  9. 
Figs.  9  and  10  show  the  change  of  fsh(9)  and  fss{9)  with 
the  thickness  of  the  mass-loading  layer  for  some  specific 


Fig.  10.  Resonant  frequency  shift  of  FBAR  based  on  c-axis  tilted  AIN 
film  vs.  thickness  of  mass  loading,  (a)  0  =  0,  46.1,  67.1,  and  90°;  (b)  6 
=  15  and  30°. 


tilt  angles.  In  the  simulation,  A1  and  Ag  are  the  materi¬ 
als  for  the  mass-loading  layer,  and  Table  IV  lists  their 
parameters  for  calculation  of  /sl(0)  and  /sS(0)  using  (35). 
As  shown  in  Figs.  9  and  10,  for  mass  loading  A1  and  Ag, 
both  /sl(0)  and  /ss(0)  linearly  decrease  when  the  thickness 
of  mass-loading  layer  is  small.  For  the  same  thickness  of 
mass  loading,  the  resonant  frequency  shift  for  Ag  is  higher 
than  Al,  because  Ag  has  higher  density  than  Al,  inducing 
higher  acoustic  loading.  The  mass  sensitivity  of  dual-mode 
FBARs  is  defined  as 


Sl(0) 


A  U9)/U9) 
Am 


(36a) 


sm  = 


A  Uff)/UQ) 
Am 


(36b) 
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TABLE  IV.  Parameters  for  Calculation  of  fsL(6)  and  fsS(0). 


Mass-loading 

layer 

Density 

(kg/m3) 

Thickness 

(nm) 

Young’s 

modulus 

(GPa) 

Poisson’s 

ratio 

Ql 

Qs 

A1 

2700 

0-500 

70 

0.35 

1000 

1000 

Ag 

10490 

0-500 

83 

0.37 

1000 

1000 

TABLE  V.  Mass  Sensitivity  of  FBAR  Based  on  c-Axis  Tilted  ZnO. 

Tilt 

MB) 

MB) 

MB)  (cm 

2/g) 

MB)  (cm2/ g) 

angle  (°) 

(GHz) 

(GHz) 

A1 

Ag 

A1 

Ag 

0 

1.533439408 

— 

-904 

-904 

— 

— 

15 

1.520757366 

0.709173538 

-896 

-894 

-896 

-896 

30 

1.496097496 

0.751143640 

-882 

-879 

-926 

-919 

43 

— 

0.767714844 

— 

— 

-911 

-909 

65.4 

1.499892844 

— 

-881 

-877 

— 

— 

90 

— 

0.690955092 

— 

— 

-896 

-896 

TABLE  VI.  Mass  Sensitivity  of  FBAR  Based  on  c-Axis  Tilted  AlN. 


Tilt 

angle  (D) 

MB) 

(GHz) 

MB) 

MB)  (cm2/ g) 

MB)  (cm2/ g) 

(GHz) 

A1 

Ag 

A1 

Ag 

0 

2.767650148 

— 

-1570 

-1558 

— 

15 

2.746620142 

1.517122380 

-1556 

-1548 

-1548 

-1539 

30 

2.695601826 

1.541801144 

-1541 

-1531 

-1570 

-1560 

46.1 

— 

1.550894306 

-1563 

-1548 

67.1 

2.588555566 

— 

-1548 

-1525 

— 

— 

90 

— 

1.507487836 

-1548 

-1535 

Am  —  (36c) 

where  Sl(0)  the  mass  sensitivity  of  longitudinal  mode 
for  FBAR  based  on  AIN  or  ZnO  with  c-axis  tilt  angle  9 
and  5g(0)  is  the  mass  sensitivity  of  shear  mode;  A/sL(0) 
and  AfsS(9)  are  the  resonant  frequency  shifts  caused  by 
mass  loading  Am;  pm,  and  fm  are  the  density  and  thick¬ 
ness  of  mass- loading  layer.  The  mass  sensitivity  of  FBARs 
in  a  linear  range  (A fs(9)/fs(9)  <  5%)  is  evaluated  through 
examining  resonant  frequency  shift  for  5  nm  A1  or  Ag 
mass  loading,  and  the  calculated  mass  sensitivity  for  some 
special  angles  according  to  (36)  is  listed  in  Tables  V  and 
VI.  It  can  be  seen  from  Tables  V  and  VI  that  Sl(0)  and 
S$(9)  of  ZnO  are  around  —900  cm2/g,  and  8^(9)  and  S$(9) 
of  AIN  are  around  —1550  cm2/g,  which  is  much  higher 
than  —14  cm2/g,  the  mass  sensitivity  of  a  6-MHz  QCM 
[36].  Because  of  the  small  value  of  a  or  C35  for  all  9 ,  the 
coupling  between  longitudinal  and  shear  mode  is  weak,  so 
the  resonator  can  be  approximated  by  a  simple  combina¬ 
tion  of  two  single  modes.  On  the  other  hand,  it  is  known 
that  the  mass  sensitivity  (Sm)  of  single  mode  resonator  is 
dependent  on  the  density  (pp)  and  thickness  (dp)  of  the 
piezo  layer  [36]. 

Sm  =  (37) 

Ppdp 

Hence,  for  ZnO  or  AIN  FBAR,  Sl(0)  and  8$(9)  have  simi¬ 
lar  values,  and  do  not  change  much  with  tilt  angle,  which 
also  can  be  seen  in  Table  V  and  VI. 


It  should  be  pointed  out  that  in  (35)  the  mass  loading 
layer  is  assumed  to  be  isotropic.  When  this  layer  is  the 
material  used  for  electrodes,  (35)  can  be  used  to  analyze 
the  effect  of  electrode  thickness,  which  actually  has  been 
shown  in  mass  sensitivity  calculations.  Moreover,  through 
replacing  the  acoustic  impedance  of  the  single  mass-load¬ 
ing  layer  with  the  acoustic  impedance  of  a  multilayer  load¬ 
ing,  (35)  can  be  used  for  analysis  of  multilayer  loading 
on  the  frequency  change  of  the  resonator,  as  long  as  the 
multilayer  is  isotropic  and  non-piezoelectrically  active. 


IV.  Conclusion 

FBARs  based  on  ZnO  and  AIN  films  with  tilted  c-axis 
orientation  have  been  theoretically  studied.  Because  of  the 
crystalline  orientation  dependence  of  the  material  proper¬ 
ties,  including  elastic  stiffness,  piezoelectric  coefficients, 
and  dielectric  constants,  the  acoustic  velocity,  electro¬ 
mechanical  coefficient,  and  electric  impedance  of  FBARs 
all  depend  on  the  crystalline  orientation.  The  equation 
for  predicting  the  electric  impedance  of  FBARs  has  been 
derived  through  the  basic  piezoelectric  equations,  which 
shows  the  coexistence  of  longitudinal  and  shear  modes 
when  the  c-axis  tilt  angle  changes.  The  simulation  results 
show  that  a  pure  longitudinal  mode  occurs  at  0°  and  65.4°, 
and  a  pure  shear  mode  occurs  at  43°  and  90°  for  ZnO;  a 
pure  longitudinal  mode  occurs  at  0°  and  67.1°,  and  a  pure 
shear  mode  occurs  at  46.1°  and  90°  for  AIN.  Two  peaks 
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of  shear  electromechanical  coefficient  are  found  at  6  = 
33.3°,  90°  for  ZnO  and  0  =  34.5°,  90°  for  AIN.  ZnO  and  AIN 
films  with  a  tilt  angle  around  these  two  peaks  are  good 
options  for  FBAR  applications  in  liquid,  considering  their 
strong  shear  resonance  and  high  electromechanical  coeffi¬ 
cients.  The  mass  sensitivity  of  FBAR-based  ZnO  and  AIN 
are  much  higher  than  that  of  QCMs,  which  shows  good 
promise  for  use  as  mass  sensors. 
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Abstract — A  quartz  thickness  shear  mode  (TSM)  bulk 
acoustic  wave  resonator  was  used  for  in  situ  and  real-time  de¬ 
tection  of  liquid  flow  rate  in  this  study.  A  special  flow  chamber 
made  of  2  parallel  acrylic  plates  was  designed  for  flow  measure¬ 
ment.  The  flow  chamber  has  a  rectangular  flow  channel,  2  flow 
reservoirs  for  stabilizing  the  fluid  flow,  a  sensor  mounting  port 
for  resonator  holding,  one  inlet  port,  and  one  outlet  port  for 
pipe  connection.  A  5-MHz  TSM  quartz  resonator  was  edge- 
bonded  to  the  sensor  mounting  port  with  one  side  exposed  to 
the  flowing  liquid  and  other  side  exposed  to  air.  The  electrical 
impedance  spectra  of  the  quartz  resonator  at  different  volu¬ 
metric  flow  rate  conditions  were  measured  by  an  impedance 
analyzer  for  the  extraction  of  the  resonant  frequency  through 
a  data-fitting  method.  The  fundamental,  3rd,  5th,  7th,  and  9th 
resonant  frequency  shifts  were  found  to  be  around  920,  3572, 
5947,  8228,  and  10300  Hz  for  flow  rate  variation  from  0  to 
3000  mL/min,  which  had  a  corresponding  Reynolds  number 
change  from  0  to  822.  The  resonant  frequency  shifts  of  differ¬ 
ent  modes  are  found  to  be  quadratic  with  flow  rate,  which  is 
attributed  to  the  nonlinear  effect  of  quartz  resonator  due  to 
the  effective  normal  pressure  imposing  on  the  resonator  sen¬ 
sor  by  the  flowing  fluid.  The  results  indicate  that  quartz  TSM 
resonators  can  be  used  for  flow  sensors  with  characteristics  of 
simplicity,  fast  response,  and  good  repeatability. 


I.  Introduction 

A  flow  sensor  is  a  device  for  sensing  the  rate  of  fluid 
flow.  It  is  a  critical  component  for  system  operation 
monitoring  and  control  in  many  industrial  and  laboratory 
applications.  Typically,  a  flow  sensor  is  the  sensing  ele¬ 
ment  used  in  a  flow  meter  or  flow  logger  to  record  the  flow 
of  fluids.  There  are  various  types  of  flow  sensors  and  flow 
meters  in  practical  applications.  Traditional  mechanical 
devices,  such  as  variable  area  flow  meters,  pitot  tubes, 
and  turbine  meters,  are  simple  and  do  not  require  power 
supply  [1];  however,  they  are  limited  by  their  physical  size 
and  mechanical  characteristics,  especially  for  applications 
that  require  integration  of  electronic  control.  Some  inva¬ 
sive  methods  such  as  ultrasonic  flow  meters  [2]  and  la¬ 
ser  doppler  anemometers  [3]  are  nonintrusive  and  highly 
sensitive  to  flow  rate,  however,  they  are  expensive  and 
require  complex  auxiliary  electronic  parts.  Additionally, 
the  widely  used  thermal  anemometers  such  as  hot-wire 
anemometer,  can  measure  the  flow  in  a  precise  location 
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and  has  extremely  high  frequency  response.  However,  it 
also  has  many  problems  in  device  design  and  applications, 
including  breakage  of  filament  stability,  calibration,  and 
auxiliary  electronics  [4]. 

Acoustic  wave  resonators  including  BAW  and  SAW 
devices  have  been  regarded  as  the  most  sensitive  sensor 
platforms  that  take  the  advantage  of  precise  frequency 
counting  in  electronic  measurement,  and  they  have  been 
extensively  studied  for  many  physical-  and  chemical¬ 
sensing  applications  [5].  Acoustic  wave  resonators  have 
been  investigated  for  flow  rate  measurement,  and  flow 
sensors  based  on  SAW  resonators  have  been  demon¬ 
strated  with  attractive  features  of  high  sensitivity,  wide 
dynamic  range,  and  direct  digital  output  [6],  [7].  This 
kind  of  SAW  flow  sensor  is  based  on  thermal  convec¬ 
tion:  a  SAW  device  with  a  heater  is  placed  in  the  flow 
path.  When  the  fluid  flow  passes  over  the  sensor  surface, 
it  carries  heat  and  changes  the  temperature  of  the  sen¬ 
sor  substrate.  Thus,  the  frequency  of  the  SAW  device 
changes  due  to  its  frequency-temperature  dependence. 
Although  the  theoretical  aspect  of  this  kind  of  sensor 
is  simple  and  well  established,  the  necessity  of  a  heater 
for  sensor  operation  may  impose  practical  difficulties  in 
which  large  power  consumption  is  usually  not  desirable 
for  compact  sensor  operation.  In  addition,  the  variation 
of  ambient  temperature  and  heater  power  supply  may 
cause  instability  of  frequency  output;  to  minimize  these 
effects,  an  additional  SAW  device  that  is  not  exposed  to 
the  fluid  flow  is  needed  as  a  reference,  which  increases 
system  complexity. 

On  the  other  hand,  to  the  best  of  our  knowledge, 
BAW  resonators  have  not  been  investigated  for  flow  sen¬ 
sors.  Among  all  the  acoustic  wave  sensor  configurations, 
BAW  resonators,  particularly  quartz  thickness  shear  mode 
(TSM)  resonators,  are  often  preferable  in  many  sensor  ap¬ 
plications  due  to  its  high  sensitivity,  simple  structure,  and 
easy  interconnection  with  electronic  measurement  systems 
[8]— [12] .  Commonly,  a  resonator  consists  of  one  thin  layer 
of  piezoelectric  material  with  2  electrodes  on  both  sides. 
For  sensor  applications,  an  external  chemical  or  physi¬ 
cal  sensitive  thin  film  is  usually  coated  on  the  top  of  the 
resonator.  When  an  ac  electric  field  is  applied  across  the 
resonator  electrodes,  acoustic  waves  are  generated  and 
propagate  in  the  device  due  to  the  piezoelectric  property 
of  quartz  crystal.  Variations  in  coating  thin-film  param¬ 
eters,  including  the  thickness,  density,  and  complex  shear 
modulus  and  interfacial  interactions,  affect  the  amplitude 
and  phase  of  the  acoustic  wave  propagation  in  the  device, 
which  changes  the  electrical  properties  of  the  resonator. 
Therefore,  the  measured  quantity  can  be  precisely  ex- 
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tracted  by  the  resonance  characteristics  and  the  electrical 
spectrum  of  the  resonator. 

The  electrical  properties  of  the  resonator,  such  as  res¬ 
onant  frequency  shift  and  mechanical  damping,  depend 
not  only  on  the  mechanical  properties  of  the  coating  layer 
and  the  surface  interactions,  but  also  on  mechanical  or 
electrical  loading  conditions  because  the  material  proper¬ 
ties  of  quartz  resonators  are  subjected  to  changes  under 
electrical  or  mechanical  fields.  For  example,  it  is  found 
that  the  resonant  frequency  of  quartz  resonators  will  shift 
due  to  nonlinear  elastic  responses  when  it  is  subjected 
to  bias  fields  [13]-  [19] .  Among  these  studies,  Ballato  et 
al.  revealed  that  the  resonant  frequency  shift  of  a  single- 
rotated  Y-cut  quartz  resonator  due  to  the  applied  force 
on  the  rim  is  dependent  on  the  force  magnitude,  crystal 
cut  angle,  azimuth  angle,  and  the  vibration  mode  [13]. 
The  force-frequency  effect  of  double-rotated  quartz  reso¬ 
nators  was  also  investigated  [14].  Ratajski  introduced  the 
force-frequency  coefficient  Kt  for  quantifying  the  force- 
frequency  effect,  and  an  empirical  chart  based  on  the 
normalized  experimental  results  was  derived,  which  can 
be  used  to  determine  the  force-frequency  behavior  of  any 
crystal  cut  in  the  singly  rotated  Y-cut  group  [15].  Lee  et 
al.  derived  2-D  governing  equations  for  theoretical  analy¬ 
sis  of  force-frequency  effect;  and  the  formula  for  predicting 
the  resonant  frequency  shift  was  developed  and  compared 
with  experimental  results  [16],  [17].  Numerical  methods 
including  finite  element  method  (FEM)  were  also  used  to 
investigate  the  force-frequency  effect  of  quartz  resonators 
[18],  [19].  When  the  acoustic  wave  resonators  are  used 
for  frequency  control  and  signal  processing  applications, 
it  is  often  desirable  to  minimize  force- frequency  effect.  On 
the  other  hand,  through  the  force-frequency  effect,  quartz 
acoustic  wave  resonators  have  also  found  applications  as 
force  sensors,  pressure  sensors,  or  accelerometers  [20]— [22] . 
Therefore,  it  is  possible  to  apply  the  resonator  for  flow 
measurement  if  the  mechanical  load  due  to  the  flow  can 
significantly  change  the  frequency  response  of  acoustic 
wave  resonators  due  to  linear  or  nonlinear  response  of  the 
device  by  choosing  an  appropriate  flow  path  and  mechani¬ 
cal  boundary  conditions. 

In  this  paper,  a  quartz  acoustic  wave  resonator  is  stud¬ 
ied  as  flow  sensor  for  in  situ  and  real-time  detection  of 
fluid  (water)  flow  rate  at  room  temperature.  Considering 
the  structural  simplicity  and  high  mechanical  quality  fac¬ 
tor,  a  quartz  TSM  resonator  is  particularly  chosen  as  the 
sensing  element.  For  flow  rate  measurement,  the  quartz 
resonator  is  edge-bonded  at  the  sensor  mounting  port  in 
a  specially  designed  flow  chamber.  First,  the  relationship 
of  the  resonant  frequency  shift  to  the  flow  rate  is  theo¬ 
retically  discussed  based  on  the  nonlinear  response  of  the 
resonator.  The  resonant  frequency  shifts  of  fundamental 
thickness  shear  mode  and  overtones  are  then  monitored 
as  the  sensor  responses  to  different  flow  rates.  The  sensi¬ 
tivity  and  repeatability  of  the  acoustic  wave  flow  sensor 
are  experimentally  studied.  Finally  a  brief  discussion  is 
given  on  the  flow  sensor  optimization  and  potential  ap¬ 
plications. 


Resonator  Mounting  Port 


Fig.  1.  Schematic  of  flow  chamber:  (a)  top  and  bottom  plate  and  (b) 
cross-sectional  view. 


II.  Material  and  Methods 
A.  Quartz  TSM  Resonator  and  Flow  Chamber 

A  5-MHz  quartz  At-cut  TSM  resonator  (Stanford  Re¬ 
search  Systems,  Inc.,  Sunnyvale,  CA)  with  25.4-mm  di¬ 
ameter  was  used  in  this  study  as  a  flow  sensor.  A  flow 
chamber  was  designed  according  to  the  method  proposed 
by  Chuang  et  al.  [23].  Fig.  1  shows  the  schematic  of  the 
chamber,  which  is  composed  of  2  parallel  acrylic  plates, 
in  which  a  rectangular  flow  channel,  inlet  and  outlet  res¬ 
ervoir,  1  sensor  mounting  port,  and  1  inlet  port  and  out¬ 
let  port  are  formed.  These  2  plates  were  held  together 
through  bolts  and  an  O-ring  between  them  to  prevent 
the  liquid  leakage.  The  top  plate  was  machined  to  form  a 
sensor  mounting  port,  and  the  quartz  resonator  was  edge- 
bonded  at  the  port  using  epoxy  of  Eccobond  45  Clear 
(Ellsworth  Adhesives,  Germantown,  WI)  such  that  one 
surface  of  the  resonator  is  exposed  to  the  flowing  liquid 
and  the  other  surface  is  exposed  to  the  air.  Two  fine  Cu 
wires  were  bonded  to  the  resonator  electrodes  through  sil¬ 
ver  paste,  which  provides  electrical  connection  for  reso¬ 
nator  impedance  measurement.  The  rectangular  flow 
channel,  inlet  reservoir,  outlet  reservoir,  inlet  port,  and 
outlet  port  were  formed  in  the  bottom  of  the  plate.  Dur¬ 
ing  flow  measurement,  fluid  entered  the  inlet  reservoir 
through  the  inlet  port,  passed  through  the  rectangular 
flow  channel,  flowed  into  the  outlet  reservoir,  and  exited 
through  the  outlet  port.  Plastic  pipes  were  connected  to 
inlet  and  outlet  ports  for  flow  delivery.  The  dimension 
of  the  flow  channel  is  165.10  mm  x  69.85  mm  x  0.51  mm 
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Fig.  2  shows  the  measurement  setup,  which  is  made  of 
a  water  pool,  Manostat  Preston  pump  (Thermo  Scientific, 
Barrington,  IL),  pulse  dampener  (Cole-Parmer  Instrument 
Company,  Vernon  Hills,  IL),  flow  chamber,  Agilent  4294A 
precision  impedance  analyzer  (Agilent  Technologies,  Palo 
Alto,  CA),  Hi-Speed  GPIB  Controller  (National  Instru¬ 
ments,  Austin,  TX),  and  computer.  The  pool  was  used 
to  store  the  fluid:  deionized  (DI)  water,  which  was  driven 
by  the  pump.  The  water  level  in  the  pool  is  a  little  higher 
than  the  flow  chamber,  which  produces  around  30  cm  of 
water  pressure  on  the  resonator.  The  pulse  dampener  was 
adopted  to  smooth  the  flow  considering  the  unstable  flow 
due  to  the  pulse  pump.  The  resonator  admittance  imped¬ 
ance  spectra  near  the  resonant  frequencies  were  recorded 
by  the  impedance  analyzer.  By  data  fitting  the  admit¬ 
tance  impedance  spectra,  the  resonant  frequencies  were 
precisely  extracted;  1st  (fundamental),  3rd,  5th,  7th,  and 
9th  series  resonant  frequencies  as  the  sensor  response  were 
monitored  and  recorded  during  flow  measurement.  Water 
flow  rates  up  to  3000  mL/min  were  used  in  this  study,  all 
of  which  are  limited  to  laminar  flow;  corresponding  Reyn¬ 
olds  numbers  were  0  to  822. 


(b) 

Fig.  3.  Schematic  cross-sections  of  the  flow  channel:  (a)  yz  plane  and 
(b)  xz  plane. 
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The  resonator  is  located  at  (0,  ft/ 2,  0).  Thus,  when  Q  is 
constant,  the  normal  pressure  difference  A P  on  the  2  sides 
of  the  resonator  can  be  obtained  as: 


C.  Method  of  Flow  Measurement 

Fig.  3  shows  the  schematic  cross  sections  of  the  flow 
channel,  and  the  Cartesian  coordinate  system  (x,  y ,  z) 
for  fluid  flow  discussion.  The  top  and  bottom  plate  is  at 
y  =  h/2  and  y  =  — ft/2.  The  lateral  walls  are  at  z  =  ±w/2. 
In  the  rectangular  flow  channel,  the  relation  of  pressure  p 
and  volumetric  flow  rate  Q  is  given  by  [23] ,  [24] : 
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where  p,  L1  w,  and  ft  are  the  viscosity  of  fluid,  the  length, 
width,  and  height  of  the  flow  channel.  Shear  stress  r  at  the 
wall  (y  =  ±ft/2)  is  given  by  [23]: 


where  Lp  and  Rp  are  the  length  and  radius  of  the  pipe 
connected  to  the  outlet  port.  The  first  term  in  (4)  repre¬ 
sents  the  pressure  drop  from  the  resonator  position  to  the 
outlet  port  in  the  flow  channel,  and  the  second  term  is  the 
pipe  pressure  drop  from  the  outlet  port  to  the  water  pool. 
Also,  The  shear  stress  ry  on  the  resonator  surface  can  be 
obtained  by  substituting  the  resonator  position  into  (2): 
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Volumetric  Flow  Rate  (mL/min) 

Fig.  4.  The  normal  pressure  difference  and  wall  shear  stress  versus  volu¬ 
metric  flow  rate. 
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Eqs.  (4)  and  (5)  indicate  that  both  A P  and  rr  are  propor¬ 
tional  to  the  flow  rate.  Fig.  4  plots  the  calculation  results 
of  AP  and  rr  verse  flow  rate. 

When  the  device  is  subjected  to  fluid  flow,  the  resona¬ 
tor  will  be  deformed  due  to  AP  and  rr.  In  Fig.  5(a)  the 
dotted  line  indicates  the  middle  plane,  which  has  a  flexure 
deflection  w.  During  the  deformation,  the  resonator  disk  is 
subjected  to  either  tension  or  compression  on  the  2  sides 
of  the  middle  plane.  According  to  the  first-order  perturba¬ 
tion  integral  [25],  this  antisymmetrical  deformation  has 
no  effect  on  the  resonant  frequencies.  However,  the  stretch 
of  middle  plane  and  the  contraction  of  thickness  cause  a 
second-order  (or  nonlinear)  frequency  shift  [26]. 

It  is  found  that  the  resonant  frequencies  are  affected 
by  zero-order  strains  through  the  second-  and  third-order 
elastic  stiffness  coefficient,  and  also  by  first  order  strain 
gradients  through  the  third-order  elastic  stiffness  [17]. 
In  other  words,  resonant  frequency  shift  will  be  induced 
by  the  mechanical  strain  or  stress.  For  a  given  flow  rate 
in  the  flow  path  of  the  chamber,  finite  element  analysis 
(FEA)  was  performed  using  ANSYS  (ANSYS  Inc.,  Can- 
onsburg,  PA),  the  results  indicate  that  the  magnitude 
of  AP  is  significantly  (around  3500  times)  higher  than 
rr.  Thus  AP  could  be  the  primary  factor  for  resonant 


Pressure 
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Fig.  5.  Deflection  of  a  quartz  resonator  disk  subjected  to  normal  pressure 
caused  by  fluid  flow:  (a)  side  view  and  (b)  top  view. 


frequency  shift  of  the  quartz  resonator,  and  the  effect  of 
Tr  on  frequency  change  is  then  ignored  in  the  following 
analysis. 

In  the  theoretical  analysis  of  force-frequency  relation¬ 
ship  for  quartz  resonators,  the  isotropic  assumption  is 
often  made  for  the  calculation  of  the  stress  distribution, 
and  the  results  were  found  to  be  consistent  with  the  ex¬ 
perimental  results  [17],  [27].  For  the  determination  of  the 
mechanical  deflection  of  quartz  resonators  subjected  to  a 
stead y-fluid  flow  rate,  it  is  assumed  that  the  quartz  reso¬ 
nator  is  isotropic  plate  with  Poisson’s  ratio  v,  Young’s 
modulus  E ,  and  thickness  hr.  Thus,  for  an  isotropic  plate 
with  clamped  edge,  the  deflection  w  due  to  the  uniform 
normal  load  AP  on  the  surface  can  be  obtained  by  [28]: 

w  =  —  (P2-r2)2,  (6) 

where  D  =  (Ehr)/(  12(1  —  v  2))  is  the  rigidity  of  the  quartz 
plate,  AP  is  the  normal  pressure  difference  across  the 
thickness  of  resonator,  R  is  radius,  and  r  is  the  radial  dis¬ 
tance  of  the  point  of  interest  in  the  middle  plane.  There¬ 
fore,  the  radius  of  the  middle  plane  will  be  elongated  from 
R  to  R  +  AP,  as  was  schematically  shown  in  Fig.  5(b), 
and  AP  is  given  by: 
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The  resonant  frequency  shift  can  be  determined  by  the 
force-induced  strain  distribution  [17],  and  the  main  terms 
for  fundamental  resonant  frequency  are  given  [26]: 
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where  the  coefficients  A2,  A2,  and  A3  are  determined  by 
the  second-order  and  third-order  elastic  constants;  ef,  £2, 
and  £3  are  the  initial  strains  of  middle  plain  due  to  the 
pressure.  Here,  for  estimation,  the  average  strains  are  used 
[26],  and  the  expression  can  be  easily  found  for  £°  and 
£3  : 
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Because  the  normal  stress  along  the  thickness  direction 
is  approximately  zero  for  a  very  thin  plate,  according  to 
strain  and  stress  relations,  £2  can  be  given: 


where  Ct]  is  the  stiffness  of  quartz  crystal. 

From  (1)  to  (10),  the  relationship  between  the  funda¬ 
mental  resonant  frequency  shift  and  flow  rate  is  found  as: 


Q2 

(11) 


It  can  be  seen  from  (4)  that  the  normal  pressure  differ¬ 
ence  is  proportional  to  the  flow  rate;  and  (8)-(10)  show 
that  the  resonant  frequency  shift  has  a  quadratic  rela¬ 
tionship  with  normal  pressure  difference.  Thus,  as  shown 
in  (11),  the  resonant  frequency  shift  of  quartz  resonator 
caused  by  fluid  flow  is  in  a  quadratic  relation  with  the  flow 
rate.  For  sensor  design  and  application,  it  is  often  desir¬ 
able  to  establish  a  linear  or  first-order  transduction  rela¬ 
tionship  between  the  measure  and  the  output  electrical 
signal.  With  frequency  as  the  output  signal,  a  nonlinear 
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Fig.  6.  Admittance  (Y  =  G  +  jB)  spectrum  near  fundamental  resonant 
frequency  in  the  air  and  water:  (a)  conductance  G  and  (b)  susceptance 

B. 


relationship  is  still  quite  convenient  for  sensor  applica¬ 
tions. 


III.  Results  and  Discussion 

The  electrical  admittance  spectra  of  TSM  quartz  reso¬ 
nator  are  first  measured  in  the  air  and  deionized  (DI) 
water  with  zero  flow  rate.  Fig.  6  shows  the  admittance- 
frequency  spectra  of  the  resonator  near  the  fundamental 
resonant  frequency.  G  is  conductance,  the  real  part  of  ad¬ 
mittance  (Y),  and  B  is  susceptance,  the  imaginary  part 
of  admittance.  It  can  be  seen  that  resonant  frequency  is 
decreased  due  to  the  liquid  load,  and  the  spectrum  is  wid¬ 
ened  and  the  peak  is  greatly  decreased  due  to  increased 
viscous  damping.  The  increased  damping  is  reflected  from 
the  change  of  motional  resistance  (R)  and  mechanical  Q  of 
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the  resonator.  The  motional  resistance  and  Q  factor  of  the 
resonator  are  found  to  be  29  H  and  39477  in  the  air,  and 
362  n  and  3077  in  the  water.  The  1st  resonance  frequency 
shift  due  to  water  load  is  around  645  Hz,  which  is  very 
close  to  the  theoretical  value  669  Hz  calculated  according 
to  the  Kanazawa  equation  [29] ,  [30] . 


=  fo'i/2(PLVL) 

(jrp  qpq) 1/2 


where  /q,  pi,  rji,  pq,  and  pq  are  1st  resonant  frequen¬ 
cy  in  the  air,  density  of  liquid,  viscosity  of  liquid,  den¬ 
sity  of  quartz,  and  piezoelectric  stiffened  shear  modu¬ 
lus;  pL  =  1000  (kgm-3),  r]L  =  8.7048  x  10-4(Pa-s), 
pq  =  2651(kgm-3),  and  pq  =  2.947  x  1010  (Nm“2)  were 
used  for  the  calculation.  The  resonant  frequency  of  reso¬ 
nator  in  the  DI  water  with  zero  flow  rate  is  used  as  the 
baseline  for  the  frequency  shift  of  device  under  different 
flow  rates. 

It  is  found  that  the  flow  variation  due  to  the  pulse 
pump  makes  the  resonator  admittance  spectrum  unstable 
under  a  certain  flow  rate,  which  induces  some  error  in  ex¬ 
tracting  resonant  frequencies.  Even  when  a  pulse  damper 
is  adopted,  the  effect  of  flow  variation  can  still  be  seen  in 
the  admittance  spectrum,  as  shown  in  Fig.  7(a).  To  reduce 
the  measurement  noise,  5-times  average  values  are  taken 
for  each  data  point  in  the  electrical  admittance-frequency 
spectrum,  which  is  an  average  measurement  algorithm 
provided  by  the  precision  impedance  analyzer.  From  Fig. 
7(b)  it  can  be  seen  the  admittance  spectrum  after  average 
becomes  stable  and  smooth,  which  improves  the  precision 
of  extracting  resonant  frequency. 

For  flow  measurement,  9  different  flow  rates  (0,  375, 
750,  1125,  1500,  1875,  2250,  2625,  and  3000  mL/min)  were 
chosen  in  the  acoustic  wave  flow  sensor  test  by  justifying 
the  pumping  speed  in  the  water  circulation  setup.  Fig.  8 
shows  the  frequency  shifts  for  fundamental  shear  mode 
and  higher  order  resonances  of  the  quartz  acoustic  wave 
resonator  under  these  flow  rates.  It  can  be  seen  that  reso¬ 
nant  frequency  increases  with  flow  rate.  In  addition,  the 
peak  value  of  conductance- frequency  spectrum  (G-f)  has 
a  different  trend-to-flow  rate  for  different  modes.  For  the 
fundamental  mode,  the  peak  decreases  for  the  first  8  flow 
rates  and  increases  at  the  9th  flow  rate.  For  the  3rd,  5th, 
and  7th  modes,  it  seems  that  the  G  peak  has  not  clear 
trend  with  the  flow  rate.  For  the  9th  mode,  the  peak  in¬ 
creases  with  the  flow  rate.  Compared  with  conductance 
peak  value,  it  is  clear  that  the  frequency  shift  is  an  indica¬ 
tion  of  flow  rate;  thus,  it  can  be  used  as  the  sensor  output 
for  the  acoustic  wave  flow  sensor. 

Fig.  9  plots  the  results  of  different  mode  resonant  fre¬ 
quency  shifts  versus  flow  rate.  It  is  found  that  the  reso¬ 
nant  frequency  shift  has  a  nearly  perfect  quadratic  rela¬ 
tion  with  flow  rate,  agreeing  very  well  with  the  prediction 
of  (11).  It  is  also  found  that  the  resonant  frequency  shifts 
of  all  the  overtones  are  basically  the  same  for  a  certain 
flow  rate  when  compared  in  terms  of  fractional  frequency 
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(b) 

Fig.  7.  Typical  electrical  admittance  (Y  =  G  +  jB)  spectra  of  quartz 
resonator  flow  sensor:  (a)  without  average  and  (b)  after  5-time  average. 
G  =  conductance;  B  =  susceptance. 


shift  (A f/f.  For  example,  the  fractional  frequency  shifts 
of  fundamental  and  overtones  (3rd,  5th,  7th,  and  9th)  for 
a  flow  rate  of  2250  mL/min  are  found  to  be  118,  127.6, 
127.7,  127.6,  and  127.2  ppm. 

The  repeatability  of  the  resonator  was  also  investigat¬ 
ed.  Nine  flow  rate  levels  were  subsequently  applied  from  0 
up  to  3000  mL/min  and  from  3000  back  to  0  mL/min.  For 
each  flow  rate,  the  duration  time  is  1  min,  and  the  reso¬ 
nant  frequencies  were  monitored  in  real  time  through  a 
Lab  VIEW  program  (National  Instruments).  Fig.  10  plots 
the  result  of  the  fundamental  shear  mode  resonant  fre¬ 
quency  shift  with  time.  It  can  be  seen  that  the  resonator 
has  an  excellent  repeatability  for  different  flow  rates.  Vari¬ 
ations  of  the  frequency  shift  with  time  are  also  observed, 
particularly  under  the  5th  and  6th  flow  rate  levels,  which 
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Fig.  8.  Frequency  shift  of  quartz  resonator  conductance  spectrum  under  different  flow  rates:  (a)  1st  mode,  (b)  3rd  mode,  (c)  5th  mode,  (d)  7th 
mode,  and  (e)  9th  mode.  Along  the  arrow  direction,  the  flow  rate  increases  from  0  to  3000  mL/min,  and  the  conductance  spectrum  moves  from  left 
to  right  accordingly. 
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Fig.  9.  Resonant  frequency  shifts  of  fundamental  mode  and  overtones 
versus  flow  rate. 


Time  (s) 

Fig.  10.  Repeatability  of  fundamental  resonant  frequency  shift  under  9 
different  flow  rates. 


may  reduce  the  output  resolution  of  the  acoustic  wave 
flow  sensor.  This  variation  is  mainly  attributed  to  the 
flow  rate  instability  in  the  flow  path  of  the  flow  chamber. 
Another  possible  reason  for  this  variation  is  the  decrease 
of  Q  due  to  liquid  load,  which  broadens  the  admittance 
spectrum  and  makes  it  more  difficult  to  identify  the  peak 
frequency  of  the  spectrum  because  the  resonant  frequency 
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is  extracted  through  searching  the  peak  of  conductance 
spectrum. 

The  experimental  results  indicate  some  attractive 
features  of  the  TSM  quartz  acoustic  wave  flow  sen¬ 
sor.  One  is  its  structural  simplicity;  no  coating  layer 
or  other  auxiliary  sensing  component  is  needed  in  the 
construction  of  the  flow  sensor  cell.  Another  feature  is 
that  the  operation  of  the  acoustic  wave  flow  sensor  is 
very  simple:  the  frequency  shift  of  the  quartz  resona¬ 
tor  is  a  convenient  output  that  can  be  easily  read  by  a 
simple  oscillator  circuit  and  frequency  counter.  In  ad¬ 
dition,  a  wide  range  of  flow  rates  can  be  measured  by 
using  the  quartz  acoustic  wave  flow  sensor.  Moreover, 
for  many  industrial  applications,  one  is  interested,  not 
in  the  resolution  of  flow  rate,  but  in  the  threshold  of 
flow  rate  or  pressure  for  monitoring  the  operational  pa¬ 
rameters  of  electromechanical  or  energy  systems.  Thus, 
compared  with  the  conventional  expensive  and  complex 
flow-sensing  devices,  the  acoustic  wave  flow  sensor  stud¬ 
ied  in  this  paper  may  provide  an  alternative  method  for 
flow-rate  monitoring  with  substantially  simple  sensor 
cell  structure  and  low  operation  cost. 

It  should  be  noticed  that,  before  the  practical  imple¬ 
mentation  of  this  acoustic  wave  flow  sensor,  some  issues 
need  to  be  considered.  Using  the  relationship  between  the 
flow  pressure  and  flow  rate  given  by  (4),  it  is  possible 
to  extract  the  flow  rate  by  the  measurement  of  pressure 
at  one  point  using  a  quartz  resonator  sensor.  However, 
the  flow  rate  measurement  is  also  related  to  the  outlet 
pressure  [second  term  in  (4)],  which  is  often  unknown  in 
practical  applications.  Thus,  to  eliminate  the  effect  of  the 
outlet  pressure,  one  approach  is  to  introduce  a  fluid  con¬ 
nection  channel  from  the  outlet  to  the  other  side  of  the 
quartz  resonator  such  that  the  other  side  of  the  device 
is  exposed  to  a  hydrostatic  pressure  that  is  equal  to  the 
pressure  at  the  outlet  port.  For  this  case,  (4)  and  (11)  will 
then  be  simplified  as: 
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Another  approach  is  to  adopt  a  dual  resonator  configu¬ 
ration:  if  2  quartz  TSM  resonators  are  installed  in  the  flow 
chamber  with  one  side  contacting  liquid  and  another  side 
exposed  to  air.  Assuming  they  are  located  at  (aq,  ft/2,  0) 
and  (aq,  ft/2,  0),  then,  according  to  (1)  and  (11),  one  can 
get: 
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where  P(xj),  A Pit  and  A/,//,  are  pressure  at  x  =  xt  in 
the  flow  channel,  pressure  difference  across  the  resonator 
thickness,  and  fractional  frequency  shift  of  the  resonator. 
One  can  find: 


P{ *2)  -  P(xi)  =  A  P2  -  A  Ph 


(16) 


Combining  (15)  and  (16): 
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From  (17),  it  can  be  found  that  by  using  the  dual  reso¬ 
nator  configuration,  flow  measurement  is  independent  of 
outlet  pressure  and  Q  can  be  obtained  by  the  frequency 
shift  of  2  resonators. 

The  sensitivity  of  the  TSM  quartz  acoustic  wave  flow 
sensor  seems  not  as  high  as  the  flow  sensors  based  on 
SAW  device,  which  has  a  frequency  shift  greater  than 
142  kHz  (1900  ppm)  for  gas  flow  rate  in  the  range  from 
0  to  1000  mL/min  [6].  The  experimental  results  in  this 
work  show  for  1000  mL/min  flow  rate  a  frequency  shift 
less  190  Hz  (40  ppm)  for  the  fundamental  shear  mode 
resonance  and  2000  Hz  (43  ppm)  for  the  9th  mode.  This  is 
because  a  nonlinear  (second-order  effect)  frequency  shift 
of  the  quartz  resonator  is  used  for  the  signal  transduction, 
which  is  usually  much  smaller  than  a  linear  effect.  The  9th 
overtone  has  a  frequency  shift  of  10300  Hz  (228  ppm)  for 
3000  mL/min  flow  rate,  which  is  equivalent  to  0.076  ppm/ 
(mL/min).  Assuming  that  the  quartz  resonator  sensor  has 
a  temperature  dependent  noise  level  of  1  ppm  on  the  fre¬ 
quency  shift,  the  calculated  resolution  for  the  flow  sensor 
is  about  13  mL/min. 

It  should  also  be  noted  that  theoretically  this  acoustic 
wave  flow  sensor  is  applicable  to  the  flow  rate  measure¬ 
ment  for  any  type  of  fluid  such  as  air,  oil,  and  others. 
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Other  conduit  materials  and  geometry  may  be  used  for 
the  fabrication  of  the  flow  sensor  structures  to  meet  the 
need  of  specific  applications,  as  long  as  the  fluid  pressure 
is  not  high  enough  to  cause  the  failure  of  thin  quartz  reso¬ 
nator  disk. 

In  addition,  because  the  resonant  behavior  can  be  af¬ 
fected  by  the  chemical  reactions  or  physical  adsorption 
occurring  on  the  resonator  electrode,  the  electrode  facing 
the  fluid  needs  to  be  protected,  particularly  when  it  is 
exposed  to  reactive  or  abrasive  fluid.  For  this  purpose,  in 
the  measurement  of  abrasive  fluids,  the  electrode  material 
should  be  anti-abrasive  or  a  protective  thin  film  needs  to 
be  coated  on  the  electrode  of  the  resonator  sensor.  Further 
study  is  also  needed  to  investigate  the  resonator  behavior 
under  higher  pressure  levels  and  determine  the  maximum 
pressure  that  the  resonator  can  endure  for  long-term  use, 
because  the  highest  pressure  in  the  experiment  study  is 
about  6  x  104  Pa,  which  is  quite  moderate  as  compared 
with  many  practical  applications.  Furthermore,  it  is  well 
known  that  the  flow  is  laminar  and  the  velocity  profile 
is  parabolic  for  low  Reynolds  numbers  (generally  under 
2000),  and  the  flow  becomes  turbulent  for  Reynolds  num¬ 
bers  over  4000.  The  Reynolds  numbers  of  the  current  ex¬ 
periments  spanned  from  0  to  822,  therefore,  experiments 
at  high  Reynolds  numbers  are  still  needed  to  investigate 
the  influence  of  the  turbulent  flow  profile  on  the  resonator 
response. 

IV.  Conclusion 

In  this  paper,  a  TSM  quartz  resonator  was  studied  as 
a  flow  sensor  to  measure  the  liquid  flow  rate.  A  special 
flow  sensor  structure  was  designed  and  constructed,  and 
the  theoretical  relationship  between  the  quartz  resonator 
frequency  shift  and  flow  rate  was  derived.  Both  theoretical 
and  experimental  results  have  indicated  that  the  frequen¬ 
cy  shift  of  quartz  resonator  caused  by  fluid  flow  can  be 
attributed  to  the  nonlinear  (quadratic)  frequency  response 
of  the  device  due  to  the  normal  pressure  imposing  on  the 
sensor  disk  through  the  fluid  flow.  The  resonant  frequen¬ 
cies  of  different  modes  can  be  used  as  the  sensor  output 
for  flow  rate  monitoring,  and  the  fundamental,  3rd,  5th, 
7th,  and  9th  mode  resonant  frequency  shifts  were  found 
to  be  around  920,  3572,  5947,  8228,  and  10300  Hz  for 
flow  rate  variations  from  0  to  3000  mL/min  or  a  Reynolds 
number  change  from  0  to  822.  To  obtain  the  flow  rate  by 
measuring  the  pressure  drop  across  the  quartz  resonator 
disk  using  a  single  resonator  sensor  configuration,  the  ef¬ 
fect  of  the  outlet  pressure  needs  to  be  eliminated.  It  is 
proposed  that  a  fluid  connection  channel  from  the  outlet 
to  the  other  side  of  the  quartz  resonator  can  be  introduced 
such  that  the  other  surface  of  the  device  is  exposed  to  a 
hydrostatic  pressure  that  is  equal  to  the  pressure  at  the 
outlet  port.  Another  approach  for  flow  rate  measurement 
is  to  use  dual-sensor  configuration,  in  which  2  quartz  reso¬ 
nators  are  installed  in  the  flow  chamber  along  the  flow 
channel.  The  relative  frequency  shifts  of  the  2  resonators 


will  be  directly  related  to  the  flow  rate,  and  independent 
of  the  outlet  pressure.  In  addition,  considering  that  the 
resonant  behavior  can  be  affected  by  the  chemical  reac¬ 
tions  or  physical  adsorption  occurring  on  the  resonator 
electrode,  the  electrode  facing  the  fluid  may  need  to  be 
protected,  particularly  when  it  is  exposed  to  reactive  or 
abrasive  fluids.  Further  experiment  studies  under  higher 
pressure  and  high  Reynolds  numbers  are  needed  to  deter¬ 
mine  the  pressure  limit  and  the  flow  profile  influence  on 
resonator  response. 
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Abstract — We  study  the  theory  of  the  basic  vibration  char¬ 
acteristics  of  a  circular  cylindrical  shell  piezoelectric  trans¬ 
ducer.  The  linear  theory  of  piezoelectricity  is  used.  Both  the 
free-vibration  solution  for  resonant  frequencies  and  modes  as 
well  as  the  electrically  forced-vibration  solution  for  admittance 
are  obtained.  Numerical  results  are  presented. 


I.  Introduction 

PIEZOELECTRIC  materials  are  frequently  used  for  mak¬ 
ing  electromechanical  transducers  of  various  shapes. 
We  are  interested  in  circular  cylindrical  shell  piezoelec¬ 
tric  transducers  to  be  attached  to  circular  cylindrical 
structures  for  exciting  and  detecting  acoustic  waves.  It  is 
part  of  our  effort  in  the  development  of  electromechanical 
transducers  for  nondestructive  evaluation  (NDE)  of  circu¬ 
lar  cylindrical  elastic  structures,  in  particular  the  circular 
cylindrical  pressure  vessels  in  nuclear  reactors.  Existing 
theoretical  results  on  circular  cylindrical  piezoelectric 
transducers  are  mostly  for  full  cylinders,  e.g.  [1]—  [14] ,  but 
the  application  we  are  interested  in  requires  knowledge  of 
vibration  characteristics  of  a  finite  portion  of  a  circular 
cylinder.  There  are  limited  results  on  vibrations  of  finite 
circular  cylindrical  piezoelectric  panels  [15],  [16]  where 
mixed  boundary  conditions  are  chosen  in  a  particular  way 
so  that  analytical  solutions  are  possible,  but  the  mixed 
mechanical  boundary  conditions  treated  are  different  from 
typical  situations  in  our  applications.  There  are  also  ex¬ 
tensive  results  on  circular  cylindrical  piezoelectric  shells 
for  smart  structural  applications  [17]  where  approximate, 
2-D  structural  theories  are  used.  Results  from  smart  struc¬ 
tural  analyses  are  for  the  so-called  low-frequency  modes  of 
extension  and  flexure,  etc.  They  are  insufficient  for  NDE 
where  the  so-called  high-frequency  thickness  modes  [18] 
are  usually  used,  for  which  analyses  based  on  the  exact 
3-D  equations  are  necessary.  In  this  paper,  we  formulate 
a  problem  of  finite  circular  cylindrical  shell  piezoelectric 
transducers  that  is  relevant  to  NDE  of  circular  cylindri¬ 
cal  structures  and  perform  a  theoretical  analysis  on  these 
transducers. 
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II.  Governing  Equations 

Consider  a  portion  of  a  circular  cylindrical  shell  with 
inner  radius  a  and  outer  radius  b  as  shown  in  Fig.  1.  The 
2  radii  going  through  the  edges  of  the  shell  form  an  angle 
2(3;  (3  is  less  than  tt  and  is  otherwise  arbitrary.  The  shell 
is  unbounded  in  the  2  direction  determined  by  the  right- 
hand  rule  from  the  x  and  y  axes  in  the  figure.  Only  a  cross 
section  of  the  shell  is  shown.  The  material  of  the  shell  is 
ceramics  poled  in  the  z  direction.  A  cylindrical  coordinate 
system  is  defined  by  x  =  rcosd,  y  =  rsinf?,  and  z  =  z.  In 
the  index  notation  below,  ( x,y,z )  correspond  to  (1,2,3). 

We  consider  the  so-called  antiplane  or  shear  horizontal 
(SH)  vibration  with  [19],  [20] 

m  =  u2  =  0,  u 3  =  u(xhx2,t),  (j)  =  (j){xhx2,t),  (1) 


where  Ui  is  the  displacement  vector  and  ip  is  the  electric 
potential.  Modes  described  by  (1)  are  allowed  by  the  3-D 
equations  of  linear  piezoelectricity.  For  these  modes,  the 
3-D  equations  of  linear  piezoelectricity  reduce  to  a  2-D 
mathematical  problem.  Corresponding  to  the  displace¬ 
ment  and  potential  fields  in  (1),  the  nonvanishing  strain 
and  electric  field  components  are 


25 13 
25  23 


=  V«, 


Ei 
E  2 


-V0, 


(2) 


where  V  =  ii9i  +  F92  is  the  2-D  gradient  operator.  For 
ceramics  poled  in  the  23  direction,  the  nontrivial  compo¬ 
nents  of  Tij  and  D%  are 


Tn 

T23 


cVa  +  eV<^, 


eVu  -  eV^>,  (3) 


where  we  have  denoted  the  relevant  elastic,  piezoelectric, 
and  dielectric  constants  by  c  =  C44,  e  =  ei5,  and  e  =  en. 
The  nontrivial  equation  of  motion  and  the  charge  equa¬ 
tion  of  electrostatics  take  the  following  form: 


cV  2u  +  eV  2^>  =  pii,  eV  2u  —  eV  2(j)  =  0.  (4) 


We  introduce  [19],  [20] 


=  <t> 


e 

-u. 

£ 


(5) 
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Fig.  1.  A  circular  cylindrical  shell  piezoelectric  transducer. 


Substitution  of  (10)  into  (9)  results  in 


d2u  1  du 
dr 2  r  dr 


d  2ip  1  dip  v2 
dr2  r  dr  r2 


V 

—  M  =  0, 
r  ) 

(11a) 

ip  =  o, 

(lib) 

where  we  have  denoted 


£  =  — • 
Vt 


(12) 


Eq.  (11a)  can  be  written  as  Bessel’s  equations  of  order 
v;  (lib)  allows  a  simple  power  function  solution.  Then, 
when  v  g  0, 


Then,  in  terms  of  u  and  ip, 


u  =  [CJ  „g  r)  +  DF„g  r)]  cos  v0, 
ip  =  [Fr  v  +  Gr  ~v]  cos  v0. 


(13a) 


T23  =  cu  2  +  etjj p,  T31  =  cm. i  +  eip}h 

Hi  =  -ei/jj,  H2  =  -£^,2, 


and 

VtV2u  =  ii,  V2i/j  =  0, 


(6)  where  Jv  and  F^  are  the  /Th  order  Bessel  functions  of  the 
first  and  second  kind.  C,  H,  F.  and  G  are  undetermined 
constants.  When  v  =  0,  we  have 


(7) 


m  =  G/0gr)  +  DF0gr), 
ip  =  F  In  r  +  G. 


(13b) 


where 


The  other  fields  are  then  given  by  [v  ^  0): 


Vt 


C  _  G  /  T  T  o  G 

— ,  c  =  c  -\ - =  c(l  +  k2),  k2  =  — . 

p  £  £  C 


In  polar  coordinates,  (7)  takes  the  form 
d2u  1  du  1  d2~ 


9 

vt 


u 

,  2  nq  2 
)  2. 


<9r2  rdr  r  2  de 

d2ip  1  dip  1  d2ip 
dr2  +  r  dr  +  r2  dd: 


=  u, 


=  0. 


(8) 


(9) 


Mechanically  the  boundaries  are  all  traction-free.  Electri¬ 
cally  the  boundaries  at  9  =  ±/3  are  unelectroded.  The 
boundaries  at  r  =  a,  b  may  be  either  electroded  or  unelec¬ 
troded.  We  will  consider  time-harmonic  modes  and  use 
the  complex  notation.  All  fields  are  with  an  exp (iot)  fac¬ 
tor  that  will  be  dropped  for  convenience. 


=  ip  H — M 
£ 

-CJ „g r)  +  -DYV  g  r)  +  Frv  +  Gr~v  cos v9, 

£  £ 

T 'rz  =  CU  r  T  61p  r 

=  \ciCJ'v  gr)  +  ct,DY'v  gr) 

+  evFr 1,-1  —  evGr  _!'~1]cosiz0, 

T6z  =  c -  u  o  +  e-ip  g 
r  r 

=  —  —  [cCJ  „g  r)  +  cD  Yv  g  r)  +  eFr  v  +  eGr  ~v\  sin  v9 , 
r 

£>r  =  —  £ip,r  =  [— £vFrv~1  +  £vGr  “y_1]cosiz0, 

Hg  =  —£-ipe  =  £  —  [Frv  +  Gr~v]smv9, 
r  ’  v 


(14a) 


or,  when  iz  =  0: 


III.  Free  Vibration 

According  to  whether  m  and  ip  are  even  or  odd  functions 
of  9 ,  we  classify  the  modes  into  symmetric  and  antisym¬ 
metric.  First,  consider  the  following  symmetric  modes: 

u  =  u(r)  cos  v9,  ip  =  ip(r)  cos  v9.  (10) 


<p  =  ip  +  -u  =  -[CJ0gr)  +  DF0gr)]  +  Fhir  +  G, 

£  £ 

F 

Trz  =  cu  r  +  eip  r  =  —cQCJ igr)  +  DY^r)]  +  e  — , 

r 

F 

D  r  =  — £ljj  r  =  —£  —  . 

r 

(14b) 
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For  traction-free  and  unelectroded  edges  at  9  =  ±/3,  we 
must  have  T$z  =  0  and  De  =  0,  which  implied  that 


_1  -vma  " ~  1 


=  0. 


(20b) 


u,e  =  0,  =  0,  0  =  ±p.  (15) 

From  (15)  we  obtain 

TY17T 

sin  v[3  =  0,  v(3  =  rmr ,  v  =  vm  =  — — ,  m  =  0, 1, 2,  •••. 

P 

(16) 


Eq.  (20b)  has  no  solution  because  m  /  0.  In  the  case  of 
(19b),  we  have  F  =  0,  which  implies  a  constant  i/).  In  prac¬ 
tice,  we  are  more  interested  in  the  case  of  an  electroded 
plate  below. 

B.  Electroded  Transducer 


Similarly,  for  antisymmetric  modes,  we  have 


cos  vfj  =  0, 


v  =  v 


m 


7 T 

P 


7T) 

m  =  0,1,2,  •••. 


(17) 


For  antisymmetric  modes  v  ^  0.  We  will  discuss  unelec¬ 
troded  and  electroded  transducers  separately  below. 

A.  Unelectroded  Transducer 


For  traction-free  and  unelectroded  boundaries  at  r  = 
a,b ,  we  must  have  Trz  —  0  and  Dr  =  0,  or,  equivalently, 


When  the  transducer  is  electroded  at  r  =  a,  b  and  the 
electrodes  are  shorted,  the  electrical  boundary  conditions 
at  r  =  a,  b  are  that  ip  =  0.  When  m  /  0,  we  have 

e-CJUm{£,a)  +  -DYVjn(t;  a)  +  Fav™  +  Ga~v "  =  0, 

c&J'vSta)  +  ctDY'vSta)  + 

evmFaVm~ 1  —  ev  mGa~Vm~1  =  0, 

- CJVm{£b )  +  -  D  YUJ£  b)  +  Fb"™+  Gb~v »  =  0, 
£cCCJ',mm  +  c£DYlm(tb)  + 

evmFbl'm~1  —  ev  mGb~Vm~l  =  0. 

(21a) 


u>r  =  0,  i/’.r  =  0,  r  =  a,b. 


(18)  When  m  =  0  we  have 


Then,  when  v  ^  0  or  m  ^  0,  we  have 

C£J'vm(£  a)  +  D£Y[,J£  a)  =  0, 
FvmaVm~1  -  Gv  =  0, 

+  D&»mm  =  o, 

Fvmb^-1-Gvmb-'«-1  =  0, 


(19a) 


or,  when  m  =  0: 


-c£[CJ  i(£a)  +  DYfta)]  +  e- 

a 

— c£[GJ  i(£b)  +  DY,m  +  ey 

b 

F 

— E  — 

a 

F 


=  0, 
=  0, 
=  0, 
=  0. 


(19b) 


For  nontrivial  solutions  of  G,  D,  F,  or  G,  the  determinant 
of  the  coefficient  matrix  of  (19a)  or  (19b)  must  vanish, 
which  gives  the  frequency  equation  that  determines  the 
resonant  frequencies.  We  note  that  lines  1  and  3  in  (19a) 
are  not  coupled  to  lines  2  and  4  in  (19a).  Therefore  the 
frequency  equation  determined  by  (19a)  can  be  factored 
into  2 


£4m(£a)  £nM(£a) 


(20a) 


-cQCJita)  +  DYUa)}  +  e-  =  0, 

a 

-c^CJ^b)  +  DY^b)}  +  ey  =  0, 

b  (21b) 

-[CJ  0(£a)  +  DY0^a)}  +  Fin  a  +  G  =  0, 

£ 

-  [CJ  o(£&)  +  DY0^b)\  +  Flub  +  G  =  0. 
e 

The  frequency  equation  is  obtained  by  setting  the  deter¬ 
minant  of  the  coefficient  matrix  of  (21a)  or  (21b)  to  zero. 

C.  Numerical  Results 

As  a  numerical  example,  for  geometric  parameters  we 
choose  a  =  0.2  m,  b  =  0.21  m,  (3  =  tt/6.  For  the  material, 
we  consider  PZT-7A  [21]  with  p  =  7600  kg/m3,  C44  =  2.53 
x  1010  N/m2,  ei5  =  9.2  C/m2,  en  =  0.407  x  10"8  C/ 
Vm. 

Fig.  2  shows  resonant  frequencies  versus  m;  m  is  effec¬ 
tively  the  wave  number  in  the  6  direction.  Therefore,  in 
Fig.  2,  if  the  dots  are  properly  connected  in  the  horizon¬ 
tal  direction,  the  resulting  curves  resemble  the  dispersion 
curves  of  SH  waves  in  plates  or  circular  cylindrical  shells, 
but  here  the  wave  number  and  frequencies  are  discrete. 
Frequencies  with  m  =  0  are  called  thickness-shear  frequen¬ 
cies  whose  modes  have  no  0  dependence.  When  m  is  posi¬ 
tive,  the  frequencies  are  called  thickness-twist  frequencies 
whose  modes  have  9  dependence. 
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Fig.  2.  Resonant  frequencies. 


(b) 

Fig.  3.  (a)  Displacement  and  (b)  electric  potential  for  m  =  0,  uj  =  0.60 
x  106  Hz. 


(a) 


(b) 

Fig.  4.  (a)  Displacement  and  (b)  electric  potential  for  m  =  0,  u>  =  1.55 
x  106  Hz. 


Figs.  3-8  show  displacement  and  electric  potential  dis¬ 
tributions  of  the  modes  whose  frequencies  are  near  the 
origin  of  Fig.  2.  Only  symmetric  modes  with  a  cos#  de¬ 
pendence  are  considered.  Antisymmetric  modes  are  of  less 
interest  because  they  are  not  electrically  excitable  by  a 
radial  electric  field,  which  will  be  seen  later  in  the  forced- 
vibration  analysis,  and  therefore  are  not  plotted.  Due  to 
symmetry,  only  the  fields  in  the  first  quadrant  of  Fig.  1 
with  0  <  9  <  f3  are  plotted.  The  figures  show  that  increas¬ 
ing  m  introduces  more  nodal  points  in  the  6  direction.  For 
a  fixed  m,  modes  with  more  nodal  points  in  the  r  direction 
have  higher  frequencies. 

IV.  Forced  Vibration 

For  forced  vibration,  we  consider  a  time-harmonic  volt¬ 
age  2  V  applied  across  the  electrodes  at  r  =  a,  b  with  the 
electric  potential  as 
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(a) 


(b) 

Fig.  5.  (a)  Displacement  and  (b)  electric  potential  for  m  =  1,  ui  =  0.05 
x  106  Hz. 


(a) 


0.22  0 

(b) 

Fig.  6.  (a)  Displacement  and  (b)  electric  potential  for  m  =  1,  a)  =  0.60 
x  106  Hz. 


—V  exp  (iwt), 
V  exp  (iwt), 


(22) 


Obviously  only  the  symmetric  modes  can  be  excited. 
Therefore,  a  general  symmetric  solution  satisfying  the 
transducer  edge  conditions  at  9  =  ±/3  is  constructed  as 

u  =  CoJ&r)  +  D0Y0(£r) 

OO 

+Y  [CmJvJ£ r)  +  DmYvj£ r)]cosi/m9, 

m= 1 

oo 

i/j  =  F0lnr  +  G0  +  Y  [ Fm.r Vm  +  Gmr ~i'm] cos v m0. 

m—1 

(23) 


The  boundary  conditions  at  r  =  a,  b  require  that  the 
conditions  in  (25)  (see  page  2333)  are  met. 

We  multiply  (25)  by  cos  vp6  and  integrate  the  resulting 
equation  from  -(3  to  (3  for  p  =  0,  1,  2,  3,  ....  Then,  from 
the  orthogonality  of  cos  vp91  the  summation  in  (25)  disap¬ 
pears.  Only  the  0-independent  terms  are  left  because  the 
applied  voltage  is  a  constant  and  does  not  vary  with  9.  To 
calculate  the  charge  on  and  the  current  flowing  into  the 
electrode  at  r  =  b,  we  need 

Dr  e  1525  rz  T  c r, 

Q  =  f&  - Dr\r=bbd9 ,  (26) 

/  =  Q  =  iuQ. 


The  stress,  electric  displacement,  and  electric  potential  Then  the  impedance  of  the  transducer  is  given  by 
needed  for  boundary  conditions  are  obtained  as  shown  in 

(24)  (see  next  page).  Z  =  2V/I.  (27) 
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i/j  H — u 
e 


=  -poJofcr)  +  D0Y0^r)}  +  F0\nr  +  G0 


(24a) 


-GmJVm{^r)  +  -DmYvJ£r)  +  Fmrv »  +  Gmr~v" 

£  £ 


cos  v  m9, 


Trz  =  CU,2  +  eip:  2 

=  -c£[CVi(£r)  +  D0Y^r)}  +  e ^ 


(24b) 


+  X]  ic€C™J'vj£r)  +  cZDmY'Vm(€r)  +  evmFmrVm  1-evmGmr  Vm  l]cosv m0, 


Dr  =  -£1p,r  =  -£—  +  ^  [~£UmFmrVm~1  + 


(24c) 


(b) 

Fig.  7.  (a)  Displacement  and  (b)  electric  potential  for  m  =  2,  uj  =  0.11 
X  106  Hz. 


(b) 

Fig.  8.  (a)  Displacement  and  (b)  electric  potential  for  m  =  2,  io  =  0.61 
x  106  Hz. 
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0 
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m= 1 


(25) 


With  (24c),  the  impedance  can  be  written  in  the  follow¬ 
ing  form: 


Z  = 


V 

iu)eF0p 


(28) 


In  the  numerical  calculation  for  the  forced  vibration,  the 
elastic  constant  C44  is  replaced  by  C44(l  +  iQ)  to  take 
viscous  damping  into  consideration.  For  ceramics,  Q  is  of 
the  order  of  10~3  to  10”2.  Q  =  0.05  is  used  in  the  calcula¬ 
tion. 

Fig.  9  shows  the  real  part  of  the  impedance  as  a  func¬ 
tion  of  the  driving  frequency.  Only  the  modes  with  m 
=  0  can  be  excited  due  to  the  uniformity  of  the  applied 
voltage.  Among  the  modes  with  m  —  0,  in  the  r  direction, 
odd  modes  are  mainly  antisymmetric  about  the  middle  of 


(a) 


0.22  0 


(b) 

Fig.  10.  (a)  Displacement  and  (b)  electric  potential  when  10  =  0.774  X 
106  Hz. 
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0.22  0 

(b) 


Fig.  11.  (a)  Displacement  and  (b)  electric  potential  when  u  =  2.32  x 
106  Hz. 

the  transducer  where  r  =  (a  - b  b)/ 2  and  even  modes  are 
mainly  symmetric.  Therefore,  the  major  peaks  in  Fig.  9 
correspond  to  the  first  and  the  third  resonant  frequencies 
with  m  =  0  in  Fig.  2. 

Displacement  and  electric  potential  distributions  at  the 
first  2  major  resonances  in  Fig.  9  are  shown  in  Fig.  10  and 
11.  What  is  shown  are  the  real  parts  of  the  complex  fields. 
These  are  modes  with  m  =  0.  Modes  with  higher  frequen¬ 
cies  have  more  variations  along  r. 

V.  Conclusion 

Exact  solutions  from  the  3-D  equations  of  linear  piezo¬ 
electricity  are  obtained  for  free  and  electrically  forced  vi¬ 
brations  of  a  finite  circular  cylindrical  panel  piezoelectric 
transducer.  The  modes  obtained  include  high-frequency 
modes  of  thickness-shear  and  thickness-twist.  These  cy¬ 
lindrical  transducers  can  be  used  to  excite  SH  modes  and 
are  useful  for  NDE  of  circular  cylindrical  structures.  When 
they  are  attached  to  such  structures,  the  modeling  prob¬ 


lem  becomes  much  more  complicated.  Numerical  proce¬ 
dures  that  may  be  needed  are  left  for  future  work. 
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AIN  thin  films  with  c-axis  orientation  have  been  investigated  for  fabricating  thin  him  bulk  acoustic 
wave  resonators  in  the  past  few  years.  Characterization  of  thin  him  material  properties  including 
density,  elastic  modulus,  and  piezoelectric  coefficient  is  essential  in  processing  study  and  for 
predicting  the  performance  of  the  acoustic  devices.  In  this  paper,  we  present  our  results  on  the 
fabrication  of  highly  c-axis  oriented  AIN  thin  hlms  on  Pt/Ti/Si  (100)  substrates  by  dc  reactive 
magnetron  sputtering  method.  The  crystalline  structure  and  the  surface  morphology  of  AIN  hlms  are 
characterized  by  x-ray  diffraction  and  scanning  electron  microscopy.  The  effective  piezoelectric 
coefficient  r/33e ff  of  the  AIN  hlms  was  measured  by  a  laser  interferometer  method  and  the 
piezoelectric  coefficient  c/33  was  estimated.  A  recently  developed  resonance  spectrum  method  is 
applied  to  characterize  the  electromechanical  properties  of  AIN  thin  hlms  based  on  the  input 
electrical  impedance  equation  derived  by  one-dimensional  transmission  line  theory  for  composite 
resonators.  Using  the  experimental  impedance  spectrum  data,  the  density  and  elastic  constant  of  the 
piezoelectric  AIN  thin  him  in  the  four-layer  composite  resonator  structure  are  evaluated.  The 
calculated  results  reveal  that  the  piezoelectric  coefficient  r/33,  density,  and  velocity  of  the  c-axis 
oriented  AIN  thin  him  are  4.19  pm/V,  3187.3  kg/m3,  and  10631  m/s,  respectively.  ©  2007 
American  Institute  of  Physics.  [DOI:  10.1063/1.2716391] 


I.  INTRODUCTION 

With  a  signihcant  progress  in  thin  him  deposition  tech¬ 
nologies  for  complex  material  systems  such  as  AIN,  ZnO, 
and  PZT  in  recent  years,  thin  him  bulk  acoustic  wave  reso¬ 
nator  (FBAR)  and  hlter  concepts  are  becoming  more  and 
more  attractive  for  microwave  frequency  control  applica¬ 
tions.  For  resonators  operating  in  the  gigahertz  range,  a  thin 
him  piezoelectric  layer  in  the  order  of  a  few  microns  with 
desirable  electromechanical  properties  is  required.  Among 
these  materials,  AIN  is  particularly  attractive  due  to  that  it 
has  a  number  of  interesting  properties  such  as  high  thermal 
conductivity,  high  electrical  insulation,  and  chemical 
stability. 1  In  addition,  it  is  piezoelectric  in  its  crystalline  form 
and  has  high  longitudinal  acoustic  velocity,  ~11  000  m/s. 
These  characteristics  make  it  possible  to  design  and  fabricate 
high  frequency  resonators  and  bandpass  hlters  for  signal  pro¬ 
cessing  and  communication  devices.  It  is  obvious  that  the 
deposition  of  high  quality  AIN  thin  him  with  desirable  ma¬ 
terial  properties  is  needed  for  the  fabrication  of  thin  him  bulk 
acoustic  wave  resonators  and  bandpass  hlters. 

In  this  paper,  we  report  our  recent  studies  on  the  depo¬ 
sition  of  highly  c-axis  oriented  AIN  thin  hlms  on  both  sap¬ 
phire  and  silicon  substrates  by  dc  reactive  magnetron  sput¬ 
tering  and  the  characterization  of  material  elastic  and 
piezoelectric  properties.  A  single  beam  laser  interferometer 
system  that  can  detect  very  small  displacement  signals  with 
resolution  down  to  10-3  A  in  the  low  frequency  range  was 
hrst  employed  to  measure  the  effective  piezoelectric  coeffi¬ 
cient  of  the  AIN  thin  him  through  converse  piezoelectric 
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effect,”  and  a  resonant  spectrum  method  was  applied  to  ex¬ 
tract  the  other  properties  including  density,  elastic  constant, 
and  acoustic  wave  velocity  of  the  AIN  hlms  based  on  the 
input  electric  impedance  of  thin  him  composite  resonators. 

Many  approaches  have  been  used  to  derive  the  input 
electric  impedance  of  a  thin  him  composite  resonator.  These 
include  the  equivalent  circuit  method  ’  and  transfer  matrix 
approach.  Sittig4  and  Meitzler  and  Sittig  used  the  transfer 
matrix  approach  to  calculate  the  impedance  of  the  multilayer 
piezoelectric  transducers  used  in  ultrasonic  digital  delay 
lines.  Zhang  et  al.  utilized  the  input  electric  impedance  spec¬ 
trum  for  a  four-layer  composite  resonator  to  characterize  the 
material  properties  of  ZnO  thin  him  on  the  quartz  substrate.6 
Here,  we  follow  a  similar  approach  by  applying  the  reso¬ 
nance  spectrum  method  to  characterize  the  material  proper¬ 
ties  of  the  AIN  thin  him  in  a  composite  resonator. 

II.  METHODS  FOR  PROPERTY  CHARACTERIZATION 

A.  Measurement  of  effective  piezoelectric 
coefficient  d33eff  of  AIN  thin  film 

In  order  to  measure  the  piezoelectric  coefficient  of  the 
AIN  thin  him,  both  the  direct  and  converse  piezoelectric  ef- 
fects  can  be  employed.”’  It  is  in  general  more  convenient 
to  use  the  converse  piezoelectric  effect  to  characterize  the 
piezoelectric  coefficient  of  thin  him  or  bulk  materials  by 
measuring  the  displacement  of  the  surface  of  the  materials 
subjected  to  the  external  ac  electric  held  or  voltage.  Since  the 
thin  him  AIN  is  deposited  on  the  substrate,  which  can  be 
regarded  as  a  thin  layer  clamped  on  the  substrate,  the  strain 
and  electric  held  ratio  d33eff=x3/E3  =  Ad/V  where  A d  is  the 
displacement  induced  by  applied  voltage  V,  does  not  repre¬ 
sent  the  real  piezoelectric  coefficient  d33  of  the  thin  him  in 
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free  mechanical  boundary  conditions  but  an  effective  value. 
The  piezoelectric  coefficient  of  AIN  is  relatively  small,  the 
resulting  deformation  excited  by  the  electric  field  is  corre¬ 
spondingly  very  small.  Thus  how  to  measure  the  small  de¬ 
formation  signal  of  the  AIN  accurately  is  important  for  the 
measurement  of  the  7/33^.  For  the  characterization  of  the 
piezoelectric  property  of  AIN  film,  a  very  sensitive  single 
beam  laser  interferometer  system  with  displacement  resolu¬ 
tion  down  to  10-3  A  in  the  low  frequency  range  was  em- 
ployed  in  the  measurement." 

B.  Resonance  spectrum  method 

In  the  past  decades,  several  methods  have  been  devel¬ 
oped  to  determine  the  coefficients  of  piezoelectric  materials 
in  order  to  better  characterize  and  describe  the  electrome¬ 
chanical  properties  and  performance  of  materials  and  de¬ 
vices.  The  measurement  methods  for  low  loss  materials  have 
been  summarized  in  IEEE  standards.10  Those  methods  have 
been  successfully  used  to  determine  the  coefficients  of  ma¬ 
terials  with  high  mechanical  quality  factors  such  as  lead  ti- 
tanate  zirconate  ceramics.11'1"  In  order  to  characterize  the 
material  properties  of  the  piezoelectric  composite  resonators, 
the  resonance  spectrum  method  has  also  been  developed, 4-0 
which  includes  the  effect  of  electrodes  and  substrate.  Using 
the  information  contained  in  the  impedance  or  admittance 
spectrum  for  piezoelectric  resonators  such  as  parallel  and 
series  resonance  frequencies,  the  elastic  and  electromechani¬ 
cal  properties  of  the  piezoelectric  thin  film  layer  can  be  ob¬ 
tained.  The  resonance  spectrum  method  is  particularly  suit¬ 
able  for  the  characterization  of  material  properties  of 
piezoelectric  thin  films  that  are  used  for  microdevice  appli¬ 
cations.  We  have  also  applied  the  transfer  matrix  approach  to 
derive  the  one-dimensional  input  electric  impedance  for 
three-layer  and  four-layer  thin  film  composite  resonators  and 
the  input  electric  impedance  has  been  used  for  the  discussion 
of  the  effect  of  electrodes  and  the  temperature  compensation 
in  composite  resonators.6'1  44  The  input  electric  impedance 
of  the  four-layer  resonator  is 

(zi  +  £2) sin  7  +  ;2(1  ~  cos  7) 

7  (z;  +z2)cos  7  +  7'(l  +  Z|Z2)sin  yj  ’ 

(1) 

where  Z\  =Zt /Z0  and  z2=Z2/Z0  are  the  normalized  acoustic 
loading  impedances  on  both  sides  of  the  piezoelectric  layer. 
Cf)= / 1  is  the  clamped  capacitor  of  the  resonator  with 
area  S  and  dielectric  permittivity  S33.  Z0=Spv  is  the  acoustic 
impedance  of  the  piezoelectric  layer  with  density  p,  v  is  the 
longitudinal  acoustic  wave  velocity  in  the  piezoelectric  layer 
along  the  direction  normal  to  the  resonator  surface,  y 
=  a>tlv  is  the  phase  delay  of  the  acoustic  wave  in  the  piezo¬ 
electric  film  in  which  o)=27r/.  It  is  noted  that  a  three-layer 
thin  film  resonator  is  a  special  case  of  four-layer  resonators 
with  zero  substrate  thickness.  Using  Eq.  (1),  the  electric  im¬ 
pedance  spectrum  of  the  thin  film  composite  resonator  can  be 
simulated  by  using  appropriate  material  property  data. 

It  is  known  that  parallel  resonances  correspond  to  maxi- 
mums  of  the  resistance  of  the  composite  resonator.  In  a  loss¬ 
less  resonator,  parallel  resonance  corresponds  to  the  infinite 


_  V _ 1_ 

/  ja>C0 


impedance  Zin.  Thus,  by  setting  the  denominator  equal  to 
zero,  we  can  get  the  determination  equation  for  parallel  reso¬ 
nance  frequencies:6 

(zi  +  z2)cos  y+  /(l  +  ZtZ2)sin  y=0.  (2) 


At  the  center  of  the  first  normal  region,  where  y~  77,  we 
can  get  the  space  between  two  adjacent  parallel  resonant  fre¬ 
quencies  which  can  be  derived  from  Eq.  (2): 


A/jv  =  A/o  1  + 


Pe  \te\  +  Pe  2fe2  +  Pf 
Pst/sb 


(3) 


where  Af0=vsh/(2tsh),  and  the  subscripts  el,  e2,  and  sb  rep¬ 
resent  the  top  electrode,  bottom  electrode,  and  substrate. 
Similarly,  at  the  center  of  the  first  transition  region,  where 
y~  7t/2,  we  can  get  the  space  between  two  adjacent  parallel 
resonant  frequencies:6 


a  r  A  r  I  ,  PslTsl/  PstPst/e2  PsbPe  I  t’sbfi' ! 

A/ 7  -  A/ol  1  +  t  +  1  + 


pv\b 


Pel^el^b 


2  2* 

P  V  tsb 


-1 


(4) 


The  density  and  the  acoustic  velocity  of  the  AIN  thin 
film  can  be  calculated  from  Eqs.  (3)  and  (4).  The  value  of  the 
stiffened  elastic  constant  of  the  film  can  then  be  calculated 
by 


cD  =  v2p. 


(5) 


C.  Nanomechanical  property  characterization 
by  nanoindentation 


Indentation  has  been  widely  used  to  measure  the  me¬ 
chanical  properties  of  thin  film  materials  and  a  host  of 
nanoindentation  methods  is  available,  which  employ  various 
tip  shapes  and  rely  on  different  model  assumptions. 16-1 8  The 
two  mechanical  properties  measured  most  frequently  using 
indentation  techniques  are  the  hardness  H  and  the  elastic 
modulus  E.  As  the  indenter  is  pressed  into  the  sample,  both 
elastic  and  plastic  deformation  occurs,  which  results  in  the 
formation  of  a  hardness  impression  conforming  to  the  shape 
of  the  indenter.  During  indenter  withdrawal,  only  the  elastic 
portion  of  the  displacement  is  recovered,  which  facilitates 
the  use  of  an  elastic  solution  in  modeling  the  contact  process. 
Thus  in  a  typical  nanoindentation  test,  data  from  a  complete 
loading  and  unloading  cycle  are  collected  to  generate  a 
force-distance  curve,  which  is  then  used  for  the  estimation  of 
nanomechanical  properties.  Nanoindentation  hardness  is  de¬ 
fined  as  the  indentation  load  divided  by  the  projected  contact 
area  of  the  indentation.  It  is  the  mean  pressure  that  a  material 
can  support  under  load.  From  the  load-displacement  curve, 
hardness  can  be  obtained  at  the  peak  load,  H=Pmax/A,  where 
A  is  the  projected  contact  area.  The  elastic  modulus  of  the 
indented  sample  can  be  obtained  from  the  initial  slope  of  the 
unloading  curve,  S=dP/dh.  A  geometry-independent  rela¬ 
tion  has  been  derived  for  contact  stiffness,  contact  area,  and 
elastic  modulus  as  follows:  19~_1 


S'  =  2/3 


(6) 


where  [3  is  a  constant  that  depends  on  the  geometry  of  the 
indenter  (/?=  1.034  for  a  Berkovich  indenter)  and  Er  is  the 
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TABLE  I.  Sputtering  parameters  for  AIN  thin  films. 


Target  (99.999%  purity) 

050  mm  aluminum 

dc  power  (W) 

200-275 

Substrate  temperature  (°C) 

450-550 

Target-substrate  distance  (mm) 

40 

Base  pressure  (torr) 

<10“6 

Sputtering  pressure  (mtorr) 

2-4 

Ar  gas  flow  rate  (SCCM) 

5-10 

N2  gas  flow  rate  (SCCM) 

12-20 

Substrate 


AIN  film 

Bottom 

electrode 


FIG.  1.  (Color  online)  The  schematic  four-layer  structure  of  the  thin  film 
composite  resonator. 


reduced  elastic  modulus,  which  accounts  for  the  fact  that 
elastic  deformation  occurs  in  both  the  sample  and  the  in¬ 
dented  Er  is  given  by 


1  -  vl 

E,.  = - 1  + 


l-vf 


(7) 


where  Er  and  Ei  are  the  elastic  modulus  of  the  thin  film 
sample  and  the  indenter,  and  Vj ■  and  u!  are  the  Poisson  ratios 
of  the  film  and  the  indenter,  respectively.  For  diamond,  Ej 
=  1141  GPa  and  e,  =  0.07. 

In  this  study,  fixed  probe  nanoindentation  method  was 
employed  to  characterize  the  nanomechanical  properties  of 
the  AIN  thin  films.  The  load-displacement  curves  were  ob¬ 
tained  at  different  load  conditions,  and  the  hardness  and  the 
elastic  modulus  of  the  thin  films  can  be  mapped  at  different 
locations  of  the  films.  Comparing  with  other  indentation 
methods  that  use  cantilever  indenters,  fixed  probe  nanoin¬ 
dentation  has  a  few  advantages.  First,  since  there  is  a  tilting 
angle  for  cantilever  indenter,  the  indentation  made  by  the 
indenter  is  usually  not  totally  perpendicular  to  the  sample 
surface.  However,  the  probe  for  the  fixed  probe  nanoindenter 
is  normal  to  the  sample  surface.  Second,  the  maximum  load 
for  a  fixed  nanoindenter  can  be  up  to  10  mN,  which  is  much 
larger  than  the  cantilever  indenter’s  maximum  load,  which  is 
normally  below  100  /cN.  It  should  be  pointed  that  in  the 
determination  of  elastic  modulus  and  hardness  of  supported 
thin  films,  complicated  substrate  influences  may  be  intro¬ 
duced  due  to  finite  sample  thickness  and/or  the  presence  of 
multiple  layers.  To  avoid  this  complication,  it  is  usually 
required  that  the  indentation  depth  is  less  than  one-tenth  of 
the  film  thickness. 


for  the  c-axis  oriented  AIN  film  growth  where  the  AIN  crys¬ 
tal  grains  grow  quasiepitaxial  on  Pt  (111)  substrate.25,26  Thus, 
before  sputtering  the  AIN  (002),  the  bottom  electrode  Pt/Ti 
was  grown  first  on  the  substrate.  Here,  Ti  thin  film  was  used 
as  the  adhesive  layer  between  Pt  thin  film  and  Si  substrate. 
The  Ti  thin  film  should  have  an  adequate  thickness  to  bond 
the  Pt  and  Si  together  but  not  affect  the  electromechanical 
properties  of  the  resonator.  Table  I  lists  the  processing  pa¬ 
rameters  used  in  the  deposition  process  for  AIN  film.  The 
deposition  time  of  the  AIN  thin  film  varied  from  3  to  6  h. 

B.  Fabrication  of  the  four-layer  composite 
resonator 

After  the  deposition  of  the  piezoelectric  AIN  layer,  hot 
H3PO4  acid  is  used  to  etch  the  AIN  film  to  expose  the  bottom 
electrode.  Then  the  top  electrode  A1  was  deposited  and  pat¬ 
terned  to  form  the  four-layer  composite  resonator  (Fig.  1). 

The  AIN  thin  film  composite  resonator  on  sapphire  sub¬ 
strate  was  used  to  characterize  the  effective  piezoelectric  co¬ 
efficient  G?33eff  of  AIN  film,  and  that  with  thin  silicon  sub¬ 
strate  (100  yum)  was  utilized  to  characterize  the  density  and 
stiffened  elastic  modulus  and  acoustic  wave  velocity  of  the 
AIN  thin  film  by  the  resonance  spectrum  method. 

C.  Measurement  of  effective  piezoelectric  coefficient 

^33eff 

The  setup  for  the  single  beam  laser  interferometer  is 
shown  in  Fig.  2.  The  laser  beam  from  the  He-Ne  laser-light 
source  is  split  into  two  equal  components  by  the  beam  split¬ 
ter.  One  of  the  beams  is  incident  on  the  reference  mirror  and 
the  other  on  the  sample  to  be  tested  where  both  lasers  travel 


III.  EXPERIMENTAL  PROCEDURES 
A.  Deposition  of  the  AIN  thin  film 

The  AIN  thin  films  were  deposited  by  using  dc  reactive 
magnetron  method  (ATC-1300  dc  and  rf  magnetron  sputter¬ 
ing  systems,  AJA  International,  MA)  from  an  A1  target  of 
high  purity  (99.999%)  in  the  high  purity  nitrogen  (99.999%) 
and  argon  (99.999%)  gas  mixture.  Before  the  sputtering,  the 
silicon  (100)  and  sapphire  (0001)  substrates  were  cleaned. 

The  (100)  direction  of  silicon  is  the  more  common  wafer 
orientation  for  microelectromechanical  devices.  Unfortu¬ 
nately,  the  c-axis  orientation  growth  of  the  AIN  film  on  the 
silicon  substrate  with  this  direction  is  very  difficult.  There 
are  also  a  large  lattice  mismatch  and  a  large  difference  in  the 
coefficient  of  thermal  expansion  between  the  AIN  film  and 
the  Si  (100)  substrate.  The  Pt  (111)  is  an  ideal  buffer  layer 


FIG.  2.  Experimental  setup  of  laser  interferometer  for  the  measurement  of 
piezoelectric  coefficient. 
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paths  of  equal  length.  There  is  a  reflected  beam  from  each 
surface  and  the  beams  are  then  combined  in  the  beam  splitter 
and  the  intensity  is  gathered  at  the  photodetector.  The  abso¬ 
lute  phase  of  each  beam  in  turn  depends  on  the  distance 
between  the  beam  splitter  and  the  respective  reflecting  sur¬ 
face.  If  the  surface  of  the  test  sample  moves  due  to  the  ex¬ 
ternal  voltage  applied,  the  absolute  phase  of  the  combined 
probe  beam  will  change.  By  detecting  the  electric  signal 
variation  induced  by  the  phase  change,  the  displacement  of 
the  sample  surface  can  be  obtained.  To  eliminate  possible 
electric  and  mechanical  noise  at  the  measurement  frequency 
(1  or  10  Hz),  a  lock-in  amplifier  was  used  to  detect  the  signal 
at  the  measurement  frequency;  thus  measurement  sensitivity 
can  be  improved. 

D.  Measurement  of  the  scattering  parameters  Sn 

When  an  electric  device  operates  at  the  low  frequency 
range,  the  dimensions  of  the  elements  in  the  circuit  are  much 
smaller  than  the  wavelength  and  thus  the  circuit  can  be 
treated  as  the  interconnection  of  lumped  elements.  The  volt¬ 
age  and  current  are  unique  at  any  point  in  the  circuit.  There¬ 
fore,  we  can  use  electric  impedance  and  admittance  param¬ 
eters  to  describe  the  circuit  network.  However,  when  the 
operating  frequency  is  up  to  the  radio  frequency  (rf)  or  mi¬ 
crowave  wave  range,  the  dimensions  of  the  elements  are 
larger  than  the  wavelength.  The  impedance  and  admittance 
are  not  well  defined  in  such  situation.  Normalized  wave  vari¬ 
ables  (a  and  b )  are  more  suitable  for  the  characterization  of 
the  circuit  at  microwave  or  higher  frequencies.  The  scatter¬ 
ing  parameters  or  S  parameters  are  used  to  express  the  rela¬ 
tionship  of  the  normalized  wave  variables.  Unlike  other  pa¬ 
rameters,  the  S  parameters  relate  to  the  traveling  waves  that 
are  scattered  or  reflected  when  a  network  is  inserted  into  a 
transmission  line  with  certain  characteristic  impedance  Z0. 
For  a  two-port  network,  there  are  totally  four  scattering  pa¬ 
rameters  Sjj  ( i,j=  1,2).  Each  S  parameter  has  two  subscripts, 
i  and  j,  where  i  refers  to  the  port  number  at  which  power 
emerges  and  j  refers  to  the  port  number  at  which  power 
enters.  Thus  ,S'| ,  and  S22  are  the  reflection  coefficients  and  S , 2 
and  S21  are  the  transmission  coefficients.  Specifically,  reflec¬ 
tion  scattering  parameter  Sn  is  defined  as 


2  power  reflected  from  the  network  input 

(8) 

S 1 1  —  _  _ 

power  incident  on  the  network  input 

It  reflects  the  impedance  Zin  and  admittance 
resonator  through  the  following  equations: 

Tin  of  the 

1  +  <Sn 

Zm-Zo,  ", 

l  - 

(9) 

1  1-Sn 

7in~z0i  +  V 

(10) 

where  Z0=50  fl  is  the  standard  normalized  impedance  for 
the  network  analyzer.  Measurement  of  the  S  parameter  is 
usually  performed  using  a  vector  network  analyzer  (VNA). 
In  the  measurement  of  the  S  parameter,  a  careful  calibration 
of  the  VNA  is  very  important  to  get  accurate  results  for  the 
devices  under  test.  VNA  calibration  is  the  process  of  mea- 


FIG.  3.  Schematic  of  a  fixed  probe  nanoindenter  with  in  situ  imaging 
capability. 

suring  devices  with  known  characteristics  and  using  these 
measurements  to  establish  the  measurement  reference  planes. 
Calibration  also  corrects  for  the  imperfections  of  the  mea¬ 
surement  system.  These  imperfections  include  not  only  the 
nonideal  nature  of  cables  and  probes  but  also  the  internal 
characteristics  of  the  VNA  itself.  A  VNA  uses  multiple  stan¬ 
dards,  typically  open  circuit,  short  circuit,  loads,  and  through 
connections,  to  perform  the  calibration. 

The  51! !  of  AIN  thin  film  composite  resonator  was  mea¬ 
sured  by  an  Agilent  network  analyzer  (Agilent  85 10C  vector 
network  analyzer,  45  MHz  to  50  GHz,  Agilent  Technologies, 
Palo  Alto,  CA).  Calibrations  were  first  performed  using 
open,  short,  and  standard  (50  (1).  In  the  measurement,  the 
four-layer  thin  film  composite  resonator  was  placed  on  the 
probe  station  (Summit  9101  manual  probe  station.  Cascade 
MicroTech,  OR),  which  contacted  the  sample  bottom  surface 
only  at  the  edge  so  that  the  resonator  could  vibrate  with  free 
boundary  condition  in  the  thickness  direction.  A  proper  SG 
probe,  which  has  two  pins  on  the  tip,  was  used  to  contact  the 
electrodes  of  the  composite  resonator  electrically.  The  other 
end  of  the  probe  was  connected  to  the  network  analyzer 
through  a  cable. 

The  reflection  scattering  parameter  5)  |  was  measured 
from  400  to  2500  MHz.  The  entire  frequency  span  of 
2100  MHz  was  divided  into  tens  of  narrow  subspans  to 
achieve  higher  resolution  in  the  resonance  spectrum  mea¬ 
surement. 


E.  Nanomechanical  test  of  AIN  thin  films 

The  nanoindentation  tests  were  performed  in  a  Hysitron 
TriboScope  (Minneapolis,  MN)  attached  to  a  Quesant  (Ago- 
ura  Hills,  CA)  atomic  force  microscope  (AFM).  Figure  3  is  a 
schematic  of  the  working  principle  for  the  Hysitron  nanoin¬ 
denter.  In  situ  imaging  capability  allows  the  indenter  probe 
to  be  positioned  to  any  desired  material  feature  within 
10  nm.  A  Berkovich  diamond  indenter  with  a  nominal  tip 
radius  of  40  nm  is  used.  AIN  films  with  thicknesses  of  0.6, 
1.2,  and  1.5  /im  deposited  on  sapphire  substrate  were  used 
for  this  study.  The  maximum  loads  are  varied  from 
500  to  4000  juN.  The  results  of  the  loads  and  displacements 
during  the  test  offer  us  a  lot  of  information  including  the 
reduced  elastic  modulus  and  hardness. 
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FIG.  4.  X-ray  diffraction  patterns  of  AIN  films  deposited  at  550  °C  on 
different  substrates:  (a)  sapphire  (0001)  and  (b)  silicon  (100). 

IV.  RESULTS  AND  DISCUSSIONS 

A.  Surface  morphology  and  crystal  structure  of  AIN 
film 


After  the  deposition  of  the  bottom  electrode  and  AIN 
thin  films,  the  thin  him  crystalline  structures  were  character¬ 
ized  by  x-ray  diffraction  (XRD)  and  scanning  electron  mi¬ 
croscopy  (SEM).  Figure  4  shows  the  XRD  patterns  of  the 
AIN  thin  films  deposited  on  different  substrates  with  Pt/Ti 
bottom  electrode  layers.  Only  AIN  (002)  and  Pt  (111)  peaks 
are  observed,  indicating  that  highly  c-axis  oriented  AIN  films 
have  been  achieved.  Figures  5(a)  and  5(b)  show  the  SEM 
images  of  the  cross  section  of  the  AIN  thin  films  deposited 
on  the  sapphire  (0001)  and  silicon  (100)  substrates,  respec¬ 
tively.  The  SEM  images  reveal  that  the  AIN  films  have  dense 
column  grains  that  are  perpendicular  to  the  surface  of  the 
substrates.  Figure  6  shows  the  surface  morphology  of  the 
AIN  thin  film  on  silicon  substrate.  The  results  indicate  that 
under  the  processing  conditions,  AIN  thin  films  deposited  are 
smooth,  crack-free,  and  with  a  good  uniformity. 

B.  Piezoelectric  coefficient  d33  of  AIN  film 

With  quasistatic  driving  voltage  at  1  Hz  and  peak  to 
peak  voltage  from  0  to  5  V,  the  voltage  induced  dynamic 


FIG.  5.  SEM  micrographs  of  AIN  films  deposited  at  550  “Con  (a)  sapphire 
(0001)  substrate  and  (b)  silicon  (100)  substrate. 

displacement  of  the  piezoelectric  AIN  thin  film  was  mea¬ 
sured  using  the  single  beam  laser  interferometer.  An  effective 
piezoelectric  strain  coefficient  of  d33eS=  2.9  pm/V  was  de¬ 
tected  for  AIN  film  in  the  multilayer  composite  resonator 
structure  (electrode/AIN  thin  film/electrode/sapphire).  Since 
the  piezoelectric  coefficient  is  measured  under  the  fact  that 
the  AIN  film  is  always  clamped  to  the  sapphire  substrate,  the 
displacement  and  driving  voltage  ratio  A d/V  does  not  reflect 
the  true  piezoelectric  coefficient  d33  of  the  AIN  film  under 
free  boundary  condition.  In  the  converse  piezoelectric  effect, 
the  effective  piezoelectric  coefficient  d33 eff  of  thin  film  ma¬ 
terials  can  be  related  to  the  true  piezoelectric  coefficient  by 
the  following  relationship” 

^33eff  =  ^33  -  2^3lH3/(-s’n  +  st 2)  »  (1 1) 

where  .v  1 2,  and  s13  are  the  mechanical  compliances  of  the 
piezoelectric  film  and  d3\  is  the  transverse  piezoelectric  co¬ 
efficient.  Since  d33~-2d3\,  we  have 


FIG.  6.  Surface  morphology  of  the  AIN  thin  film  on  silicon  substrate. 


Downloaded  21  Sep  2010  to  150.212.9.124.  Redistribution  subject  to  AIP  license  or  copyright;  see  http://jap.aip.org/about/rights_and_permissions 


084103-6  Chen,  Qin,  and  Wang 


J.  Appl.  Phys.  101,  084103  (2007) 


TABLE  II.  Properties  of  electrode  and  substrate  materials  used  in  the  cal¬ 
culation. 


Electrode 

Density  (kg/m3) 

Acoustic  velocity  (m/s) 

Mechanical  Q 

A1 

2  700 

6418 

2000 

Pt 

21  450 

2823 

500 

Si 

2  320 

8270 

3000 

t/33  -  ^33eff(ifi  +  ^l{s\x  +  sn  +  s13)-  (12) 

Substituting  the  values  of  s12,  and  s13  into  Eq.  (12),  the 
true  piezoelectric  coefficient  <y33 =4.17  pm/V  can  be  calcu¬ 
lated,  which  is  in  fairly  good  agreement  with  the  reported 
values  in  literatures.28  31 

C.  Density  and  acoustic  velocity  of  AIN  film 

Based  on  the  one-dimensional  input  electric  impedance 
[Eq.  (1)],  the  electric  impedance  spectra  for  a  four-layer  thin 
film  composite  resonator  and  a  three-layer  resonator  can  be 
simulated.  Listed  in  Tables  II  and  III  are  the  material  prop¬ 
erties  used  in  the  calculation.  Figures  7(a)  and  7(b)  show  the 
typical  simulation  results  of  the  electric  impedance  spectra 
for  composite  resonators,  in  which  material  properties  of 
each  layer  are  plugged  in  the  equation  for  calculation.  It 
should  be  mentioned  that  a  small  imaginary  part  of  the 
acoustic  velocity  for  each  material  has  been  introduced  in  the 
calculation,  which  can  stand  for  the  mechanical  quality  fac¬ 
tor  of  the  material.  The  mechanical  quality  factor  Q  in  the 
table  represents  half  of  the  real  part  to  imaginary  part  ratio  of 
the  acoustic  velocity.13,14  The  impedance  of  the  composite 
resonator  follows  a  hyperbolic  decrease  that  is  the  presenta¬ 
tion  of  the  static  capacitance  C0  of  the  piezoelectric  layer. 
There  are  multiple  sharp  peaks  corresponding  to  the  reso¬ 
nance  “modes.”  At  a  resonance  mode,  a  standing  wave  is 
established  in  the  thickness  direction  of  the  composite  reso¬ 
nator.  Since  the  substrate  is  much  thicker  than  the  piezoelec¬ 
tric  layer,  the  mode  frequencies  are  mainly  determined  by  the 
substrate  and  modulated  by  the  piezoelectric  layer.  There¬ 
fore,  both  the  piezoelectric  layer  and  the  substrate  will  affect 
the  distribution  of  multiple  resonance  modes  of  the  resonator. 
By  analyzing  the  space  of  the  adjacent  resonance  frequency 
at  special  regions,  for  instance,  the  first  normal  and  transition 
regions,  we  can  derive  the  material  properties  of  the  piezo¬ 
electric  layer. 

The  experimental  results  for  the  .Sj ,  parameter  and  the 
converted  impedance  spectrum  with  three  consecutive  peaks 
for  a  four-layer  composite  resonator  are  shown  in  Figs.  8(a) 
and  8(b).  By  measuring  the  parallel  resonance  frequency  of 
these  modes,  we  can  obtain  the  space  of  two  adjacent  parallel 
resonance  frequencies  A fP  over  the  measurement  frequency 


0.5  1.0  1.5  2.0  2.5 

Frequency  (GHz) 


FIG.  7.  Simulated  electric  impedance  spectra  of  (a)  a  four-layer  AIN  thin 
film  composite  resonator  and  (b)  a  three-layer  AIN  thin  film  composite 
resonator. 

range.  Figure  9  shows  the  experiment  results  and  the  simu¬ 
lation  results  of  A fP.  Although  the  experimental  data  are  a 
little  dispersive,  the  periodic  shape  of  the  curve  is  in  good 
agreement  with  the  simulation  results.  By  averaging  the  mea¬ 
surement  data  around  the  center  of  the  first  normal  and  tran¬ 
sition  regions,  we  can  obtain  values  of  A fr  and  A fN  of  40.81 
and  39.86  MHz,  respectively. 

Using  Eqs.  (3)  and  (4),  the  density  and  acoustic  velocity 
of  the  thin  film  are  attained  as  3187.3  kg/m3  and  10631  m / s. 
From  Eq.  (5),  the  stiffened  elastic  constant  of  the  AIN  thin 
film  is  341.2  GPa. 

D.  Elastic  modulus  and  hardness  of  the  c-axis 
oriented  AIN  film 

Using  the  nanoindentation,  the  load-displacement  curves 
were  obtained.  Figure  10  shows  the  nanoindentation  load- 
displacement  curves  for  the  AIN  thin  films  with  thicknesses 
of  0.6  and  1.5  /im  that  are  deposited  on  sapphire  substrates. 
Maximum  loads  Pmax  for  the  indentation  are  varied  from 


TABLE  III.  Properties  of  the  piezoelectric  thin  film  used  in  the  simulation  in  the  calculation. 


Piezoelectric 

film 

Density 

(kg/m3) 

Acoustic 

velocity 

(m/s) 

if 

(%) 

Thickness 

(H 

Area 

(m2) 

Q 

sr 

AIN 

3300 

11  000 

6.25 

1.9 

2.5  X  10“9 

500 

8.5 
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FIG.  8.  (a)  Scattering  parameter  Su  measured  by  the  network  analyzer  and 
(b)  impedance  spectrum  for  a  four-layer  AIN  composite  resonator. 

500  to  4000  /zN  during  the  tests.  All  the  load-displacement 
curves  are  smooth  and  regular  without  discontinuities,  which 
indicate  that  no  cracking  occurred  in  the  penetration  process. 
From  the  load-displacement  curves  for  various  applied  maxi¬ 
mum  loads,  the  reduced  elastic  modulus  Er  and  hardness  H 
can  be  obtained. 

Figures  11(a)  and  11(b)  show  the  reduced  elastic  modu¬ 
lus  Er  and  hardness  of  AIN  thin  films  with  different  thick- 


FIG.  9.  (Color  online)  Space  of  the  parallel  resonance  frequency  as  a  func¬ 
tion  of  parallel  resonance  frequency. 


FIG.  10.  Nanoindentation  load-displacement  curves  on  AIN  thin  films  with 
different  thicknesses  and  under  different  load  conditions:  (a)  thickness  of  the 
AIN  thin  film  tf-  0.6  /rm;  (b)  thickness  of  the  AIN  thin  film  tf— 1.5  /cm. 

nesses  under  different  maximum  loads  at  room  temperature. 
No  significant  difference  is  observed  for  values  of  the  elastic 
modulus  and  hardness  of  AIN  thin  films  with  different  thick¬ 
nesses.  For  all  samples  under  test,  the  thickness  of  AIN  film 
is  quite  thick  as  compared  with  the  indentation  depth;  the 
effect  of  the  substrate  is  thus  insignificant.  The  variation  of 
the  properties  could  be  due  to  the  surface  conditions  such  as 
roughness  at  different  indentation  locations  on  the  films.  In 
addition,  indentation  at  or  close  to  the  grain  boundaries  may 
lead  to  lower  values  of  elastic  modulus  and  hardness  as  can 
be  expected.  The  average  reduced  modulus  of  all  thin  film 
samples  is  about  250  GPa  and  the  average  hardness  of  the 
thin  films  is  about  18  GPa. 

The  elastic  modulus  is  related  to  the  reduced  elastic 
modulus  by  Eq.  (7).  Substituting  the  Poisson  ratio  of  0.2  for 
AIN,  "  the  elastic  modulus  and  Poisson  ratio  for  diamond 
indenter,  and  with  the  average  reduced  elastic  modulus  into 
the  equation,  the  average  elastic  modulus  of  the  AIN  thin 
films  can  be  estimated  to  be  about  307  GPa. 

V.  SUMMARY 

Smooth  piezoelectric  AIN  thin  films  were  grown  on  the 
sapphire  and  silicon  substrates  by  the  dc  reactive  magnetron 
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Load  <nN> 

FIG.  11.  (a)  Reduced  elastic  modulus  and  (b)  hardness  of  the  AIN  thin  films 
with  different  maximum  loads  at  room  temperature. 

sputtering.  SEM  and  XRD  characterization  results  show  that 
the  films  show  highly  c-axis  orientation.  Single  beam  laser 
interferometer  system  was  used  to  measure  the  effective  pi¬ 
ezoelectric  coefficient  of  the  AIN  thin  film  deposited  on  sap¬ 
phire  substrate.  The  effective  and  true  piezoelectric  strain 
coefficients  <^33eff  and  c/33  °f  the  AIN  thin  films  are  2.9  and 
4.17  pm/V,  respectively.  Thin  film  composite  resonators 
have  also  been  fabricated  on  silicon  substrate  for  the  charac¬ 
terization  of  electromechanical  properties.  The  .S’-paramcter 
(5„)  spectrum  of  the  composite  resonator  was  first  measured 
at  a  wide  frequency  range  and  converted  into  electric  imped¬ 
ance  spectrum.  Using  the  space  of  two  adjacent  parallel  reso¬ 
nance  frequencies  at  the  first  normal  and  transition  regions  in 
the  electric  impedance,  the  density  and  acoustic  velocity  of 
the  AIN  thin  film  were  calculated  to  have  values  of 
3187.3  kg/m3  and  10631  m/s,  respectively.  Consequently 
the  stiffened  elastic  modulus  cD  of  the  AIN  film  is 


341.2  GPa.  The  results  obtained  from  the  experiments  are  in 
good  agreement  with  those  for  bulk  AIN  crystals. 
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